SECURITY  CLASSIFICATION  Or  THIS  RAGE  Dmtm  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  OOVT  ACCESSION  NO. 

SIO  REFERENCE  85-33 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (tnd  Submit) 

DEEP  SEAFLOOR  ACOUSTIC  BACKSCATTERING 
MEASUREMENTS  USING  SEA  BEAM 

S.  TYPE  OF  REPORT  ft  PERIOD  COVEREO 

Summary 

t.  PERFORMING  ORG.  REPORT  NUMBER 

MPL-U- 54/85 

7.  AUTHOR*"*; 

Christian  de  Moustier 

».  CONTRACT  OR  GRANT  NUMBER*-*; 

N00014-79-C-0472 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  California,  San  Diego,  Marine 
Physical  Laboratory,  Scripps  Institution  of 
Oceanography,  San  Diego,  CA  92152 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 

AREA  &  WORK  UNIT  NUMBERS 

1  1.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research,  Department  of  the 

Navy,  800  North  Quincy  Street,  Arlington 

Virginia  22217 

12.  REPORT  DATE 

December  1985 

13.  NUMBER  OF  PAGES 

142 

14.  MONITORING  AGENCY  NAME  *  ADDRESS^//  different  from  Controlling  Office) 

15.  SECURITY  CLASS,  (ol  Ihlt  report) 

UNCLASSIFIED 

15a.  DECLASSI  FI  CATION/ DOWNGRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (o  1  thle  Report) 

Document  cleared  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (o l  the  ebetrect  entered  in  Block  20,  it  different  from  Report) 

18.  supplementary  NOTES 

19.  KEY  WOROS  (Continue  on  teeetee  eide  It  neceeeery  end  identify  by  block  number) 

multibeam  echo  sounders,  sea  beam,  bathymetric  surveys,  acoustic 
boundary  mapping,  acoustic  backscatter 

20.  ABSTRACT  ( Contlnu e  on  reveree  elde  1/  n«ctf««ry  and  /danUfy  by  bloc*  numbrnr) 


Multibeam  echo-sounders  such  as  Sea  Beam  allow  investigators  to  carry 
out  detailed  bathymetric  surveys  of  large  areas  of  the  seafloor.  However, 
bathymetry  only  reveals  the  shape  of  seafloor  features  to  the  resolution 
o  the  sounding  system,  and  in  order  to  make  geological  interpretations  one 
needs  to  characterize  the  nature  of  the  seafloor  surveyed. 

Because  bottom  roughness  and/or  variations  in  bottom  substrate  cause 
fluctuations  in  the  backscattered  acoustic  signal  received  by  an  echo-sounder 
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such  characteristics  can  be  inferred  in  part  by  analyzing  the  structure 
and  the  variations  of  this  signal  over  several  transmission  cycles. 

The  approach  taken  has  been  to  record  digitally  the  detected 
echo  envelopes  of  Sea  Beam's  16  narrow  beams  over  a  variety  of  seafloor 
environments,  and  process  these  data  to  determine  whether  the  acoustics 
held  enough  information  to  differentiate  between  bottom  types. 

Significant  results  derived  from  these  acoustic  data  concern  (1) 
the  Sea  Beam  system's  performance  (2)  its  potential  for  mapping 
acoustic  boundaries  and  (3)  the  display  of  the  echoes  received  in  a 
side  looking  sonar-like  picture. 

The  system  was  found  to  work  well  under  most  circumstances  as  a 
contour  mapping  tool,  but  it  occasionally  suffers  echo  processing 
malfunctions  producing  artifacts  in  the  contoured  bathymetric  output 
which  can  induce  the  unwary  investigator  to  make  geological  interpre¬ 
tation  errors  [de  Moustier  and  Kleinroce  (1985)  J.G.R.  in  press]. 

Acoustic  boundary  mapping  has  been  successfully  carried  out  over  a 
manganese  nodule  mining  site  in  the  Northeastern  Tropical  Pacific  where 
a  first  order  assessment  of  nodule  coverage  validated  with  bottom 
photographs  proved  the  techniques  feasible  [de  Moustier  (1985) 

Geophysics,  V.  50,  p.  989-1001].  Comparison  of  results  from  environments 
as  diverse  as  a  lava  sheet  flow  on  the  crest  of  the  East  Pacific  Rise, 
a  sedimentary  basin  offshore  Southern  California,  and  a  manganese 
nodule  field  show  marked  differences  in  the  overall  backscattered 
acoustic  energy  level  as  well  as  in  the  inferred  roughness  character 
between  these  areas,  [de  Moustier  J.A.S.A.,  submitted].  Although 
these  simple  characteristics  are  insufficient  to  determine  the  nature 
of  the  seafloor  surveyed,  they  represent  additional  clues  to  help 
understand  the  geological  processes  under  investigation. 

Further  insight  into  these  processes  is  afforded  by  the 
qualitative  measure  of  acoustic  backscatter  inherent  to  a  gray  scale 
side  looking  sonar  display.  Such  displays  have  been  created  with  Sea 
Beam  acoustic  data.  They  contain  a  wealth  of  information  not  available 
in  the  bathymetric  contours  and  are  expected  to  become  a  very  valuable 
addition  to  the  Sea  Beam  system. 
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ABSTRACT  OF  THE  DISSERTATION 


Deep  seafloor  acoustic  backscattering 
measurements  using  Sea  Beam 

by 

Christian  de  Moustier 
Doctor  of  Philosophy  in  Oceanography 
University  of  California,  San  Diego,  1985 
Professor  Fred  N.  Spiess,  Co-Chairperson 
Associate  Professor  Robert  C.  Tyce,  Co-Chairperson 

Multibeam  echo-sounders  such  as  Sea  Beam  allow  investigators  to  carry  out 
detailed  bathymetric  surveys  of  large  areas  of  the  seafloor.  However,  bathymetry  only 
reveals  the  shape  of  seafloor  features  to  the  resolution  of  the  sounding  system,  and  in 
order  to  make  geological  interpretations  one  needs  to  characterize  the  nature  of  the 
seafloor  surveyed. 

Because  bottom  roughness  and/or  variations  in  bottom  substrate  cause 
fluctuations  in  the  backscattered  acoustic  signal  received  by  an  echo-sounder,  such 
characteristics  can  be  inferred  in  part  by  analyzing  the  structure  and  the  variations  of 
this  signal  over  several  transmission  cycles. 

The  approach  taken  has  been  to  record  digitally  the  detected  echo  envelopes 
of  Sea  Beam’s  16  narrow  beams  over  a  variety  of  seafloor  environments,  and  process 
these  data  to  determine  whether  the  acoustics  held  enough  information  to  differentiate 
between  bottom  types. 

Significant  results  derived  from  these  acoustic  data  concern  (1)  the  Sea  Beam 
system’s  performance  (2)  its  potential  for  mapping  acoustic  boundaries  and  (3)  the 
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This  dissertation  presents  the  results  of  an  experiment  to  measure  seafloor 
acoustic  backscatter  with  the  multibeam  echo-sounder  Sea  Beam.  With  this  system, 
sound  transmission  and  reception  are  achieved  by  two  adjacent  acoustic  arrays  so 
that  the  echo  signals  received  correspond  to  the  portion  of  the  acoustic  energy  which 
is  scattered  back  towards  the  sound  source  in  a  direction  opposite  that  of  the  incident 
sound.  Interest  in  such  measurements  stems  from  the  fact  that  changes  in  the  compo¬ 
sition  or  the  roughness  structure  of  the  seafloor  cause  variations  in  the  acoustic  signals 
scattered  by  the  seafloor. 

Sea  Beam  processes  the  bottom  echo  signals  it  receives  to  determine  the 
bathymetric  data  necessary  for  contour  mapping  and  then  discards  these  signals.  A 
parallel  acoustic  data  acquisition  system  is  therefore  required  to  preserve  them  for 
further  processing  and  analysis.  The  approach  taken  here  consists  in  recording  digi¬ 
tally  the  bottom  echo  signals  received  by  Sea  Beam  over  a  variety  of  seafloor  environ¬ 
ments.  As  a  first  step  in  the  processing  and  analysis  of  these  data,  the  recorded 
acoustic  signals  are  used  to  assess  the  performance  of  the  Sea  Beam  system.  The  next 
step  concerns  use  of  the  amplitude  variations  in  the  recorded  signals  to  produce  sides- 
can  images  of  the  seafloor.  Further  analysis  of  the  structure  and  variations  of  these 
echo  signals  over  several  transmission  cycles  is  used  to  relate  information  contained  in 
the  acoustics  to  seafloor  properties.  In  the  scope  of  this  dissertation,  this  analysis  is 
limited  to  areas  of  the  seafloor  which  are  nearly  flat.  However,  acoustics  alone  do  not 
give  a  unique  answer  about  the  nature  of  the  seafloor,  and  independent  measurements 
(e.g.  bottom  photographs,  core  samples,  etc.)  provide  the  necessary  ground  truth. 

Following  the  analysis  sequence  outlined  above,  this  dissertation  is  divided 
into  four  chapters  and  an  appendix.  Each  chapter  is  self-contained  and  corresponds  to 
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a  paper  that  has  been  published  by,  or  submitted  to  a  refereed  scientific  journal.  As  a 
result,  some  redundancy  exits  between  chapters  whenever  the  Sea  Beam  system  is 
described. 

The  first  chapter  reviews  the  operational  characteristics  of  the  Sea  Beam  sys¬ 
tem  and  uses  the  recorded  acoustic  data  to  take  an  in-depth  look  at  the  system’s  per¬ 
formance  as  a  mapping  tool.  Although  the  system  works  well  under  most  cir¬ 
cumstances,  several  bathymetric  artifacts  have  been  observed  in  the  contoured  out¬ 
puts.  These  artifacts  are  analyzed  and  their  causes  are  identified  by  examining  the 
character  of  the  bottom  echoes  received  and  processed  by  the  system.  This  analysis 
leads  to  recommendations  for  improving  bathymetric  data  quality.  This  chapter 
represents  a  paper  co-authored  by  Martin  C.  Kleinrock  who  contributed  examples  of 
bathymetric  artifacts  found  in  the  data  contoured  by  Sea  Beam,  and  who  assessed 
their  geological  implications.  This  paper  has  been  accepted  by  the  Journal  of  Geophy¬ 
sical  Research. 

In  the  second  chapter,  the  variations  in  amplitude  of  the  echo  signals 
received  by  Sea  Beam  are  used  to  produce  a  sidescan  image  of  the  seafloor  surveyed. 
This  is  a  new  application  of  the  Sea  Beam  system  which  is  expected  to  enhance  its 
usefulness  during  geological  investigations  of  the  seafloor,  by  complementing  the  stan¬ 
dard  bathymetric  contouring  with  a  qualitative  measure  of  bottom  acoustic  back- 
scatter  which  is  related  to  the  texture  of  the  seafloor.  This  chapter  corresponds  to  a 
short  invited  paper  submitted  to  the  American  Society  of  Mechanical  Engineers  for  its 
symposium  on  Ocean  Engineering.  As  such,  it  also  reviews  the  approaches  currently 
taken  by  the  oceanographic  community  to  measure  acoustic  backscattering  from  the 
deep  seafloor. 

The  third  chapter  is  a  first  attempt  at  relating  the  character  of  the  acoustic 
backscatter  to  the  geological  processes  that  cause  it.  The  variations  of  the  peak 
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amplitude  of  Sea  Beam’s  near-specular  return  along  the  ship’s  track  are  compared 
with  the  density  of  coverage  of  manganese  nodules  at  a  mining  site  in  the 
Northeastern  Tropical  Pacific.  A  98%  qualitative  agreement  is  found  between  man¬ 
ganese  nodule  coverage  inferences  obtained  with  Sea  Beam’s  acoustics  and  with  bot¬ 
tom  photographs  taken  by  the  Deep  Tow  instrument  package  in  the  same  area.  The 
paper  corresponding  to  this  chapter  is  published  in  Geophysics  [V.50,  989-1001,  1985]. 

The  fourth  chapter  expands  the  analysis  of  Sea  Beam’s  acoustics  to  nonspec- 
ular  beams.  Measurements  from  three  geologically  different  areas  are  compared. 
Acoustic  backscattering  boundaries  are  mapped  both  along  and  across  the  ship’s  track 
by  following  trend  changes  in  the  total  energy  of  the  acoustic  returns  over  many 
transmission  cycles.  A  partial  angular-dependence  function  of  backscattering  shows 
marked  differences  in  the  overall  acoustic  energy  level  recorded  at  each  site.  Applica¬ 
tion  of  Rician  statistics  to  the  distribution  of  peak  amplitude  in  the  near-specular 
direction  gives  clues  to  the  roughness  character  of  each  area,  and  provides  a  means  to 
differentiate  them.  This  chapter  represents  a  paper  which  has  been  submitted  for  pub¬ 
lication  to  the  Journal  of  the  Acoustical  Society  of  America. 

Because  of  editorial  constraints  imposed  by  publishers,  the  method  used  to 
record  Sea  Beam’s  acoustic  data  is  mentioned  only  briefly  in  these  chapters.  For  this 
reason,  a  brief  technical  description  of  the  Sea  Beam  system  and  a  detailed  descrip¬ 
tion  of  the  method  used  to  record  the  acoustic  signals  it  receives  are  included  in  the 
appendix. 

Possible  extensions  of  this  work  include  application  of  the  Sea  Beam  results 
to  other  multibeam  echo-sounders  and  refinement  of  the  acoustic  data  processing  tech¬ 
niques  to  account  for  bottom  slopes  effects  in  areas  of  large  relief. 


Chapter  I 


BATHYMETRIC  ARTIFACTS  IN  SEA  BEAM  DATA 
HOW  TO  RECOGNIZE  THEM,  WHAT  CAUSES  THEM. 

1.1.  ABSTRACT 

Sea  Beam  multibeam  bathymetric  data  have  greatly  advanced  understanding 
of  the  deep  seafloor.  However,  several  types  of  bathymetric  artifacts  have  been 
identified  in  Sea  Beam’s  contoured  output.  Surveys  with  many  overlapping  swaths 
and  digital  recording  on  magnetic  tape  of  Sea  Beam’s  sixteen  acoustic  returns  made  it 
possible  to  evaluate  actual  system  performance.  The  artifacts  are  not  due  to  the 
contouring  algorithm  used.  Rather,  they  result  from  errors  in  echo  detection  and 
processing.  These  errors  are  due  to  internal  factors  such  as  sidelobe  interference, 
bottom-tracking  gate  malfunctions,  or  external  interference  from  other  sound  sources 
(e.g.  3.5  kHz  echo-sounders  or  seismic  sound  sources).  Although  many  artifacts  are 
obviously  spurious  and  would  be  disregarded,  some  (particularly  the  ‘omega’  effects 
described  in  this  paper)  are  more  subtle  and  could  mislead  the  unwary  observer. 
Artifacts  observed  could  be  mistaken  for  volcanic  constructs,  abyssal  hill  trends, 
hydrothermal  mounds,  slump  blocks  or  channels,  and  could  seriously  affect  volcanic, 
tectonic  or  sedimentological  interpretations.  Misinterpretation  of  these  artifacts  may 
result  in  positioning  errors  when  seafloor  bathymetry  is  used  to  navigate  the  ship. 
Considering  these  possible  geological  misinterpretations,  a  clear  understanding  of  the 
Sea  Beam  system’s  capabilities  and  limitations  is  deemed  essential. 
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1.2.  INTRODUCTION 

The  Sea  Beam  bathymetric  survey  system  is  a  multibeam  echo-sounder 
developed  by  the  General  Instrument  Corporation  to  produce  near-real-time  high 
resolution  contoured  swath  charts  of  the  seafloor  down  to  maximum  ocean  depth  (11 
km).  Since  1977  when  the  first  system  became  operational  aboard  the  French  R/V 
Jean  Charcot ,  nine  other  systems  have  been  installed  aboard  research  vessels  from  the 
United  States,  Germany,  Japan  and  Australia. 

Sea  Beam  systems  have  proven  extremely  useful  in  the  study  of  the 
geomorphology  of  the  ocean  floor,  and  have  made  possible  striking  discoveries  of 
features  which  would  not  have  been  detected  with  conventional  single  point  depth 
sounders  |e.g.  Macdonald  and  Fox,  1983;  Lonsdale,  1983].  However,  after  three  years 
of  experience  with  the  system  installed  aboard  the  R/V  Thomas  Washington  of  the 
Scripps  Institution  of  Oceanography  (SIO),  we  have  discovered  a  number  of  artifacts 
in  Sea  Beam’s  contoured  output.  Their  artificial  nature  has  been  demonstrated  by 
comparing  overlapping  Sea  Beam  swaths  and  by  analyzing  digitized  raw  acoustic 
data.  During  four  cruises  aboard  the  R/V  Thomas  Washington  we  used  a  data 
acquisition  system  developed  by  SIO’s  Marine  Physical  Laboratory  (MPL)  to  record 
digitally  the  acoustic  returns  from  Sea  Beams’  sixteen  preformed  beams  on  magnetic 
tape.  This  data  set  has  enabled  us  to  determine  the  causes  of  these  artifacts.  They 
do  not  stem  from  the  vagaries  of  the  contouring  algorithm  used,  rather  they  are  the 
result  of  errors  in  echo  detection  and  processing.  In  our  experience,  such  detection 
errors  are  mostly  related  to  some  characteristics  of  the  ocean  bottom  (e.g.  type  of 
substrate  or  sudden  change  in  slope)  or  to  interference  from  other  sound  sources 
running  in  parallel  with  Sea  Beam  (mostly  subbottom  profilers:  3.5  kHz  echo-sounder 
and  seismic  sources). 

This  paper  describes  several  artifacts  discovered  in  Sea  Beam  data  and 
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discusses  the  associated  possible  geological  misinterpretations.  We  suggest  a  number 
of  solutions  to  improve  data  quality.  We  also  consider  the  existence  of  related 
artifacts  in  similar  multibeam  echo-sounders  (e.g.  the  Sonar  Array  Sounding  System 
(SASS)  [Glenn,  1970]). 

1.3.  BACKGROUND:  THE  SEA  BEAM  SYSTEM 

Before  going  into  a  detailed  explanation  of  the  problems  found  in  the  SIO 
system,  we  briefly  review  Sea  Beam’s  general  framework  for  the  reader  unfamiliar 
with  the  system.  Further  discussion  of  the  system  are  found  in  Renard  and  Allenou 
1979i  and  in  Farr  1980  .  Because  a  clear  understanding  of  Sea  Beam’s  acoustic 
geometry  and  echo  processing  methods  is  a  prerequisite  to  analyze  its  bathymetric 
output,  we  have  included  in  the  appendices  relevant  technical  information  not 
available  in  the  literature. 

As  illustrated  by  the  simplified  block  diagram  in  Figure  1.1,  the  Sea  Beam 
system  uses  a  multibeam  narrow  beam  echo-sounder  and  an  echo  processor  (EP)  to 
generate  in  near-real  time,  contour  maps  of  the  ocean  floor.  A  20-element  projector 
array  mounted  along  the  ship's  keel  sends  out  a  7  millisecond  pulse  of  12.158  kHz  at 
intervals  that  are  integral  multiples  of  one  second.  The  transmission  period  is  usually 
determined  by  an  analog  graphic  recorder.  The  receiving  unit  lies  athwartships  and 
consists  of  40  line  hydrophone  arrays  whose  long  axes  are  oriented  fore-aft.  The 
resulting  transmit/receive  geometry  is  illustrated  in  Figure  1.2  In  this  figure,  vertical 
cross  sections  of  theoretical  beam  patterns  are  shown  for  both  the  projector  (Fig.  1.2a) 
and  the  receiver  (Fig.  1.2b)  arrays.  The  transmit  beam  pattern  spans  54  degrees 
athwartships  by  2  2/3  degrees  in  the  fore-aft  direction.  It  is  pitch  stabilized  within  a 
range  of  ±  10  degrees  of  pitch.  There  is  no  pitch  compensation  for  the  receive  beam 
pattern,  instead  it  spans  20  degrees  in  the  fore-aft  direction  to  accommodate  pitch 
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Figure  1.1  Block  diagram  of  the  Sea  Beam  system  showing  the  Narrow  Beam  Echo  Sounder  and  the  Echo  Processor .  The  position 
of  the  MPL  data  acquisition  system  is  shown  for  reference. 
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Figure  1.2  Sea  beam  transmit/ receive  geometry.  Computed  beam  pattern  cross-sections  in  the  athwartships  vertical  plane  cen¬ 
tered  on  the  array  and  in  the  vertical  plane  passing  through  the  ship’s  fore-aft  axis  are  shown  for  (a)  the  projector  array,  (b) 
the  receiver  array.  The  effect  of  Dolph-Chebyshev  amplitude  shading  is  also  illustrated,  (c)  is  a  summary  cartoon  showing  the 
angular  relationship  between  the  mainlobe  of  the  transmitted  beam  pattern  and  those  of  the  sixteen  preformed  beams. 
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angles  of  ±  10  degrees.  The  athwartships  beam  width  is  2  2/3  degrees.  Sea  Beam 
receives  with  sixteen  fixed  preformed  beams  obtained  by  electronically  steering  this  20 
x  2  2/3  degrees  beam  at  athwartships  intervals  of  2  2/3  degrees  between  ±  20  degrees 
of  incidence.  In  this  configuration  there  is  no  beam  along  the  ship’s  vertical  axis, 
rather  two  of  the  beams  point  at  1  1/3  degrees  on  either  side  of  this  axis. 

The  acoustic  energy  received  at  the  ship  comes  from  the  intersection  of  the 
transmit  and  receive  beam  patterns.  This  appears  in  Figure  1.2c  as  sixteen  squares  2 
2/3  degrees  on  a  side.  This  figure  is  only  meant  to  illustrate  the  angular  relationship 
between  the  mainlobe  of  the  transmitted  beam  pattern  and  the  mainlobes  of  the 
sixteen  preformed  beams.  Actual  footprints  are  not  rectangles  or  squares;  they  are 
ellipses  whose  areas  increase  away  from  vertical  incidence.  Since  depths  are  ideally 
determined  at  the  center  of  each  of  the  preformed  beams,  the  maximum  swath  width 
corresponds  to  73%  of  the  water  depth. 

The  beam  forming  operation  described  above  generates  sixteen  acoustic 
signals.  These  are  sent  to  the  EP  receivers  where  they  are  filtered,  rectified,  amplified 
and  transferred  to  the  Sea  Beam  computer  (Fig.  1.1).  Figure  1.3  shows  a  typical 
output  of  the  sixteen  EP  receivers.  Each  waveform  corresponds  to  a  preformed  beam 
and  is  accordingly  numbered  from  the  center  out  (T8)  on  port  and  starboard.  In  this 
figure  the  ridge  of  synchronous  returns  (see  arrow  labelled  "sidelobe" ),  corresponds  to 
energy  from  the  near-specular  direction  (tallest  return)  entering  the  sidelobes  of  the 
other  beams.  The  sixteen  bottom  return  signals  form  a  parabola,  indicating  a  flat 
portion  of  seafloor.  These  data  have  been  digitized  and  recorded  on  magnetic  tape 
with  a  separate  data  acquisition  system  (Fig.  1.1)  built  around  an  LSI  11/23 
minicomputer,  in  an  experiment  conducted  by  MPL  to  measure  acoustic  backscatter 
from  the  deep  seafloor.  They  have  proven  invaluable  to  evaluate  the  performance  of 
the  EP  because  Sea  Beam  only  retains  depths  and  cross-track  horizontal  distances. 
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Figure  1.3  Acoustic  signal  envelopes  of  the  sixteen  preformed  beams  at  the  output  of  the  Sea 
Beam  Echo  Processor  receivers .  The  time-axis  represents  seconds  after  transmission.  The 
vertical-axis  in  volts  represents  the  voltage  equivalent  of  the  sound  pressure  level  at  the 
receiver  array,  corrected  for  acoustic  transmission  losses  in  the  water  column  by  a  time-varied 
gain  (TVG).  No  roll  compensation,  recording  gain  or  receiver  gain  corrections  have  been 
applied  to  the  data  at  this  stage.  Such  data  are  recorded  digitally  on  magnetic  tape  every 
transmission  cycle,  along  with  time,  TVG  and  ship’s  roll. 
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In  the  Sea  Beam  computer,  sixteen  such  waveforms  are  simultaneously 
digitized  at  a  frequency  of  300  Hz  per  waveform.  This  corresponds  to  one  digitization 
cycle  every  3.33  ms  or  2.5  m  of  slant  range  assuming  a  sound  velocity  of  1500  m/s. 
Consequently,  slant  range  and  therefore  depth  determination  resolution  is  limited  to 
2.5  m.  While  it  digitizes  the  acoustic  data,  the  computer  also  performs  several  echo 
processing  tasks.  For  each  digitization  cycle  these  tasks  are:  receiver  gain  correction, 
refraction  correction,  roll  compensation,  detection  threshold  level  computation  and 
echo  detection.  Automatic  bottom-tracking  gates  (one  for  each  beam)  determine  a 
time  window  during  which  a  return  is  expected  on  any  one  beam  based  on  previous 

sounding  history.  A  return  is  detected  if  it  falls  within  the  gates  and  lies  above  the 
threshold. 

In  general,  the  threshold  level  is  computed  to  ride  above  the  noise,  above  the 
sidelobe  response  to  a  strong  specular  return  and  above  potential  noise  bursts 
interfering  with  bottom  echo  detection.  A  manual  threshold  can  also  be  entered  bv 
the  Sea  Beam  operator.  As  we  shall  see  in  the  following  sections,  thresholding  and 
gating  are  two  critical  operations  in  the  echo  processing.  A  more  detailed  description 
of  Sea  Beam’s  echo  processing  may  be  found  in  Appendix  B. 

For  each  roll  compensated  beam  having  sufficient  signal  to  noise  ratio,  a 
slant  range  R  is  calculated  by  computing  the  center  of  mass  of  all  the  detected  signal 
samples  for  that  beam,  and  by  multiplying  the  corresponding  arrival  time  by  750  m/s. 
Knowing  the  slant  range  R  and  the  stabilized  beam  angle  'k,  a  simple  calculation 
yields  the  depth  Z  and  the  cross-track  horizontal  distance  Y: 

Q 

Z  =  R  cos  'I'  and  Y  =  — —  R  sin  'I', 

Cn 

where  Ga  is  the  mean  sound  velocity  obtained  by  averaging  the  values  of  the  sound 
velocity  profile  from  the  surface  to  the  average  bottom  depth  (in  uncorrected  meters) 
and  C n  is  the  nominal  sound  velocity  in  water  (1500  m/s).  The  depth  Z  is  given  in 
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uncorrected  meters  referenced  to  a  sound  velocity  of  1500  m/s.  The  cross-track 
horizontal  distance  Y  is  a  true  distance  because  it  is  corrected  for  both  refraction  and 
travel  time. 

Finally,  the  (Z,Y)  coordinates  for  each  validated  beam  are  output  as  a  cross¬ 
track  bottom  profile  on  the  EP  cathode  ray  tube  (CRT)  display.  The  depths  are  also 
used  to  update  the  bottom-tracking  gates  on  each  beam  for  the  next  transmission 
cycle. 

The  depths  and  cross-track  horizontal  distances  for  each  transmission  cycle 
are  logged  on  a  magnetic  storage  medium  (disk  or  tape)  along  with  time  and  ship’s 
heading,  as  well  as  output  in  near-real  time  (  ~  1  minute  delay)  on  paper  as  a  contour 
chart  by  an  11"  digital  swath  plotter.  In  the  following,  we  will  refer  to  the  (Z,Y)  data 
as  the  raw  Sea  Beam  data. 

The  Sea  Beam  echo  processing  sequence  outlined  here  will  vary  depending  on 
the  EP  mode  chosen  by  the  operator.  Three  modes  are  available.  Mode  1  is 
essentially  a  start-up  mode  during  which  no  data  logging  or  contour  plotting  are 
performed.  The  EP  displays  the  vertical  beam  depth  and  the  CRT  shows  unprocessed 
echoes  on  the  sixteen  preformed  beams.  The  detection  threshold  used  in  mode  1  is  the 
highest  of  the  noise  threshold,  the  sidelobe  threshold  or  the  threshold  entered  by  the 
operator.  Mode  2  is  a  semi-automatic  EP  operation  with  data  logging  and  contour 
plotting.  The  CRT  displays  processed  data  in  the  form  of  a  cross-track  bottom  depth 
profile,  but  the  operator  controls  the  tracking  gates’  width  and  center.  Mode  3  is  a 
completely  automatic  version  of  mode  2.  It  is  the  mode  in  which  the  EP  usually 
operates  during  bathymetric  survey  work.  A  very  important  and  poorly  documented 
difference  exists  between  the  detection  threshold  level  determination  of  mode  1  and 
that  of  modes  2  and  3.  In  modes  2  and  3,  a  non-zero  threshold  level  input  by  the 
operator  supercedes  any  other  threshold  computation.  It  is  therefore  imperative  that 
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the  manual  threshold  be  set  to  zero  v/hen  in  mode  2  or  3.  Failure  to  do  so  results  in 
the  EP  tracking  the  sidelobe  response  any  time  a  specular  return  is  present  on  one  of 
the  sixteen  preformed  beams,  or  if  the  manual  threshold  is  set  high  enough,  in  loss  of 
data. 

On  most  ships  equipped  with  Sea  Beam  the  bathymetry  data  are  merged 
with  ship  navigation  (transit  satellite  navigation  and  dead  reckoning,  or  NAVSTAR 
Global  Positioning  System  navigation  when  available)  by  another  computer  aboard 
the  ship  and  recontoured  along  the  ship’s  track  on  a  30"  digital  plotter  with  a  delay 
time  of  about  two  minutes.  This  gives  surveyor  the  ability  to  effectively  control  ship 
navigation  and  track  spacing  by  looking  at  the  contours  plotted.  A  second  stage  of 
data  post-processing,  done  on  ship  or  ashore,  consists  of  adjusting  the  navigation  to  fit 
corresponding  contour  lines  on  adjacent  tracks,  and  of  regridding  the  entire  data  set 
to  produce  a  map.  When  navigation  comes  from  the  Global  Positioning  System  there 
is  virtually  no  need  for  adjustments.  These  operations  usually  smooth  the  raw  Sea 
Beam  soundings  by  averaging  along  track  over  a  certain  number  of  transmission 
cycles  (often  5)  to  produce  more  even  grid  spacing  along  versus  across  the  ship’s  track, 
thus  removing  most  of  the  jitter  apparent  on  near-real  time  swath  plots.  However 
when  system  errors  cause  bad  soundings,  the  resulting  fictitious  bathymetry  will  often 
not  average  out  as  we  shall  show  in  the  following  sections.  Therefore  in  order  to 
assess  the  validity  of  suspicious  Sea  Beam  bathymetry,  an  investigator  needs  to  refer 
to  the  raw  data  and  use  any  corroborating  information  available.  When  only  the  raw 
data  are  available,  as  is  often  the  case,  such  assessment  requires  a  clear  understanding 
of  the  processing  performed  by  the  Sea  Beam  computer  on  the  digitized  acoustic 
signals. 


14 


1.4.  SEA  BEAM  BATHYMETRIC  ARTIFACTS:  Examples,  Explanations 
and  Geological  implications 

Sea  Beam  has  been  used  extensively  in  the  past  several  years  to  study  the 
morphology,  tectonics,  volcanology  and  sedimentology  of  the  seafloor.  So  many  Sea 
Beam  surveys  have  been  run  that  a  complete  list  is  too  large  for  inclusion  here, 
therefore  we  only  reference  some  of  the  more  recent  works.  Bathymetric  charts 
produced  from  Sea  Beam  data  have  been  used  as  base  maps  for  more  detailed  studies 
using  deeply  towed  instrument  packages  such  as  MPL’s  Deep-Tow  [  Spiess  and 
Lonsdale,  1982;  Spiess  et  al.,  1984;  Hey  et  al.,  1985a  and  manned  submersibles  such  as 
Woods  Hole  Oceanographic  Institute’s  (WHOVs)DSRV  Alvin  and  DSRV  Cyana  of  the 
Institut  Franpais  de  Recherche  pour  PExploitation  de  la  Mer  (IFREMER)*  [e.g. 
Ballard  and  Francheteau,  1983;  Francheteau  and  Ballard,  1983].  Many  surveys 
covering  fairly  large  areas  (hundreds  of  square  km)  with  nearly  total  coverage  have 
lead  to  valuable  insights  into  the  processes  at  spreading  centers  [e.g.  Hey  et  al.,  1985a; 
Crane  et  al.,  1985;  Macdonald  et  al.,  1984;  Mammerickx,  1984;  Lonsdale,  1983], 
transform  faults  fe.g.  Gallo  et  al.,  1984;  Detrick  et  al.,  1984]  trenches  fe.g.  Shipley  and 
Moore,  1985;  Lewis  et  al.,  1984],  microplates  [e.g.  Hey  et  al  1985b,  Naar  and  Hey, 
1985  ,  seamounts  [e.g.  Fornari  et  al.,  1984],  and  submarine  canyon  systems  [e.g.  Lewis 
et  al.,  1984].  Sea  Beam’s  regional  depiction  of  the  seafloor  in  these  areas  has  been 
extremely  useful. 

In  most  cases,  the  finer  scale  Sea  Beam  bathymetry  is  dependable  and  is 
reproducible  on  overlapping  swaths.  However,  given  the  existence  of  the  bathymetric 
artifacts  discussed  in  this  paper,  investigators  should  be  cautious  when  studying 
bathymetric  details  on  the  scale  of  hundreds  to  thousands  of  meters.  This  is 
particularly  important  for  surveys  of  large  areas  where  bathymetric  and  structural 


*formeriy  Centre  National  pour  PExploitation  des  Oceans  ( CNEXO ) 
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data  are  interpolated  between  widely  separated  Sea  Beam  swaths.  Misinterpreting 
any  of  these  artifacts  as  true  bathymetric  features  could  also  result  in  positioning 
errors  when  the  vessel  is  navigated  by  comparison  of  real-time  bathymetry  with 
compiled  charts. 

In  the  following  we  discuss  three  types  of  bathymetric  artifacts  resulting 
from  echo  processing  errors.  These  errors  are  due  to  internal  factors  such  as  sidelobe 
interference  or  malfunction  of  the  bottom-tracking  gates  or  to  external  interference 
from  other  sound  sources.  In  each  case  we  present  evidence  of  artifacts  through  Sea 
Beam  data  samples,  explain  their  cause,  and  indicate  their  geological  implications. 

1.4.1  Sidelobe  interference 

Renard  and  Allenou  1979]  recognized  sidelobe  interference  as  a  potential 
problem  in  Sea  Beam  (e.g.  Fig.  21  in  their  paper).  The  interference  is  characterized  by 
small  apparent  slope  fluctuations  on  seafloor  dipping  perpendicular  to  the  ship’s  track 
and  it  typically  affects  only  a  few  beams.  A  more  serious  problem  occurs  when  the 
seafloor  surveyed  is  relatively  flat  and  Sea  Beam  renders  it  as  a  trough  (the  Tunnel’ 
effect  [Smith,  1983]).  This  is  seen  as  a  concave-up  arc  on  the  cross-track  CRT  bottom 
profile  in  Figure  1.4.  To  understand  this  artifact,  consider  the  Sea  Beam  acoustic 
data  shown  in  Figure  1.5  which  is  identical  in  format  to  Figure  1.3.  Note  that  the 
sidelobe  level  is  much  higher  in  this  figure  than  in  Figure  1.3.  If  the  EP  is  in  mode  3 
and  a  non-zero  manual  threshold  level  has  been  entered  by  the  Sea  Beam  operator, 
the  system  does  not  calculate  a  noise  or  a  sidelobe  threshold.  It  therefore  tracks  the 
sidelobe  response  when  present  and  when  above  the  manual  threshold  level.  Arrival 
times  are  then  synchronous  on  all  beams  as  if  coming  from  a  concave-up  horizontal 
half  cylinder.  Figure  1,6  shows  an  example  of  the  resulting  bathymetry. 

The  apparent  relief  of  the  Tunnel’  walls  in  this  example  ranges  from  40  m  to 
100  m,  although  theoretically  it  may  be  as  much  as  6%  of  the  water  depth. 
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Figure  1.4  Echo  Processor  CRT  display  showing  a  cross-track  bottom  profile  (solid  trace) 
characteristic  of  a  1 tunnel '  effect .  The  dashed  traces  represent  the  upper  and  lower  positions  of 
the  automatic  bottom-tracking  gates.  The  vertical  scale  is  200  m  per  division.  The  horizontal 
scale  is  compressed  to  accommodate  a  reception  beam  width  spanning  80  g  (40  *  actual  beam 
width  with  ±  20  Q  for  roll). 
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Figure  1.5  Acoustic  signal  envelopes  of  the  sixteen  preformed  beams  at  the  output  of  the  Sea 
Beam  EP  receivers.  The  format  is  identical  to  that  of  Figure  1.3.  The  ridge  of  synchronous 
returns  due  to  sidelobe  response  is  much  more  pronounced  in  this  figure  than  in  Figure  1.3 
because  it  is  due  to  a  stronger  near-specular  return  (starboard  beam  #1,  which  is  clipped  in 
this  figure) 
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Figure  1,0  The  * tunnel '  effect.  Four  portions  of  Sea  Beam  swaths  adjusted  for  navigation  are  shown  here  to  illustrate  the 
effects  of  a  non-zero  manual  threshold  when  the  Echo  Processor  runs  in  mode  2  or  3.  Sea  Beanes  rendition  of  the  bathymetry  is 
seen  as  a  trough  approximately  centered  on  the  ship’s  track.  This  trough  persists  through  changes  in  the  ship’s  course,  thereby 
indicating  its  artificial  nature.  Fictitious  gullies  also  appear  on  slopes  updip  as  well  as  downdip  (best  seen  in  the  upper  right 
section  of  this  figure).  Contour  interval  is  10  m  and  tick  marks  point  downhill. 
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The  actual  seafloor  morphology  in  this  area  is  not  precisely  known  because  the  MPL 
acoustic  data  acquisition  system  was  not  available  during  this  survey.  This  area  is 
believed  to  be  generally  flat  with  indications  of  roughly  North-South  abyssal  hill 
trends.  Such  ‘tunnels’  might  be  mistaken  for  troughs  between  abyssal  hills  or 
submarine  channels,  but  investigators  would  recognize  them  as  artificial  because  the 
trough  axes  follow  the  ship  track,  independent  of  course  changes. 

The  tunnel  effect  can  also  occur  when  a  zero  manual  threshold  has  been 
entered,  even  though  the  system  computes  a  noise  and  a  sidelobe  threshold.  In  the 
example  given  in  Figure  1.5,  we  identify  two  processes  which  combine  to  defeat  the 
sidelobe  rejection  scheme  outlined  in  the  appendices.  First  a  very  strong  specular 
return  was  received  at  the  hydrophones,  indicating  a  highly  reflective  seafloor.  Second 
the  EP  receiver  outputs  were  found  to  saturate  at  8.5  volts  rather  than  the  specified 
maximum  output  of  10  volts  [de  Moustier,  1985a],  As  a  result  the  peak  amplitude  on 
the  specular  return  is  clipped  (starboard  beam  in  Fig.  1.5).  The  sidelobe  threshold 
level  computed  on  a  clipped  peak  only  partially  removes  the  sidelobe  response,  and  the 
remaining  portions  of  sidelobe  response  bias  the  center  of  mass  calculation  in  their 
direction  given  the  comparatively  low  signal  to  noise  ratio  of  the  real  backscattered 
bottom  returns.  Eventually,  the  system  tracks  the  sidelobe  response  instead  of  the 
bottom,  creating  a  trough-like  feature.  The  limiting  case  is  that  of  a  mirror-like  hard 
surface  from  which  there  is  no  backscatter.  In  this  case,  one  would  see  only  a  strong 
specular  return  and  a  synchronous  ridge  of  sidelobe  returns.  However,  most  of  the 
time  the  bottom  offers  some  roughness  on  the  scale  of  Sea  Beam’s  12  cm  acoustic 
wavelength,  and  the  signal  to  noise  ratio  of  the  backscattered  returns  is  sufficient  to 
track  the  bottom  correctly. 

It  is  important  to  note  that  the  prerequisite  for  sidelobe  interference  is  a 
strong  near-specular  return  on  any  one  of  the  preformed  beams.  The  bottom  does  not 
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necessarily  have  to  be  flat  ^e.g.  Renard  and  Allenou’s  example].  In  cases  where  the 
sidelobe  response  is  well  separated  from  the  bottom  return  (e.g.  port  beam  #8  in  Fig. 
1.5),  it  usually  falls  outside  the  tracking  gates.  When  the  sidelobe  and  the  actual 
bottom  returns  are  close  together  or  overlapping,  as  is  usually  the  case  on  returns 
adjacent  to  the  near-specular  return  (e.g.  port  beam  #  1  in  Fig.  1.5),  the  system  has 
no  way  of  differentiating  between  sidelobe  response  and  bottom  return.  Rejecting  the 
sidelobe  response  will  most  likely  cancel  some  of  the  bottom  return,  resulting  in  a 
slightly  erroneous  depth  determination.  Likewise  the  computed  depth  is  in  error  if  the 
sidelobe  response  is  not  rejected.  The  errors  are  small  (~  5  m)  for  beams  oriented  in 
the  near-specular  directions,  and  increase  away  from  specular  incidence  due  to  the 
lengthening  of  the  backscattered  return  signal  duration  (pulse  stretching)  with  both 
beam  angle  and  depth. 

1.4.2  Interference  from  external  sound  sources 

External  sound  sources  interfere  with  the  Sea  Beam  system  when  they 
transmit  while  Sea  Beam  is  receiving  echoes  from  the  seafloor.  Figure  1.7  shows  an 
example  of  a  3.5  kHz  echo-sounder  interference  as  seen  in  the  acoustic  data.  It 
appears  as  a  synchronous  ridge  across  the  sixteen  preformed  beams.  This  is  a  classical 
example  of  a  noise  burst.  As  for  the  sidelobe  interference,  the  dynamic  thresholding 
used  to  reject  such  noise  bursts  has  side  effects  which  produce  fictitious  bathymetry, 
examples  of  which  can  be  seen  in  Figure  1.8.  The  portion  of  the  noise  burst  which  is 
well  separated  from  the  actual  bottom  returns  are  effectively  rejected  by  dynamic 
thresholding  or  by  gating.  However  where  signal  and  noise  burst  overlap,  canceling 
the  noise  burst  also  cancels  part  of  the  signal  and  skews  the  center  of  mass  calculation 
for  that  return.  As  the  noise  burst  slowly  progresses  through  the  reception  cycle  over 
a  number  of  pings,  bathymetric  peaks  appear  on  the  contoured  output  in  the  direction 
of  the  beams  which  point  away  from  specular  incidence.  The  near-specular  directions 
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Figure  1.7  Evidence  of  a  3.5  kHz  echo-sounder  transmitting  during  a  Sea  Beam  reception  cycle. 
The  corresponding  noise  burst  appears  as  a  synchronous  ridge  across  all  sixteen  preformed 
beams.  The  format  is  the  same  as  that  of  Figure  1.3.  A  noise  burst  ridge  differs  from  a 
sidelobe  response  ridge  in  that  the  levels  of  the  peaks  are  more  or  less  constant  for  the  former 
while  a  marked  difference  in  level  exists  between  the  specular  return  and  its  corresponding 
sidelobe  response  (Figures  1.3  and  1.5).  The  differences  in  level  seen  in  this  figure  are  due  to 
differences  in  receiver  gains  which  were  not  corrected. 
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Figure  1.8  Examples  of  contoured  swath  plots  showing  the  results  of  external  sound  source 
interference.  Artifacts  can  be  recognized  as  individual  peaks  on  one  or  both  sides  of  the  ship’s 
track  (center  line  in  all  3  plots;.  Contour  interval  is  10  m  and  tick  marks  point  downhill. 


23 


are  not  affected  as  much  since  pulse  stretching  is  minimal,  thereby  reducing  the 
margin  for  error. 

The  bathymetric  peaks  are  typically  short  wavelength  (hundreds  of  meters) 
and  may  vary  in  amplitude  from  tens  to  hundreds  of  meters.  The  more  pronounced  of 
these  peaks  are  often  clearly  spurious  and  extremely  steep  and  sharp;  such  features 
are  geologically  unlikely  and  investigators  will  readily  disregard  them.  Smaller 
amplitude  artifacts  are  less  obvious,  and  might  be  mistaken  for  small  volcanic  cones 
or  large  hydrothermal  mounds.  One  common,  although  not  ubiquitous  characteristic 
of  these  artifacts  is  the  simultaneous  occurrence  of  more  than  one  in  different  parts  of 
the  swath.  Investigators  aware  of  the  potential  for  these  phenomena  are  unlikely  to 
misinterpret  them.  Seismic  sound  sources  such  as  waterguns  produce  similar  effects, 
but  no  observable  interferences  have  been  reported  with  airguns,  probably  because 
they  do  not  output  enough  acoustic  energy  in  the  12  kHz  frequency  band  [Smith, 
1983]. 

A  special  case  of  interference  from  external  sound  sources  exists  for  12  kHz 
bottom  transponders.  Figure  1.9  shows  an  example  of  such  interference  with  evidence 
of  a  transponder  trace  on  the  corresponding  analog  center  beam  depth  profile.  The 
fiat  sedimentary  bottom  over  which  this  data  was  taken  illustrates  the  progression  of 
the  interference.  The  interference  enters  the  outer  beams’  tracking  gates  while  falling 
outside  those  of  the  near-specular  beams.  This  is  evidenced  in  Figure  1.9a  by  a 
central  ridge  followed  by  two  small  mounds  on  either  side  of  the  ship’s  track.  The 
small  mounds  would  be  difficult  to  identify  as  artificial,  were  it  not  for  evidence  from 
the  analog  record  (Fig.  1.9b)  which  shows  the  transponder  trace  intersecting  the 
center  beam  depth  profile  at  the  corresponding  time.  Due  to  their  small  size,  these 
artifacts  would  probably  not  be  considered  very  significant,  although  some  might 
mistake  them  for  satellite  cones  or  hydrothermal  mounds. 
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Figure  1.9  Interference  from  a  12  kHz  bottom-moored  transponder  during  a  Deep  Tow  survey. 
(a)  Sea  Beam  near-real  time  contoured  output.  Contour  interval  is  10  m  and  tick  marks  point 
downhill.  The  center  line  represents  the  ship’s  track.  On  this  line,  short  ticks  above  the  ship’s 
track  are  spaced  2  minutes  apart,  long  ticks  refer  to  information  at  the  top  of  the  plot  [time 
(hour/min/sec),  ship’s  heading,  contour  interval  in  meters]  and  short  ticks  centered  on  the 
track  refer  to  center  beam  depth  in  meters  indicated  at  the  bottom,  (b)  analog  graphic 
recorder  output  displaying  Sea  Beam’s  center  beam  depth  profile  and  the  trace  of  the  12  kHz 
transponder.  The  horizontal  scale  is  matched  in  time  to  that  of  Figure  (a)  The  artifact  can  be 
seen  in  (a)  as  the  two  small  mounds  on  either  side  of  the  ship’s  track. 
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The  situation  of  this  example  is  uncommon  because  the  ship  was 
maneuvering  at  about  1.5  knots  over  a  bottom  transponder  network  while  towing 
Deep  Tow  instrument  package.  However,  it  may  become  more  common  with  the  Sea 
Beam  system  installed  on  WHOI’s  R/V  Atlantis  II,  the  mother  ship  for  the  manned 
DSRV  Alvin  which  is  often  navigated  using  12  kHz  transponders.  At  normal  survey 
speeds  (~  10  knots)  this  artifact  would  be  greatly  reduced.  Similar  artifacts  due  to 
interference  from  the  direct  or  the  bottom  bounced  signal  of  a  12  kHz  pinger  have  also 
been  noted  during  dredging  or  coring  operations. 

1.4.3  ‘Omega’  effects  and  data  gaps 

Most  Sea  Beam  users  are  aware  of  the  possibility  of  sidelobe  or  external 
sound  source  interference  in  the  system.  A  lesser  known  and  more  insidious  artifact 
has  been  found  to  occur  on  sloping  bottoms  producing  contours  resembling  the  capital 
Greek  letter  omega  (fi)  [Kleinrock  et  al,  1984  or  data  gaps.  They  are  generally 
characterized  by  an  arcuate  plateau  followed  by  a  steep,  curvilinear  scarp.  They 
occur  within  a  single  Sea  Beam  swath,  commonly  near  the  center,  and  have  lateral 
dimensions  of  hundreds  to  thousands  of  meters.  The  plateaus  may  be  peaked  (Fig. 
l.lOe)  or  flat  (Fig.  l.lOd).  The  scarp  may  be  semicircular  (Fig.  l.lOd  --  the  classic  H 
shape)  or  irregular  (Fig.  l.lOf).  'Omegas’  sometimes  evolve  into,  or  are  associated 
with  data  gaps  (Figs.  1.10b  and  1.10g),  and  can  be  created  on  sides  of  seamounts  (Fig. 
1  10c)  as  well  as  relatively  straight  scarps. 

Figure  1.10a  shows  the  problem  clearly.  In  this  case,  the  same  portion  of 
seafloor  has  been  surveyed  in  3  different  directions.  The  arrows  indicate  the  direction 
of  ship  travel.  In  Figures  l.lOa.l  and  10a. 2  Sea  Beam’s  rendition  of  the  bathymetry  is 
nearly  identical  for  opposite  ship  courses  in  the  along  slope  direction.  The  bathymetry 
is  markedly  different  when  the  ship  track  runs  downdip  (  across  a  slope  in  the  down 
hill  direction)  (Fig.  1.10a. 3).  In  all  the  colored  contour  plots  shown,  contour  lines  have 
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Figure  1.10  'Omegas’  and  gaps.  Tick  marks  point  downhill.  Contoured  sections  not  starred 
are  original  Sea  Beam  data.  Contoured  sections  marked  with  a  star  are  the  result  of  repro¬ 
cessing  the  acoustic  data  recorded  digitally  with  the  MPL  system.  Our  simplified  echo  process¬ 
ing  technique  does  not  include  ray-bending  corrections  and  arrival  times  are  determined  by  the 
first  arrival  above  a  preset  threshold.  The  threshold  level  is  selected  after  visual  inspection  of 
the  roll-compensated  acoustic  data.  Recomputed  depths  and  cross-track  distances  are  there¬ 
fore  in  uncorrected  meters  referenced  to  a  sound  velocity  of  1500  m/s.  Although  crude,  this 
processing  method  suffices  to  prove  the  fictitious  character  of  Sea  Beam’s  contoured  bathy¬ 
metry  shown  in  (a. 3),  (b),  (d),  (e),  (f),  and  (g)  We  do  not  show  a  recomputed  version  of  (c) 
because  the  corresponding  acoustic  data  was  only  recorded  every  five  pings.  This  was  enough 
to  confirm  the  ‘omega’  effect  but  contour  resolution  was  seriously  degraded  by  the  five  ping 
decimation. 
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been  smoothed  by  averaging  over  five  transmission  cycles.  While  inspecting  the  raw 
Sea  Beam  data,  we  noticed  unrealistic  variations  in  depth  from  one  ping  to  the  next 
as  well  as  missing  soundings  in  the  data  of  Figure  1.10a. 3.  No  evidence  of  external 
interference  was  found  in  the  acoustic  data.  After  computing  depths  and  cross-track 
distances  from  this  acoustic  data,  we  contoured  them  using  the  same  post-processing 
software  used  throughout  Figure  1.10.  The  resulting  bathymetry  shown  in  Figure 
1.10a. 4  matches  that  seen  in  Figures  1  10a. 1  and  1  10a. 2.  Sea  Beam  was  clearly  in 
error  when  the  ship  track  ran  downdip. 

A  combination  of  three  factors  may  be  responsible  for  this  artifact.  First  of 
all,  the  automatic  bottom-tracking  gates  do  not  open  fast  enough  upon  a  sudden 
change  of  bottom  slope.  As  a  result,  data  are  lost  for  points  falling  outside  the  gates. 
Secondly,  when  going  downdip  across  a  slope  the  fore-aft  transmit  beam  pattern 
geometry  (Fig.  1.2. a)  is  such  that  acoustic  energy  from  the  sidelobes  may  ensonify  the 
slope  in  the  specular  direction.  Although  this  transmitted  energy  is  about  25  dB  lower 
than  that  transmitted  in  the  mainlobe  in  the  true  vertical  direction,  it  becomes 
significant  due  to  the  angular  dependence  of  backscattering.  Measurements  have 
shown  that  one  can  expect  a  drop  of  10  to  15  dB  in  the  acoustic  backscatter  between 
normal  (specular)  incidence  and  20  degrees  incidence  [Patterson  1969,  Urick  1983]. 
For  a  flat  bottom  the  acoustic  backscatter  due  to  sidelobe  transmitted  energy  is 
negligible  compared  to  the  returns  due  to  the  main  beam  in  the  specular  direction. 
This  is  no  longer  true  when  the  bottom  slope  is  such  that  the  sidelobes  of  the  transmit 
beam  pattern  point  in  the  specular  direction  as  in  this  case.  Because  the  receiving 
beam  pattern  is  20  degrees  wide  in  the  fore-aft  direction  (Fig.  1.2b),  early  returns  are 
received  as  seen  in  Figure  1.11.  Given  the  proper  threshold  and  sufficiently  narrow 
tracking  gates  (which  is  the  case  when  the  bottom  has  remained  relatively  flat  for 
some  time),  the  system  tends  to  track  these  early  returns,  creating  a  plateau-like 
feature.  At  some  point  the  signal  to  noise  ratio  of  these  returns  becomes  too  low  and 
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Figure  1.11  Evidence,  of  transmit  beam  pattern  sidelobe  interference  in  the  digitized  acoustic 
data.  Early  arrivals  corresponding  to  near-specular  returns  from  transmited  sidelobe  energy 
are  best  seen  on  starboard  beams  1-3.  These  data  correspond  to  the  ‘omega’  effect  shown  in 
Figure  l.lOe  The  format  is  the  same  as  that  of  Figure  1.3. 
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the  system  fails  to  detect  an  echo  until  the  tracking  gates  open  wide  enough  to 
recover  the  real  bottom.  Hence  a  sharp  drop  in  depth  results  at  the  end  of  the 
plateau  as  shown  in  the  cartoon  of  Figure  1.12a. 

Consider  two  possible  tracking  gate  conditions:  a  normal  condition  where  the 
instantaneous  bottom  profile  is  contained  within  the  gates  (Fig.  1.12b),  and  a 
condition  where  part  of  the  profile  falls  outside  the  gates  (Fig.  1.12c).  The  latter 
produces  a  data  gap  as  seen  in  Figures  1.10b  and  1.10g  where  the  onset  of  an  ‘omega5 
effect  immediately  precedes  the  gap.  Apparently,  the  dip  of  the  bottom  increased  too 
rapidly  for  the  ‘omega5  effect  to  fully  develop,  and  a  gap  appeared  because  the  gates 
simply  could  not  open  fast  enough.  Such  gaps  exist  in  Sea  Beam  data  on  updip  as  well 
as  downdip  ship  tracks,  however  in  our  data  we  have  seen  ‘omega5  effects  only  for 
downdip  ship  tracks.  This  was  confirmed  at  sea  when  an  observed  ‘omega5  effect  on  a 
downdip  track  was  immediately  resurveyed  updip,  and  no  ‘omega5  was  detected.  The 
most  likely  explanation  for  this  asymmetry  comes  from  the  fact  that  the  gates  are 
always  lagging  upon  a  sudden  change  in  bottom  slope.  Downdip,  the  gates  track  from 
the  left  in  Figure  111  and  they  are  therefore  likely  to  track  early  returns.  Updip  the 
gates  track  from  the  right  in  Figure  1.11  so  they  have  a  better  chance  to  track 
bottom  returns  instead  of  early  arrivals.  Also  the  gates  have  more  time  to  adjust  at 
the  base  of  rising  slopes  due  to  the  accumulation  of  talus.  We  cannot  specify  a  slope 
range  for  which  an  ‘omega5  effect  occurs  because  this  effect  varies  with  ship’s  speed 
and  depends  on  the  sidelobe  level  on  the  transmit  beam  pattern  of  the  Sea  Beam 
system  considered.  As  ship  speed  is  reduced,  the  tracking  gates  have  more 
transmission  cycles  to  adjust  to  a  sudden  drop  in  slope  and  the  ‘omega5  effect  is  less 
likely.  Our  data  shows  ‘omega5  effects  on  slopes  between  30  and  45  degrees  for  ship 
speeds  of  about  10  knots,  but  similar  though  subtler  artifacts  seem  to  appear  on 
gentler  slopes,  perhaps  as  low  as  15  degrees. 
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Figure  1.12  Cartoons  of  bottom  profiles  associated  with  ‘omegas’  and  gaps .  (a)  Along  track 
bottom  profile  of  an  ‘omega’  artifact.  Sea  Beam’s  rendition  of  the  bathymetry,  shown  by  the 
dashed  line,  is  a  plateau  followed  by  a  steep  scarp.  (b,c)  Bottom-tracking  gates  conditions  as 
would  be  seen  on  the  Echo  Processor’s  CRT.  In  (b)  the  athwartships  bottom  profile  (solid  line) 
lies  inside  the  tracking  gates  (dashed  lines).  This  is  the  normal  mode  of  operation  when  the 
Echo  Processor  is  in  mode  3.  It  is  also  what  one  would  see  were  an  ‘omega’  effect  present.  In 
(c)  the  athwartships  bottom  profile  lies  partially  outside  the  tracking  gates  and  the 
corresponding  data  points  are  lost.  This  situation  is  characteristic  of  data  gaps. 
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The  tracking  gates  and  the  transmit/receive  acoustic  geometry  are  the  two 
main  factors  contributing  to  the  ‘omega7  effect.  A  third  factor  is  related  to  the  half- 
hour  calibration  of  the  EP  receivers.  In  several  instances,  we  found  that  this 
calibration  occurred  immediately  prior  to  an  ‘omega7  effect.  Figure  1.13  shows  the 
onset  of  a  receiver  gain  calibration  sequence  just  at  the  end  of  a  reception  cycle.  The 
following  transmission  cycle  showed  only  the  calibration  signals.  Inspection  of  the  raw 
Sea  Beam  data  showed  that  no  data  had  been  logged  for  these  two  cycles  even  though 
one  would  have  expected  the  first  (Fig.  1.13)  to  have  been  processed  by  the  EP.  As  a 
consequence,  data  are  lost  for  two  transmission  cycles  every  half  hour.  Moreover  the 
tracking  gates  are  not  updated  during  this  time.  A  coincidental  increase  in  bottom 
slope  puts  the  EP  in  a  difficult  bottom-tracking  situation  which,  given  the  appropriate 
slope  angle,  ship  direction  and  signal-to-noise  ratio  in  the  acoustic  returns,  generates 
an  ‘omega7  effect. 

Of  all  the  artifacts  discussed  here,  ‘omega7  effects  are  the  most  likely  to 
mislead  investigators  because  they  often  appear  as  geologically  plausible  volcanic, 
tectonic  or  even  depositional  (mass-wasting)  features.  An  ‘omega7  on  the  side  of  a 
seamount  as  in  Figure  1.10c,  could  possibly  be  mistaken  for  a  flank  or  satellite 
construct.  Irregular  ‘omega7-type  artifacts  (Fig.  l.lOf)  appearing  on  what  are  actually 
relatively  straight  scarps  might  be  misinterpreted  as  changes  or  variations  in 
structural  trend.  This  could  result  in  errors  in  determining  the  tectonic  character  and 
evolution  of  an  area.  Other  ‘omegas7  (Fig.  LlOe)  might  be  mistaken  for  volcanic 
constructional  features  and  incorrect  conclusions  could  be  reached  regarding  post-  or 
syn-tectonic  volcanism  on  scarps  such  as  fracture  zone  walls,  rift  valley  walls, 
pseudofault  walls,  abyssal  hill  scarps,  caldera  walls,  etc.  In  addition,  some  ‘omegas7 
might  be  mistaken  for  serpentinite  diapirs,  while  others  (for  example  at  trenches  or 
submarine  canyons)  could  be  erroneously  identified  as  slump  blocks  or  other  mass- 
wasting  deposits.  When  dealing  on  scales  of  hundreds  of  meters  to  several  kilometers, 


32 


EARLY 


ARRIVALS 


Seco/ia 


BOTTOM 
RETURNS 


CALIBRATION 

SIGNAL 


Figure  1.13  Receivers  gam  calibration.  The  raised  signals  seen  at  the  end  of  each  of  the  six¬ 
teen  preformed  beam  returns  represent  calibration  signals  injected  into  each  receiver  to  deter¬ 
mine  its  gain  setting.  Here  one  can  appreciate  the  usefulness  of  the  Echo  Processor  gain 
correction  which  brings  all  receivers  to  a  common  gain  level.  The  format  of  this  figure  is 
identical  to  that  of  Figure  1.3. 
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failing  to  recognize  ‘omegas’  as  artifacts  could  lead  to  errors  in  geologic 
understanding,  because  they  might  suggest  unexpected  volcanism  or  tectonism  in 
supposedly  inactive  areas.  The  implications  of  these  possible  misinterpretations  are 
very  important. 

Suspicious  features  which  have  the  characteristic  shape  of  ‘omegas’  have  been 
observed  in  data  from  every  Sea  Beam  survey  we  have  investigated  thus  far.  Many 
geophysical  surveys  are  run  orthogonal  to  the  tectonic  fabric  because  important 
variations  in  magnetic,  seismic,  gravity  and  bathymetric  data  often  are  found  in 
cross-strike  profiles.  Unfortunately,  because  ‘omegas’  are  found  on  downdip  tracks, 
this  type  of  survey  pattern  increases  the  probability  of  occurrence  of  these  artifacts. 
In  an  effort  to  quantify  this  probability,  we  have  analyzed  data  from  such  a  survey 
(Fig.  1.14).  Figure  L14  shows  Sea  Beam  data  from  the  propagating  rift  at  95.5  °W  on 
the  Cocos-Nazca  spreading  center  [Hey  et  al.,  1985a].  In  this  figure,  all  areas  where 
the  gradient  of  the  bathymetry,  as  detected  by  Sea  Beam,  exceeds  a  specified  slope  are 
darkened.  Figure  1.14a  is  included  mainly  to  show  the  data  density  and  the  overall 
tectonic  structure.  The  ‘omega’  effect  was  discovered  while  analyzing  this  dense  data 
set  with  several  overlapping  swaths  and  the  ‘omegas’  shown  in  Figure  1.10  are 
examples  of  artifacts  that  Hey  and  coworkers  removed  from  their  data.  In  eight  cases 
for  which  we  initially  suspected  the  ‘omega’  bathymetry  to  be  false,  and  then  studied 
the  acoustic  data,  our  suspicions  were  confirmed.  By  checking  the  raw  Sea  Beam  data, 
we  have  identified  eight  others.  We  then  estimated  the  probability  of  encountering 
‘omegas’  on  downdip  tracks  over  fairly  steep  slopes.  We  have  visually  examined  the 
computer-generated  Sea  Beam  20  m  contour  plots  and  identified  all  the  ‘omegas’ 
which  we  feel  confident  are  artifacts  (many  questionable  examples  were  also  found, 
but  not  included  in  the  exercise).  Triangles  in  Figure  1.14b  mark  the  locations  where 
features  on  this  plot  are  associated  with  ‘omegas’  seen  in  the  contour  plots.  Knowing 
the  direction  of  ship’s  travel,  we  were  able  to  distinguish  on  this  figure  those  downdip 
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slope  crossings  which  have  ‘omega’  artifacts  and  those  which  do  not.  For 
approximately  50%  of  all  occurrences  where  the  ship  steamed  downdip  across  slopes 
which  Sea  Beam  detected  as  being  greater  than  30  degrees,  ‘omega’  artifacts  are 
present.  Though  we  realize  these  estimates  are  rough,  the  salient  point  is  clear: 
‘omega’  artifacts  can  be  very  common  in  conventional  surveys  which  run 
perpendicular  to  the  tectonic  fabric.  As  mentioned  above,  the  frequency  of  occurrence 
of  these  artifacts  may  vary  with  each  Sea  Beam  system  depending  on  the  sidelobe 
level  of  its  transmit  beam  pattern,  as  well  as  with  ship’s  speed. 

Our  detailed  analysis  of  these  artifacts  has  concentrated  on  the  SIO  Sea 
Beam  system  because  of  the  availability  of  its  acoustic  data.  But  we  have  observed 
‘omega’-like  features  in  bathymetric  data  collected  with  systems  aboard  the  R/V’s 
Conrad  (Lamont-Doherty  Geological  Observatory),  Surveyor  (National  Oceanographic 
and  Atmospheric  Administration),  and  Jean  Charcot  (IFREMER).  Sea  Beam 
investigators  who  see  suspicious  features  with  characteristics  similar  to  the  ‘omegas’ 
shown  here  would  be  prudent  to  survey  the  sites  with  crossing  Sea  Beam  swaths  for 
confirmation  before  attributing  them  great  significance  or  planning  higher  resolution 
studies.  If  such  features  are  not  recognized  as  potentially  important  until  after  the 
survey,  the  raw  Sea  Beam  data,  should  be  checked  looking  for  unrealistic  depth 
changes  from  one  ping  to  the  next,  and  for  missing  data  points  which  indicate  that  the 
system  lost  tracking  of  the  bottom  during  that  time. 

Although  we  have  only  analyzed  data  from  the  Sea  Beam  system,  we  believe 
similar  multibeam  echo-sounders  might  output  the  same  artifacts.  The  U.S.  Navy 
SASS  system  has  been  in  operation  since  1965,  and  some  of  its  data  has  been 
declassified  for  use  on  the  Mid-Atlantic  Ridge  Rift  Valley  Phillips  and  Fleming,  1978; 
Ballard  and  van  Andel,  1977]  and  on  the  Galapagos  Rift  at  86  0  W  [van  Andel  and 
Ballard,  1979;  Crane  and  Ballard,  1980].  Comparison  of  SASS  bathymetry  with  Deep 
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Tow  bathymetry  [Crane,  1978;  van  Andel  and  Ballard,  1979]  seems  to  indicate  that 
‘omega’-like  artifacts  exist  in  SASS  data.  Data  gaps  and  onsets  of  ‘omegas’  similar  to 
those  of  Figure  1.10b  are  also  apparent  in  Phillips  and  Fleming  [1978]  (e.g.  Figure  3D 
in  their  paper). 

L5.  POSSIBLE  CORRECTIONS 

Depths  and  cross-track  horizontal  distances  cannot  be  recomputed  as  a  post¬ 
processing  operation  unless  the  acoustic  data  is  digitized  and  recorded  on  tape  as  was 
done  for  our  data.  Therefore,  investigators  discovering  fictitious  bathymetry  in  their 
data  have  no  alternative  but  to  disregard  the  portion  of  data  affected.  Also,  these 
artifacts  occur  too  infrequently  to  warrant  recording  of  the  acoustic  data  on  a 
routine  basis.  Rather  than  relying  on  data  reprocessing,  it  seems  more  sensible  to  deal 
with  the  problems  at  their  source.  In  the  following  we  suggest  a  number  of  solutions 
to  the  problems  discussed  in  the  previous  section. 

L5.1  Possible  Corrections  for  Sidelobe  interference 

When  a  ‘tunnel’  effect  develops  during  a  survey,  it  is  common  practice  to 
switch  the  EP  from  mode  3  to  mode  2  [Smith,  1983].  As  a  result  the  automatic 
tracking  gates  open  to  their  maximum  (upper  and  lower  limits  of  the  CRT  display) 
according  to  parameters  set  by  the  Sea  Beam  operator.  A  subsequent  return  to  mode 
3  resumes  automatic  tracking.  This  method  has  proven  effective  in  dealing  with  the 
tunnel  effect  which  is  essentially  sidelobe  response  related.  However  it  relies  on  the 
vigilance  of  an  operator  and  it  will  not  correct  the  effect  pointed  out  by  Renard  and 
Allenou  [19791.  In  general  the  current  sidelobe  response  suppression  technique  suffers 
from  the  saturation  in  the  EP  receivers.  A  simple  modification  of  the  detection 
amplifiers  may  solve  the  saturation  problem,  but  use  of  logarithmic  detection 
amplifiers  to  increase  the  EP  receivers’  dynamic  range  seems  desirable  [de  Moustier, 


37 


18854  At  presenl  there  is  no  ‘he  performance  of  the  amplitude 

shading  in  the  receiver  array.  As  an  example,  we  changed  the  shading  coefficient  by 

30  o  on  four  array  elements  in  the  computed  beam  pattern  of  Figure  1.2b.  It  brought 
the  sidelobe  level  from  30  to  23  dB  below  the  mamlobe.  This  may  not  appear  as 
significant  since  the  sidelobe  threshold  computation  is  based  on  a  value  of  12  dB  below 
a  peak  amplitude,  which  is  approximately  the  sidelobe  level  of  an  unshaded  array  (Fig. 
1-b).  However,  we  believe  that  ensuring  optimum  performance  of  the  array 
amplitude  shading  can  only  benefit  any  subsequent  sidelobe  response  rejection  scheme. 

As  the  EP  works  on  the  rectified  envelope  of  the  return  signals,  it  has  no  way 
of  differentiating  between  sidelobe  response  and  actual  bottom  return  when  the  two 
overlap.  To  tell  them  apart  requires  phase  information  which  is  not  available  to  the 
EP  in  the  current  mode  of  operation.  One  way  to  deal  with  this  problem  would  be  to 
erodyne  (multiply  by  an  external  oscillator  frequency  and  filter  in  the  desired 
frequency  band  (Clay  and  Medwin,  1977J)  each  of  the  sixteen  preformed  beam 
channels  to  obtain  16  channels  of  complex  data  (32  channels  of  real  data).  These 
could  be  digitized  and  processed  as  currently  done  the  EP.  It  would  then  be 
possible  to  apply  advanced  adaptive  filtering  techniques  [MeCool  and  Widrow,  1977]  to 
effectively  cancel  sidelobe  response  as  well  as  noise  bursts  while  retaining  the  real 

bottom  return  signal.  Of  course,  such  a  scheme  might  be  hampered  by  the  processing 
time  required  and  it  needs  to  be  tested. 

1.5.2  Possible  Corrections  for  Interference  from  external  sound  sources 

Whenever  the  3.5  kHz  echo-sounder  is  run  in  conjunction  with  Sea  Beam 
aboard  the  R/  V  Thomas  Washington,  the  analog  graphic  recorder  is  used  to  display 
both  the  Sea  Beam  center  beam  profile  and  the  3.5  kHz  ech«ounder  outgoing  pulse. 
Interference  occurs  when  the  corresponding  signal  traces  intersect  To  prevent  this,  it 
is  necessary  to  phase  delay  the  3.5  kHz  outgoing  pulse  enough  to  keep  the  two  traces 
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separated  [Smith,  1983].  This  method  is  not  entirely  reliable  since  it  requires  an 
operator.  A  more  reliable  method  consists  of  using  a  simple  electronic  circuit  which 
gates  out  3.5  kHz  transmission  whenever  a  Sea  Beam  reception  cycle  is  in  progress. 
Such  a  device  reportedly  works  well  on  the  R/V  Conrad  rTyce,  1984j.  Unfortunately, 
this  device  will  not  prevent  interference  from  12  kHz  transponders,  pingers  or  seismic 
sound  sources.  The  latter  usually  cannot  be  phase  delayed  for  mechanical  reasons 
(constant  pulse  energy  requirement)  as  well  as  data  post-processing  reasons  (constant 
firing  rate  requirement).  A  solution  which  takes  into  account  the  transmission  rate 
requirements  of  all  possible  underwater  sound  sources  available  aboard  a  ship  can  be 
implemented  on  the  shipboard  computer.  With  a  knowledge  of  the  water  depth,  the 
computer  would  decide  the  best  firing  sequence  necessary  to  keep  the  sound  sources 
from  interfering  with  each  other  (J.L.  Abbott,  SIO:  personal  communication). 

1.5.3  Possible  Corrections  forlOmega’  effects  and  data  gaps 

Omegas'  and  data  gaps  are  dealt  with  in  the  same  way  the  ‘tunnel’  effect  is, 
by  switching  the  EP  from  mode  3  to  2  and  back  again.  However,  there  is  no  way  to 
detect  an  ‘omega’  effect  in  real-time  since  the  cross-track  bottom  profiles  on  the  CRT 
appear  to  be  within  the  gates,  and  by  the  time  evidence  of  it  is  seen  on  the  swath  plot 
it  is  usually  too  late  to  correct  anything.  The  automatic  tracking  gate  software  was 
modified  by  the  manufacturer  in  October  1984  on  the  system  installed  aboard  the 
German  Rf  V  Polarstern  (W.  Capell,  General  Instrument  Corp:  personnel 
communication).  The  changes  consisted  of  increasing  the  minimum  allowable  width 
for  each  gate  and  enabling  a  faster  rate  of  change  of  the  gates  from  ping  to  ping.  Sea 
Beam  systems  installed  since  October  1984  benefit  from  this  modification  which  has 
proven  effective  in  substantially  reducing  the  problems  of  data  gaps  and  ‘omega’ 
effects.  However,  the  widening  of  the  bottom-tracking  gates  tends  to  decrease  the 
depth  determination  accuracy  on  the  outer  beams  because  the  remaining  sidelobe 
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response  is  no  longer  gated  out  on  those  beams. 

Because  of  the  side  effects  of  the  EP  receivers’  calibration,  mentioned  in 
Section  1.4.3,  we  recommend  that  an  additional  change  be  made  to  allow  the  gates  to 
widen  during  a  calibration  cycle.  This  way  a  coincidental  increase  of  bottom  slope 
will  easily  be  accommodated  by  the  EP  upon  return  to  a  normal  reception  cycle.  Also, 
it  would  be  useful  to  have  the  half  hour  calibration,  which  is  triggered  upon  interrupt 
from  the  Sea  Beam  computer  clock,  wait  for  the  completion  of  the  transmission  cycle 
in  progress  and  avoid  situations  such  as  that  of  Figure  1.13.  Data  would  then  be  lost 
for  only  one  transmission  cycle  and  the  updating  of  the  tracking  gates  would  be  more 
reliable. 

Finally,  as  for  the  receiver  array,  some  measure  of  the  performance  of  the 
projector  array  amplitude  shading  seems  necessary.  At  present  the  system  tests  the 
performance  of  the  power  amplifiers  on  an  all  or  nothing  (blown  fuse)  basis.  Tyce 
1984  reported  deviations  from  the  manufacturer’s  specifications  by  as  much  as  40% 
on  the  outputs  of  four  projector  elements  for  the  system  installed  aboard  the  Rj  V 
Conrad.  For  comparison  a  change  of  40%  in  the  shading  coefficients  of  four  elements 
in  the  computed  beam  pattern  (Fig.  1.2a)  moved  the  siaelobe  level  from  30  dB  to  22 
dB  below  the  mainlobe.  Such  levels  will  definitely  enhance  the  signal  to  noise  ratio  of 
early  specular  returns  discussed  in  Section  1.4.3  (Fig.  1.11),  increasing  the  probability 
of  ‘omegas’. 

1.6.  CONCLUSIONS 

In  conclusion,  we  would  like  to  stress  the  importance  of  a  clear  understanding 
of  the  capabilities  and  limitations  of  the  Sea  Beam  system  when  analyzing  its  output. 
We  fully  recognize  the  value  of  Sea  Beam  bathymetry  in  conducting  a  survey  and  in 
describing  abyssal  morphology.  This  paper  attempts  to  make  the  scientific  community 
aware  of  a  number  of  bathymetric  artifacts  observed  in  Sea  Beam  data  which,  if 
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unrecognized,  might  lead  to  geological  misinterpretations.  We  have  shown  that 
artifacts  due  to  external  sound  sources  (e.g.  subbottom  profilers)  or  internal  sidelobe 
interference  can  usually  be  clearly  identified  as  resulting  from  spurious  data.  In  most 
of  these  cases,  corrective  action  can  be  taken  in  real-time  while  surveying.  We  also 
discussed  a  more  insidious  artifact  (the  ‘omega1  effect)  which  is  virtually  impossible  to 
detect  in  real-time  for  lack  of  warning.  In  addition,  such  artifacts  commonly  found 
when  steaming  downdip  over  slopes  greater  than  30  degrees  may  appear  as 
geologically  plausible  volcanic,  tectonic  or  sedimentary  features.  When  navigation  is 
based  on  seafloor  morphology,  failure  to  recognize  bathymetric  artifacts  may  lead  to 
positioning  errors. 

In  order  to  explain  the  causes  of  these  artifacts,  we  have  analyzed  Sea 
Beam’s  echo  detection  and  processing  techniques.  Errors  have  been  found  to  relate  to 
the  methods  of  sidelobe  rejection,  automatic  bottom-tracking  and  automatic  receiver 
gain  calibration.  These  errors  result  in  incorrect  depth  determinations  which  cause 
the  artifacts  observed.  Because  Sea  Beam  only  retains  depths  and  cross-track 
horizontal  distances  from  the  received  acoustic  signals,  investigators  have  no 
alternative  but  to  disregard  the  bathymetric  artifacts  they  identify.  A  number  of 
corrections  are  proposed  to  prevent  such  data  disposal:  improved  sidelobe  control  in 
the  transmit/ receive  acoustic  geometry,  extension  of  the  EP  receivers1  dynamic  range, 
sidelobe  and  noise  burst  rejection  through  advanced  adaptive  filtering  techniques, 
improved  bottom-tracking  gate  operation,  delayed  receivers1  gain  calibration  to  allow 
for  completion  of  the  reception  cycle  in  progress,  and  computer  coordinated  signal 
transmission  for  all  active  sound  sources  during  a  survey. 

Recently,  presentation  formats  for  Sea  Beam  data  have  extended  beyond 
contour  maps  to  include  gray-tone  and  color  shaded  relief  maps  [Edwards,  et  ah, 
1984],  and  bathymetric  gradient  charts  [Hey  et  al.,  1985a].  These  formats  are  very 
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valuable  in  interpreting  the  data,  however  ‘omegas’  and  other  bathymetric  artifacts 
will  persist,  because  the  errors  are  in  the  raw  Sea  Beam  data,  not  in  the  contouring 
algorithm  employed. 
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1.8.  APPENDIX  A:  Sea  Beam  Acoustic  geometry 

In  the  following,  the  beam  widths  are  calculated  at  the  half  power  point  of 
the  beam  patterns.  The  transmitted  beam  pattern  spans  54  degrees  athwartships  by 
2  2/3  degrees  in  the  fore-aft  direction.  It  is  pitch  stabilized  to  ensure  vertical 
projection  by  phasing  the  outputs  of  the  twenty  power  amplifiers  relative  to  a  pitch 
angle  supplied  by  the  vertical  reference  gyroscope  (Fig.  1.1)  within  the  limits  of  ±  10 
degrees  of  pitch.  As  shown  in  the  computed  beam  pattern  (Fig.  1  2a)  the  projector 
array  is  designed  for  sidelobe  attenuation  30  dB  down  from  the  mainlobe  and  grating 
lobes  appear  at  55  degrees  on  the  fore-aft  axis.  The  sidelobe  level  is  controlled  by 
amplitude  shading  the  output  of  the  twenty  power  amplifiers  (Fig.  1.1)  using  the 
Dolph-Chebyshev  amplitude  shading  method  for  acoustic  arrays  [Dolph,  1946,  Riblet 
and  Dolph,  1947  .  Since  the  array  is  contained  in  a  housing,  the  actual  sidelobe  level 
is  25  to  26  dB  down  from  the  mainlobe  [Dolph  1946;  Renard  and  Allenou  1979]. 
Proper  control  of  the  sidelobes  on  the  transmitted  beam  in  the  fore-aft  direction  is 
crucial  for  adequate  performance  of  the  system  when  the  ship’s  track  runs  downdip 
(across  a  slope  in  the  downhill  direction).  In  this  geometry,  weakly  attenuated 
sidelobes  ensonify  the  slope  at  near-normal  incidence  in  the  fore-aft  direction.  The 
corresponding  bottom  returns  are  received  earlier  than  those  due  to  vertical 
projection  in  the  mainlobe,  and  they  disrupt  the  echo  processing  and  bottom-tracking 
functions. 

The  design  of  the  receiving  array  yields  a  beam  pattern  which  is  2  2/3 
degrees  athwartships  by  20  degrees  in  the  fore-aft  direction  (Fig.  1.2b).  The  20 
degrees  beam  width  in  the  fore-aft  direction  is  meant  to  accommodate  pitch  angles  of 
±  10  degrees  as  no  pitch  stabilization  is  performed  on  the  receiving  beams.  Sea  Beam 
generates  sixteen  preformed  beams  fixed  with  respect  to  the  ship’s  vertical  by 
electronically  steering  such  2  2/3  degrees  beams  at  intervals  of  2  2/3  degrees 
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athwartships  from  20  degrees  incidence  on  port  to  20  degrees  on  starboard.  Dolph- 
Chebyshev  amplitude  shading  of  the  output  of  the  forty  preamplifiers  (Fig.  1.1) 
attenuates  the  side  lobes  30  dB  below  the  mainlobe  (Fig.  1.2b).  For  the  same  reasons 
given  for  the  projector  array,  the  actual  sidelobe  level  may  be  somewhat  higher. 
Renard  and  Allenou  [1979]  measured  a  value  of  28  dB  on  two  preformed  beams.  The 
acoustic  data  we  have  recorded  indicates  a  mean  sidelobe  level  of  25  dB  below  the 
mainlobe  on  ten  preformed  beams  for  the  SIO  system  [de  Moustier,  1985b],  Proper 
sidelobe  level  control  is  important  for  the  receiving  array,  because  each  of  the 
preformed  beams  has  sidelobes  oriented  in  the  direction  of  the  mainlobe  of  all  the 
other  beams.  A  strong  return  coming  into  the  mainlobe  of  a  particular  beam  will 
therefore  be  received  by  all  the  sidelobes  pointing  in  the  same  direction. 

1.9.  APPENDIX  B:  Sea  Beam  Echo  processing 

Information  concerning  Sea  Beam’s  echo  processing  is  contained  in  the  Sea 
Beam  software  technical  manual  [General  Instrument  Corp.,  1981].  In  the  following 
we  give,  with  the  manufacturer’s  permission,  an  overview  of  the  main  features  of  the 
echo  processing  software.  We  emphasize  the  features  important  to  understand  the 
causes  of  the  bathymetric  artifacts  discussed  in  this  paper. 

During  each  transmission  cycle,  16  bottom  returns  (e.g.  Fig.  1.3)  are  digitized 
by  Sea  Beam  s  analog-to-digital  converter  at  a  frequency  of  300  Hz  for  each  beam. 
For  each  conversion  cycle  the  Sea  Beam  computer  performs  the  following  operations: 
receiver  gain  correction,  refraction  correction,  roll  compensation,  detection  threshold 
level  determination,  and  signal  detection  for  each  of  the  roll  compensated  beams. 

The  receiver  gam  correction  consists  of  multiplying  the  digitized  signal 
voltage  for  each  beam  by  an  amplitude  multiplication  factor  to  compensate  for 
differences  among  receivers.  Because  the  roll  compensation  involves  interpolation 
between  beams,  it  is  important  that  all  16  signals  have  a  common  gain  at  any  one 
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time.  The  gains  of  the  individual  receivers  are  automatically  calibrated  by  the  EP 
every  half  hour  by  inputting  a  common  voltage  through  the  beam  line  drivers  (Fig. 
1.1),  and  digitizing  the  output  of  the  receivers. 

The  refraction  correction  uses  values  of  a  sound  velocity  versus  water  depth 
profile,  entered  at  the  beginning  of  a  survey  by  the  Sea  Beam  operator,  and  Snell’s  law 
to  calculate  the  reception  angle  0  for  each  beam  with  respect  to  the  ship’s  vertical. 
The  sound  velocity  profile  is  measured  with  an  expendable  bathythermograph  cast  for 
the  first  few  hundred  meters  and  extended  to  the  maximum  bottom  depth  in  the 
survey  area  using  values  from  Carter’s  tables  of  sound  velocity  in  the  ocean  [Carter, 
1980]. 

The  roll  compensation  uses  the  ship’s  instantaneous  roll  angle  given  by  the 
vertical  reference  gyroscope  to  reference  the  reception  angle  ©  to  the  true  vertical: 
4>  =©+/?.  A  set  of  stabilized  beam  angles  spaced  2  2/3  degrees  apart  are  then 
created.  The  amplitudes  of  the  stabilized  beams  are  linearly  interpolated  between 
those  of  the  two  adjacent  preformed  beams  with  corrected  reception  angle  4>,-  and  4>t*+1. 
This  yields  fifteen  stabilized  beams  each  2  2/3  degrees  wide,  fixed  in  a  vertical  plane 
athwartships  with  one  beam  aligned  with  the  true  vertical.  As  provision  has  been 
made  for  ±  20  degrees  of  ship’s  roll,  there  are  thirty  one  possible  stabilized  beam 
positions  between  ±  40  degrees.  Occasionally,  one  of  the  preformed  beam  angles  4> 
will  lie  on  the  true  vertical  ('I'  =0)  and  there  may  be  sixteen  stablized  beams. 

A  set  of  bottom-tracking  gates  determines  the  detection  time  window  during 
which  a  bottom  echo  is  expected  based  on  previous  sounding  history.  The  tracking 
gate  is  an  essential  feature  of  the  EP  because  it  conditions  proper  echo  signal 
detection  and  therefore  reliable  depth  determination.  Each  beam  has  its  own  tracking 
gate.  It  is  centered  on  the  average  depth  for  that  beam  using  depth  history  over  the 
last  five  transmission  cycles  (pings)  weighted  decreasingly  into  the  past.  The  gate 
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width  is  determined  by  the  observed  ping-to-ping  depth  fluctuations  with  allowance 
for  variations  in  signal  duration  due  to  beam  angle,  bottom  slope  and  beam  width.  As 
a  result,  the  gates  are  narrower  for  the  near-vertical  beams  than  for  the  outer  beams. 

A  constant  value  (20  m)  is  added  to  the  width  of  each  gate  as  a  safety  margin  to 
ensure  that  the  echo  signal  does  not  fall  outside  the  gate.  Bottom  echoes  falling 
outside  the  tracking  gates  are  not  taken  into  account  by  the  EP  which  usually  will  not 
compute  a  depth  and  a  cross-track  distance  for  the  corresponding  beams  for  lack  of 
signal  to  noise  ratio.  This  situation  creates  a  data  gap.  Since  only  fifteen 
(occasionally  sixteen)  of  the  possible  thirty  one  stabilized  beams  bear  data,  gate 
settings  for  null  or  unused  beams  are  interpolated  or  extrapolated  from  those  of 
adjacent  beams.  Finally  the  gate  settings  are  smoothed  across  all  thirty  one  possible 
beams.  The  analog-to-digital  conversion  starts  at  the  onset  of  the  gate  with  the 
shallowest  setting  (earliest  time).  The  conversion  stops  when  the  deepest  gate  has 
been  reached. 

The  detection  threshold  level  determination  is  a  very  critical  operation  in  the 
echo  processing.  It  is  adjusted  every  conversion  cycle,  and  is  therefore  a  dynamic 
process  taking  several  parameters  into  consideration:  1)  the  manual  threshold  level 
input  by  the  Sea  Beam  operator,  2)  the  background  noise  level  of  the  receivers,  3)  the 
receivers’  sidelobe  response  and  4)  potential  noise  bursts  interfering  with  the  bottom 
echo  detection.  In  general,  the  threshold  level  is  computed  to  ride  above  the  noise  and 
above  the  sidelobe  response.  For  reference,  the  noise  level  measured  on  data  similar 
to  that  of  Figure  1.3  is  usually  around  20  mv.  The  sidelobe  threshold  is  computed  as 
l/4t^1  the  amplitude  (12  dB  down)  of  the  highest  of  the  sixteen  signals  at  any  one  time 
(Fig.  1.3).  A  noise  burst  appears  as  a  synchronous  ridge  similar  to  the  sidelobe  ridge 
of  Figure  1.3,  but  the  amplitudes  of  the  individual  peaks  are  more  or  less  constant  on 
all  beams  (e  g.  Fig.  1.7).  By  comparing  the  maximum  amplitude  with  the  median 
amplitude  across  all  beams  at  any  one  time,  the  software  is  able  to  recognize  a  noise 
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burst.  When  a  noise  burst  is  detected,  and  when  the  corresponding  threshold  level  is 
higher  than  both  the  noise  and  the  sidelobe  thresholds,  the  sixteen  amplitudes  are 
rejected.  Otherwise  the  higher  of  the  noise  or  the  sidelobe  thresholds  will  be  used  as 
the  detection  threshold.  With  this  method  however,  canceling  sidelobe  response  or 
noise  bursts  when  they  overlap  with  a  bottom  return  results  in  cancellation  of  the 
corresponding  part  of  the  bottom  return.  Also,  because  of  saturation  in  the  EP 
receivers’  amplifiers,  sidelobe  rejection  is  only  partially  achieved  in  cases  when  the 
specular  return  is  clipped.  This  results  in  both  echo  detection  and  depth  computation 
errors.  Finally  for  each  conversion  cycle,  a  signal  sample  is  detected  if  it  is  above  the 
detection  threshold  and  within  the  bottom-tracking  gates. 

Once  the  analog-to-digital  conversion  sequence  has  been  completed  on  all 
beams  the  next  set  of  echo  processing  operations  is  done  once  per  transmission  cycle. 
The  signal  level  of  each  detected  beam  is  integrated  over  the  duration  of  the  detected 
return  (within  the  gates  and  above  the  threshold).  If  the  resulting  energy  in  the 
return  is  below  a  prescribed  minimum,  the  beam  is  deemed  invalid  due  to  poor  signal 
to  noise  ratio,  [Farr,  1980j .  For  a  valid  beam,  a  slant  range  is  calculated  by 
computing  the  center  of  mass  of  all  the  detected  signal  samples  for  that  beam,  and  by 
multiplying  the  corresponding  arrival  time  by  750  m/s.  Depth  and  cross-track 
horizontal  distances  are  then  calculated  as  described  in  Section  1.3  of  the  text. 


47 


1.10.  REFERENCES 


Ballard,  R.  D.  and  T.  H.  vanAndel,  Morphology  and  tectonics  of  the  inner  rift  valley 
at  lat.  36°50N  on  the  Mid-Atlantic  Ridge,  88,  507-580,  1977. 

Ballard,  R.  D.  and  J.  Francheteau,  Geologic  processes  of  the  mid-ocean  ridge  and  their 
relation  to  sulfide  deposition,  in  Hydrothermal  processes  at  sea  floor  spreading 
centers,  pp.  17-26,  Plenum  Press,  New  York,  1983. 

Carter,  D.  J.  T.,  Echo-sounding  correction  tables,  Hydrographic  Department  Ministry 
of  Defense,  Taunton,  1980. 

Clay,  C.  S.  and  H.  Medwin,  Acoustical  oceanography:  Principles  and  applications, 
p.  128,  John  Wiley  and  Sons,  1977. 

Crane,  K.,  Structure  and  tectonics  of  the  Galapagos  Inner  Rift,  86  0  10'W,  J.  Geology, 
86,  715-730,1978. 

Crane,  K.  and  R.  D.  Ballard,  The  Galapagos  Rift  at  86  °  W’,  Morphological  waveforms, 
4,  structure  and  morphology  of  hydrothermal  fields,  J.  Geophys .  Res.,  85,  1443-1454, 
1980. 

Crane,  K.,  F.  Aikman,  R.  Embley,  S.  Hammond,  A.  Malahoff,  and  J.  Lupton,  The 
distribution  of  geothermal  fields  on  the  Juan  de  Fuca  Ridge,  J.  Geophys.  Res.,  90, 
727-744,  1985. 

Detrick,  R.  S.,  P.  J.  Fox,  K.  Hastens,  W.  B.  F.  Ryan,  L.  Mayer,  J.  Karson,  and  A  Sea 
Beam  survey  of  the  Kane  Fracture  Zone  and  the  adjacent  Mid- Atlantic  Ridge  Rift 
Valley,  Trans.  Am.  Geophys.  Union,  65,  1006,  1984. 

Dolph,  C.  L.,  A  current  distribution  of  broadside  arrays  which  optimizes  the 
relationship  between  beam  width  and  side-lobe  level,  Proc.  Inst.  Radio  Engrs.,  8 4, 
335-348,  1946. 

Edwards,  M.  H  ,  R.  E.  Arvidson,  and  E.  A.  Guiness,  Digital  image  processing  of  Sea 
Beam  bathymetric  data  for  structural  studies  of  seamounts  near  the  East  Pacific 
Rise,  J.  Geophys.  Res.,  89,  11108-11116,  1984. 

Farr,  H.  K.,  Multibeam  bathymetric  sonar:  Sea  Beam  and  Hydrochart,  Marine 
Geodesy,  4,  77-93,  1980. 

Fornari,  D.  J.,  W.  B.  F,  Ryan,  and  P.  J.  Fox,  The  evolution  of  craters  and  calderas  on 
young  seamounts:  Insights  from  Sea  MARC  I  and  Sea  Beam  sonar  surveys  of  a  small 
seamount  group  near  the  axis  of  the  East  Pacific  Rise  at  10  °N,  J.  Geophys .  Res., 
89,  11069-11084,  1984. 

Francheteau,  J.  and  R.  D.  Ballard,  The  East  Pacific  Rise  near  21  0  N,  13  °  N  and  20  0  S; 
inferences  for  along  strike  variability  of  axial  processes,  Earth  and  Planetary  Science 
Letters,  64,  93-116,  1983. 


48 


Gallo,  D.  J.,  P.  J.  Fox,  and  J.  A.  Madsen,  The  morphotectonic  signature  of  fast- 
slipping  ridge-transform-ridge  systems:  A  synthesis  of  Sea  Beam  bathymetry, 
Trans.  Am.  Geophys.  Union,  65,  1103,  1984. 

General  Instrument  Corporation,  Sea  Beam  bathymetric  survey  system  technical 
manual,  volume  2,  1981 

Glenn,  M.  F.,  Introducing  an  operational  multi-beam  array  sonar,  International 
Hydrographic  Review,  XLVII( 1),  35-39,  Monaco,  1970 

Hey,  R.  N.,  M.  C.  Kleinrock,  S.  P.  Miller,  T.  M.  Atwater,  and  R.  C.  Searle,  Sea 
Beam/Deep-Tow  investigation  of  an  active  oceanic  propagating  rift  system,  J 
Geophys.  Res  (in  press),  1985a. 

Hey,  R.  N.,  D.  F.  Naar,  M.  C.  Kleinrock,  W.  J.  Phipps  Morgan,  E.  Morales,  and  J.  G. 
Schilling,  Microplate  tectonics  along  a  superfast  seafloor  spreading  system  near 
Easter  Island,  Nature  (in  press),  1985b. 

Kleinrock,  M.  C.,  R.  N.  Hey,  and  C.  de  Moustier,  The  ’Omega’  deception  in  Sea  Beam 
data,  Trans.  Am.  Geophys.  Union,  65,  1103,  1984. 

Lewis,  S.  D.,  J.  W.  Ladd,  T.  R.  Bruns,  D.  E.  Hayes,  and  R.  VonHuene,  Growth 
patterns  of  submarine  canyons  and  slope  basins,  Eastern  Aleutian  Trench,  Alaska, 
Trans.  Am.  Geophys  Union,  65,  1104,  1984. 

Lonsdale,  P.  F.,  Overlapping  Rift  Zones  at  the  5.5  °S  offset  of  the  East  Pacific  Rise,  J. 
Geophys.  Res.,  88,  9393-9406,  1983. 

Macdonald,  K.  C.  and  P.  J.  Fox,  Overlapping  spreading  centers:  A  new  kind  of 
accretion  geometry  on  the  East  Pacific  Rise,  Nature,  801,  55-58,  1983. 

Macdonald,  K.  C.,  J.  C.  Sempere,  and  P  J.  Fox,  East  Pacific  Rise  from  Siqueiros  to 
Orozco  fracture  zones:  along-strike  continuity  of  axial  neo-volcanic  zone  and 
structure  and  evolution  of  overlapping  spreading  centers,  J.  Geophys.  Res.,  89, 
6049-6069,  1984. 

Mammerickx,  J.,  Morphology  of  propagating  spreading  centers:  new  and  old,  J. 
Geophys.  Res.,  89,  1817-1828,  1984. 

McCool,  J.  M.  and  B.  Widrow,  Principles  and  applications  of  adaptive  filters:  A 
tutorial  review,  NOSC,  NUC  TP  580,  1977. 

de  Moustier,  C.,  Inference  of  manganese  nodule  coverage  from  Sea  Beam  acoustic 
backscattcring  data,  Geophysics  ,  50,  989-1001,  1985a. 

de  Moustier,  C.,  Beyond  Bathymetry,  mapping  acoustic  backscattering  from  the  deep 
sea  floor  with  Sea  Beam,  J  Acoust.  Soc.  Amer.,  submitted,  1985b. 

Naar,  D.  F.  and  R.  N.  Hey,  Fast  rift  propagation  along  the  East  Pacific  Rise  near 
Easter  Island,  J.  Geophys.  Res.  (in  press),  1985. 


49 


Patterson,  R.  B.,  Relationships  between  acoustic  backscatter  and  geological  charac¬ 
teristics  of  the  deep  ocean  floor,  ,/.  Acoust.  Soc.  Am.,  J6,  756-761,  1969. 

Phillips,  J.  D.  and  A.  S.  Fleming,  Multibeam  sonar  study  of  the  MAR  Rift  Valley  36- 
37  °  N,  Geol.  Soc.  of  Amer.,  1978.  Map  Series  MC-1 9. 

Renard,  V.  and  J.  P.  Allenou,  SEA  BEAM  multi-beam  echo-sounding  in  "Jean  Char¬ 
cot".  Description,  evaluation  and  first  results,  International  Hydrographic  Review, 
LVI(l),  35-67,  Monaco,  1979. 

Riblet,  II.  E.  and  C.  L.  Dolph,  Discussion  on  a  current  distribution  of  broadside  arrays 
which  optimizes  the  relationship  between  beam  width  and  side-lobe  level,  Proc.  Inst . 
Radio  Engrs.,  35,  489-492,  1947. 

Shipley,  T.  11.  and  G.  F.  Moore,  Sediment  accretion  and  subduction  in  the  Middle 
America  Trench,  OJI  International  Seminar  on  the  formation  of  ocean  margins,  1985. 

Smith,  S.  M.,  Sea  Beam  Operator  Manual,  SIO  Reference  83-7,  1983. 

Spiess,  F.  N  and  P.  F.  Lonsdale,  Deep- Tow  rise  crest  exploration  techniques,  Marine 
Tech.  J 16,  67-75,  1982. 

Spiess,  F.  N.,  R.  Messier,  G.  Wilson,  M.  Weydert,  and  P.  Rude,  Echo  I  cruise  report, 
SIO  Reference  84-3,  1984. 

Tyce,  R.  C.,  Sea  floor  mapping  aboard  R/V  CONRAD  -  The  first  year  of  Sea  Beam, 
Trans.  Am.  Geophys .  Union,  65,  1103,  1984. 

Urick,  R.  J.,  Principles  of  Underwater  Sound,  third  edition,  McGraw-Hill  Book  Co., 
1983. 

vanAndel,  Tj.  H.  and  R.  D.  Ballard,  The  Galapagos  Rift  at  86°  W:  2,  volcanism  struc¬ 
ture  and  evolution  of  the  Rift  valley,  J .  Geophys.  Res.,  84,  5390-5406,  1979. 


Chapter  II 


APPROACHES  TO  ACOUSTIC  BACKSCATTERING 
MEASUREMENTS  FROM  TIIE  DEEP  SEAFLOOR. 

H.l.  ABSTRACT 

Because  the  average  ocean  depth  is  four  kilometers,  seafloor  investigations 
are  mostly  remote  sensing  operations.  The  primary  means  to  determine  the  morphol¬ 
ogy,  the  structure,  and  the  texture  of  the  seafloor  are  acoustic.  This  paper  considers 
the  current  seafloor  remote  sensing  approaches  involving  acoustic  backscattering.  The 
physical  constraints  imposed  by  the  ocean  as  a  propagation  medium,  by  the  seafloor 
as  a  backscattering  boundary,  and  by  the  measuring  instruments  are  briefly  reviewed. 
The  sonar  systems  currently  used  by  the  oceanographic  community  for  deep  seafloor 
acoustic  backscattering  measurements  deal  with  these  constraints  differently  depend¬ 
ing  on  their  specific  application  and  on  whether  they  are  towed  behind  a  ship  or 
mounted  on  her  hull 

Towed  sidescan  systems  such  as  Gloria  II  (U.K.),  the  Sea  Mapping  and 
Remote  Characterization  (Sea  MARC)  1  and  II,  the  Deep  Tow  system  of  the  Marine 
Physical  Laboratory  (MPL),  and  hull-mounted  systems  such  as  Swathmap  all  give  a 
qualitative  measure  of  backscattering  by  converting  echo  amplitudes  to  gray  levels  to 
produce  a  sidescan  image  of  the  seafloor,  A  new  approach  is  presented  which  uses  a 
Sea  Beam  multibeam  echo-sounder  to  produce  similar  sidescan-like  images. 

Quantitative  measurements  of  backscattering  have  been  attempted  in  recent 
experiments  using  the  Deep  Tow  system  and  Sea  Beam.  Such  measurements  provide 
some  insight  into  the  geological  processes  responsible  for  the  acoustic  backscatter, 
with  useful  applications  for  geologists  as  well  as  designers  and  operators  of  bottom¬ 
interacting  sonars. 
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II.2.  INTRODUCTION 


With  an  average  depth  of  four  kilometers,  the  ocean  floor  is  only  accessible 
by  special  purpose  submersibles  from  which  observations  are  limited  by  the  aperture 
of  a  view  port.  For  this  reason,  most  seafloor  investigations  arc  remote  sensing  opera¬ 
tions  which  use  underwater  sound  as  their  primary  tool  and  core  or  grab  sample  as 
well  as  bottom  photography  as  ground  truth.  Seafloor  acoustic  measurements  are  com¬ 
monly  divided  in  two  broad  categories:  low  frequency  (<100  II/)  seismic  measurements' 
in  reflection  or  refraction  work,  and  high  frequency  (>  3  kHz)  reflectivity  or  back- 
scattering  measurements.  This  paper  focuses  on  the  high  frequency  seafloor  acoustic 
backseatter  which  is  on  one  hand  a  noise  background  against  which  active  sonars 
must  operate,  and  on  the  other  hand  a  source  of  information  for  geologists  because  its 
variations  are  caused  by  changes  in  bottom  type  or  bottom  microroughness. 

The  systems  currently  used  by  the  oceanographic  community  to  measure 
acoustic  backseatter  from  the  deep  seafloor  are  most  often  sidescan  sonars  or  echo- 
sounders,  and  occasionally  specially  designed  multifrequency  arrays.  Sidescan  sonars 
give  a  qualitative  measure  of  backscattering  by  converting  echo  amplitudes  to  gray 
levels  in  the  process  of  forming  an  acoustic  image  of  the  seafloor,  and  their  resolution 
depends  on  a  combination  of  parameters  such  as  frequency,  pulse  length,  beamwidth  ... 
[Somers  and  Stubbs,  1984].  The  same  parameters  apply  to  echo-sounders  whose  pri¬ 
mary  output  is  numerical  bathymetry.  However,  broad-beam  (>30°)  echo-sounders 
have  been  used  to  obtain  quantitative  normal-incidence  measurements  in  an  effort  to 
classify  bottom  types  [Breslau,  1957]  or  bottom  microroughness  [Stanton,  1984).  Based 
on  the  author’s  work,  this  paper  shows  that  multi-narrowbeam  echo-sounders  such  as 
Sea  Beam  [Farr,  1980;  Renard  and  Allenou,  1979]  provide  a  means  to  obtain  both 
quantitative  seafloor  acoustic  backscattering  measurements  at  discrete  angles  of 
incidence  and  qualitative  measures  in  an  acoustic  imaging  mode  which  is  new  to  such 
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systems. 

In  the  following,  the  physical  constraints  which  condition  the  performance  of 
sonar  systems  are  reviewed  briefly.  The  approaches  taken  in  the  design  of  several 
operational  systems  are  compared,  and  recent  experiments  involving  rnultifrequency 
arrays  are  presented. 

IL3.  SONAR  DESIGN  CONSTRAINTS 

This  section  briefly  reviews  the  constraints  imposed  by  the  ocean  as  a  pro¬ 
pagating  medium,  by  the  seafloor  as  a  backscattering  boundary  and  by  the  physical 
characteristics  of  the  sonar  on  the  design  of  systems  to  measure  seafloor  acoustic 
backscatter.  More  exhaustive  treatments  on  the  subject  may  be  found  in  Urick  1983 
and  in  Clay  and  Medwin  (1977]. 

As  sound  waves  propagate  through  the  water  column,  they  are  attenuated 
due  to  spherical  spreading  and  absorption.  The  former  increases  as  the  square  of  range 
and  the  latter  increases  roughly  as  the  square  of  frequency.  This  frequency  dependence 
is  usually  imbedded  in  a  logarithmic  absorption  coefficient  a  in  decibels  per  meter  (e.g. 
a  ^lCT5  dB/m  at  10  kHz)  which  sets  limits  on  the  operating  frequency  of  a  sonar  for  a 
given  range.  The  deep-water  ambient-noise  level  is  also  frequency  dependent  and 
decreases  5  to  6  dB  per  octave  in  the  interval  1-100  kHz.  At  long  ranges,  attenuation 
dominates  and  adversely  counterbalances  improvements  in  ambient-noise  characteris¬ 
tics  obtained  by  using  higher  frequencies.  Expressed  in  decibels,  the  attenuation  terms 
add  into  a  transmission  loss  : 

TL  =20 logl0  r  +  ar  (1) 

where  the  slant  range  r  is  in  meters.  The  propagation  medium  also  imposes  horizontal 

range  limitations  on  sound  waves  as  they  are  refracted  due  to  variations  in  the  sound 

velocity  with  temperature  and  pressure.  In  practice,  the  range  limitation  due  to  ray- 


bending  depends  on  water  depth  and  bottom  relief.  At  oceanic  depths  (>4  km)  the 
maximum  horizontal  range  attainable  with  a  sonar  at  the  sea  surface  is  between  30 
and  40  km  [Andrews  and  Humphrey,  1980 

Bottom  relief  enters  as  a  geometrical  parameter  in  the  range  limitation,  but 
for  a  given  relief,  bottom  type  and  bottom  microroughness  condition  the  backscatter- 
ing  process.  Whether  a  body  backscatters  sound  effectively  depends  on  how  its  density 
and  compressibility  differ  from  those  of  the  surrounding  ocean,  and  how  its  roughness 
scale  compares  to  the  acoustic  wavelength.  There  is  no  simple  theory  to  predict  the 
level  of  seafloor-backscattered  sound  waves,  and  one  relies  on  reported  measurements 
(mostly  done  in  shallow  water  [e.g.  Wong  and  Chesfcrrnan,  1968;  McKinney  and 
Anderson,  1964])  and  geoacoustic  models  Hamilton,  1980  to  derive  empirical  expres¬ 
sions  for  bottom  backscattering  strength  as  a  function  of  grazing  angle  and/or  fre¬ 
quency  for  various  types  of  substrates. 

In  practice,  the  bottom  backscattering  strength  BS,  in  decibels,  is  expressed 
as  the  sum  of  a  backscattering  strength  per  unit  area  (lm2)  SB  and  an  effective 
scattering  area  A  : 

BS  =  SB  +  A  with  SB  =  50  +  10/o</10  sin^a  and  A  =  H)log  10( •  )  (2) 

2  cost* 

where  as  indicated  in  Figure  2.1,  r  is  the  slant  range,  0  the  horizontal  beam  width  of 
the  transducer,  t*  the  grazing  angle,  r  is  the  transmitted  pulse  length,  c=1500  m/s  the 
sound  speed,  and  S0  is  a  scattering  coefficient  independent  of  grazing  angle.  For  most 
sidescan  sonar  or  echo-sounder  applications,  a  scpiare  law  dependence  on  siim  ( K  =2) 
is  a  reasonable  match  to  existing  shallow-water  measurements  and  by  extrapolation  to 
the  few  measurements  existing  for  the  deep  seafloor.  The  dependence  of  backscatter¬ 
ing  on  bottom  type  is  taken  into  account  in  the  scattering  coefficient  S0  which 
increases  by  nearly  25  dB  from  clay,  through  silt  and  sand  to  rock;  and  bottom  rough¬ 
ness  can  cause  variations  of  several  dB  in  S0  for  the  same  bottom  type. 
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Figure  2.1  Sidescan  geometry .  A  transducer  array  at  depth  D  and  elevation  H  above  the 
seafloor  has  a  beam  pattern  with  angular  dimensions  0  (horizontal)  and  (vertical).  It 
ensonifies  a  strip  of  seafloor  out  to  horizontal  range  R  at  right  angle  to  the  direction  of  travel 
(arrow).  A  portion  of  the  strip  (broken  line)  is  not  ensonified  as  it  is  shadowed  by  the  hill. 
Backscattered  acoustic  intensities  received  at  the  transducer  are  converted  to  gray  levels  and 
mapped  at  increasing  times  of  arrival  (or  slant  range  r).  In  this  case,  the  first  return  is 
received  H/750  sec.  after  transmission.  On  the  acoustic  image,  specular  returns  appear  darker 
than  the  normal  background  and  shadow  zones  are  white.  A/?  is  the  cross-track  resolution 
and  o  is  the  grazing  angle. 


As  mentioned  above,  the  effects  of  bottom  roughness  depend  on  the  operating- 
characteristics  of  the  sonar.  To  sample  the  small  scale  roughness  rather  than  bottom 
slopes  requires  a  system  with  narrow  horizontal  beamwidth  and  large  bandwidth. 
These  parameters  are  constrained  by  the  physical  characteristics  of  the  transducer 
used. 

The  bandwidth  capability  of  the  transducer  is  usually  limited  to  10-15%  of 
the  center  frequency  [Kosalos,  1983].  Since  the  range  resolution  of  a  system  with 
bandwidth  w  is  AR=c/2v>,  the  incentive  is  to  use  high  frequencies  for  greater- 
bandwidth  and  higher  range  resolution,  but  attenuation  imposes  horizontal  range  limi¬ 
tations  and  it  becomes  necessary  to  operate  close  to  the  bottom.  The  most  common 
and  simplest  way  of  achieving  range  resolution  in  sonars  is  to  use  short  pulsed  con¬ 
tinuous  wave  signals  (CW)  but  the  shortest  available  pulse  length  r  is  limited  to 
r= l/u>.  An  alternative  approach  is  to  use  long  frequency-modulated  (FM)  pulses  with 
ivt>1  and  process  the  returns  through  matched  filters.  The  advantages  of  this  method 
are  the  possibility  to  use  lower  frequency  and  still  maintain  adequate  range  resolution, 
and  a  theoretical  signal-to-noise  power  gain  over  short  CW  pulses  with  equal  source 
level  of  10 logl0(wr)  [Somers  and  Stubbs,  1984].  However  implementation  of  the  correla¬ 
tion  processor  adds  complexity  and  cost  to  the  overall  sonar  system  design. 

The  beamwidth  requirements  are  constrained  by  the  size  of  the  transducer 
array  and  by  its  operating  frequency.  Defining  the  beamwidth  as  the  width  of  the 
mainlobe  of  the  radiation  pattern  3  dB  down  from  its  on-axis  response,  the  beamwidth 
9  in  radians  is  roughly  the  reciprocal  of  the  number  of  wavelength  across  the  effective 
aperture  of  the  array  i.e.  0=X/L  with  X=c//.  L  is  the  length  of  the  aperture  in  meters, 
X  is  the  acoustic  wavelength  in  meters  and  f  is  the  acoustic  frequency  in  Hertz.  Conse¬ 
quently,  for  a  given  frequency,  the  longer  the  array  the  narrower  the  beamwidth. 
However  the  optimum  size  of  a  transducer  is  a  compromise  between  its  directionality 


and  its  acoustic-power  output  capability. 

The  maximum  strain  bearable  by  the  transducer  and  cavitation  limit  the 
maximum  power  to  which  a  transducer  can  be  excited.  These  constraints  set  a  lower 
limit  on  the  size  of  the  effective  aperture  necessary  to  keep  the  radiated  power  per 
unit  area  below  the  cavitation  threshold.  The  power  P  per  unit  area  necessary  to  pro¬ 
duce  cavitation  increases  as  the  square  of  the  ambient  pressure  Pa.  P~P%/2pc 
watts  /cm2  where  pc=  1.5  105£  is  the  acoustic  impedance  of  water  [Albers,  1965]. 

A  substantial  gain  in  cavitation  threshold  is  therefore  achieved  by  operating  an  array 
a  few  decameters  below  the  surface  with  the  added  advantage  of  reduced  refraction 
effects  as  one  goes  below  the  thermocline.  This  technique  is  used  in  most  shallow- 
towed  sidescan  arrays. 

In  the  case  of  sidescan  sonar  arrays  mounted  on  either  side  of  a  towed  vehicle 
for  port-starboard  coverage,  an  additional  constraint  is  imposed  on  the  choice  of 
operating  frequency  by  the  existence  of  mutual  interference  (cross-talk)  between  the 
arrays  due  to  radiation  of  sound  from  the  back  of  each  array.  As  a  result,  a  mirror 
image  of  returns  on  one  side  is  mapped  on  the  other  side  (for  examples,  see  Belderson 
et  al.  [1972]).  In  general,  this  problem  is  alleviated  by  using  slightly  different  operating 
frequencies  on  each  side. 

Finally,  an  estimate  of  the  echo-to-noise  ratio  EN  measurable  with  a  given 
sonar  system  is  obtained  through  the  sonar  equation  which  combines  the  various 
aforementioned  design  parameters.  In  decibel  units,  this  equation  has  the  form: 

EN  =SL  -2TL  -  (AL  +  10 loglQw)  +  BS  +20loglob{$)  (3) 

where  SL  is  the  source  level,  2TL  accounts  for  the  round-trip  transmission  loss  (Eq.  1) 

between  the  sonar  and  the  bottom,  (NL-\-10loglow)  is  the  noise  level  in  the  bandwidth 

w ,  BS  is  the  backscattering  strength  (Eq.  2)  and  201og106(<t>)  accounts  for  the 


transducer’s  vertical  directivity. 


57 


For  most  sidcscan  systems  the  amplitude  corresponding  to  the  echo  level  is 
converted  to  a  gray  level  and  output  on  a  linescan  recorder  to  create  an  acoustic 
image  of  the  seafloor  surveyed.  In  the  process,  it  is  customary  to  multiply  the  ampli¬ 
tudes  by  a  time  varied  gain  to  compensate  for  transmission  losses  and  keep  the 
returns  within  the  dynamic  range  of  the  recording  instrument. 

IL4.  CURRENT  SYSTEMS 

The  design  constraints  outlined  above  are  dealt  with  differently  depending  on 
the  intended  application  of  the  sonar  system.  In  this  section,  several  systems  currently 
used  by  the  oceanographic  community  are  compared  on  the  basis  of  their  design 
approach,  and. of  the  characteristics  of  the  acoustic  measurements  obtained. 


System 

Gloria 

Swathmap 

Sea  MARC  H 

Sea  MARC  I 

Sea  Beam 

Deep  Tow 

F requency  (kHz) 

P  6.6,  S  6.7 

3.5 

P  11,  S  12 

P  27,  S  30 

12.168 

110 

FM  sweep 

FM  sweep 

cw 

CW 

CW 

CW 

Pulse  Length  r 

2,  4  sec 

.26-10  ms 

.  15-3.2  ms 

7  ms 

.2,  .5  ms 

Pulse  rep.  (sec) 

20-40 

40-48 

1,  2,  4,  8,  16 

.5,  1,  2  or  4 

1-22 

1 

Bandwidth  w 

100  Hz 

2  kHz— 000  Hz 

5  kHz —200  Hz 

226  Hz 

10  kHz 

Beamwidth  (deg.)  9  ,  $ 

2.7  ,  30 

2  ,  40 

17  ,  60 

T  2.7  ,  54 

R  16(20  ,  2.7) 

.75  ,  60 

Cross-track  resolution  A R 

20  m 

5  m 

.6  m 

.15  m 

depth  D 

Tow 

altitude  H 

30-60  m 

hull  mount 

60-100  m 

up  to  6  km 

hull  mount 

up  to  7.5  km 

12-1250  m 

10-100  m 

Horizontal  range  (km)  R 

15-30 

up  to  36 

.5,  1,  2.6,  6 

.26,  .6,  1,  2.5 

3/8  water  depth 

.5 

Speed  (knots) 

up  to  11 

up  to  20 

up  to  10 

1-3 

up  to  15 

1-2 

length  (m) 

Array 

width  (m) 

6.33 

3.8 

1.6 

T  2.8,  R  2.8 

1.25 

.18 

.2 

.2 

T  .16,  R  .4 

.08 

length  (m)  ' 

Vehicle 

width  (m) 

7.75 

6.6 

3 

2 

.8 

1.3 

1.2 

.7 

Table  2.1  Operational  characteristics 


Six  systems  covering  a  wide  spectrum  of  ranges  in  deep-ocean  work  have  been 
chosen  for  this  comparison.  They  are  the  British  system  Gloria  II  [Somers  et  al,  1978; 
Laughton,  1981],  Swathmap  [Andrews  and  Humphrey,  1980a-b;  Andrews  et  al,  1977], 
the  Sea  Mapping  and  Remote  Characterization  systems  (Sea  MARC1  I  and  II) 
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[Blackington  et  al,  1983;  Ilussong  and  Fryer,  1983;  Kosalos  and  Chayes,  1983|,  Sea 
Beam  and  the  Deep  Tow  instrument  package  [Spiess  and  Lonsdale,  1982;  Spiess  and 
Tyce,  1973],  The  operating  characteristics  of  these  systems  are  listed  in  Table  2.1 
where  symbols  appearing  in  Section  II. 3  and  Figure  2.1  have  been  repeated  for  ease  of 
correspondence. 

Gloria  II,  Sea  MARC  I  and  II  and  Deep  Tow  are  original  sidescan  sonar  sys¬ 
tems  whereas  Sea  Beam  is  a  multibeam  echo-sounder  and  Swath  map  is  a  sidescan 
application  of  the  SQS-26  ASW  sonar  used  aboard  U.S.  Navy  frigates.  The  table  is 
incomplete  for  Swathmap  because  some  of  the  technical  details  of  the  system  are 
classified.  In  the  Swathmap  design,  a  beam  is  steered  to  one  side  of  the  ship’s  track 
only  with  a  maximum  horizontal  range  of  about  36  km  and  a  cross-track  resolution  of 
several  hundred  meters.  The  high  ship  speed  results  in  a  poor  along  track  resolution  as 
transmit  cycles  are  spaced  400  to  500  m  apart  at  20  knots  (~10m/s).  This  system  is 
primarily  a  reconnaissance  tool  designed  to  map  at  a  rapid  rate  the  large  scale 
seafloor  relief  such  as  seamounts  or  fracture  zones,  and  the  intermediate  relief  (50- 
100m)  typical  of  abyssal  hills.  Although  the  acoustic  backscatter  measured  with 
Swathmap  is  modulated  by  bottom  texture,  it  is  dominated  by  slope  effects,  and  the 
records  are  qualitative  representations  of  bottom  slopes  and  regional  trends. 

Similar  measurements,  with  greater  along  and  across-track  resolution,  arc 
obtained  with  Gloria  II  which  records  acoustic  backscattering  from  both  sides  of  a 
shallow-towed  vehicle.  Different  frequencies  are  used  in  the  port  and  starboard  arrays 
to  avoid  cross-talk  between  them  To  enhance  its  signal  to  noise  characteristics,  this 
system  uses  a  long  transmit  pulse  and  processes  the  echoes  by  match-filtering  tech¬ 
niques.  The  records  are  automatically  corrected  from  slant  range  to  horizontal  range 
by  projecting  slant  range  onto  a  horizontal  plane  at  the  mean  depth  below  the  vehicle 
[Searle,  personal  communication].  Measurements  are  qualitative  and  recorded  acoustic 


returns  are  usually  dominated  by  bottom  slopes  with  marginal  indication  of  textural 
changes  [Teleki  et  a].,  1981).  The  strength  of  this  systems  is  its  ability  to  map  with 
sufficient  detail  large  swaths  of  seafloor  (30  km  or  more)  at  an  average  tow  speed  of  8 
knots. 

The  Sea  MARC  11  system  is  also  a  shallow-tow  vehicle  with  port-starboard 
coverage  and  a  tow-speed  capability  comparable  to  that  of  Gloria  II.  It  operates  at 
higher  frequencies  and  uses  short  transmit  pulses  resulting  in  better  cross- track  resol u 
tion  but  reduced  horizontal  ranges  (Table  2.1).  In  order  to  maintain  even  spacing 
between  data  points  in  the  sidescan  image,  the  bottom  returns  are  sampled  more 
rapidly  at  close  ranges  than  at  far  ranges.  As  with  Gloria  11,  the  records  are  automat¬ 
ically  corrected  for  slant  range  by  assuming  a  fiat  plane  at  the  mean  depth  below  the 
vehicle.  Due  to  its  finer  resolution,  the  qualitative  bottom  backscattering  measure¬ 
ments  obtained  from  Sea  MARC  II  images  arc  useful  for  bottom  slope  and  texture 
determination  on  a  scale  of  acoustic  wavelengths  (-^12  cm),  and  for  large  scale 
regional  trend  determination.  Sea  MARC  II  uses  a  pair  of  transducer  arrays  (e.g.  Fig. 
2.1)  on  each  side  of  the  vehicle.  By  measuring  the  phase  difference  between  the  outputs 
of  the  two  arrays,  it  is  possible  to  determine  the  angle  of  arrival  of  a  given  return  and 
therefore  compute  the  corresponding  depth  and  cross-track  distance.  In  a  post¬ 
processing  operation,  approximately  one  hundred  such  pairs  of  depths  and  cross-track 
distances  are  computed  for  every  transmission  cycle  and  used  to  produce  a  contour 
map  of  the  swath  of  seafloor  surveyed.  This  bathymetry  is  an  important  element  in 
the  interpretation  of  backscattering  measurements  made  with  the  system,  because  it 
offers  the  potential  to  remove  bottom  slope  effects  from  the  data  while  in  principle 
retaining  the  effects  of  bottom  composition  and  microroughness. 

The  Sea  MARC  I  system  is  a  deep-towed  version  of  Sea  MARC  II  with  a 
lesser  swath  width  (5  km  maximum).  Its  use  of  higher  frequencies  and  short  pulse 


lengths  yields  sub-meter  cross-track  resolution;  and  its  inherent  slow  speed  over  the 
bottom  yields  an  along-track  resolution  between  .5  and  5  m,  depending  on  the  pulse 
repetition  rate  (Table  2.1).  Although  this  system  uses  the  same  transducer  array  pair 
configuration  as  Sea  MARC  II,  the  phase  measurement  technique  is  not  implemented. 
Without  numerical  bathymetry,  it  is  therefore  not  possible  to  correct  the  backscatter- 
mg  measurements  for  bottom  slope.  Nontheless,  the  fine  resolution  of  the  sonar  allows 
detection  of  changes  in  bottom  texture  over  areas  of  constant  slope  between  10  and 
100  m  in  extent  depending  on  the  swath  width  chosen. 

Through  a  combination  of  frequency,  pulse  length  and  tow  speed,  the  sides- 
can  sonars  of  the  Deep  Tow  instrument  package  achieve  the  highest  resolution  of  all 
the  systems  listed  in  Table  2.1.  The  lateral  coverage  to  port,  and  starboard  is 
inherently  limited  to  a  swath  about  1  km  wide.  This  system  is  primarily  intended  for 
fine  scale  studies  of  the  seafloor  down  to  depths  in  excess  of  7  km.  The  backscatter 
measurements  are  qualitative  gray  scale  displays  uncorrected  for  slant  range.  Changes 
in  bottom  texture  are  readily  observable  and  micro-relief  (<1  m)  such  as  small  fissures 
is  resolvable  on  such  displays.  However,  the  limited  coverage  of  the  sonars  make  the 
system  marginally  useful  for  regional  trend  assessment  as  it  requires  long  and  costly 
surveys.  The  strength  of  the  Deep  Tow  instrument  package  is  its  versatility  as 
environmental  (e.g.  temperature,  conductivity)  and  other  geophysical  (e.g.  magnetics, 
subbottom  profiles)  data  can  be  collected  simultaneously  with  the  sidescan  data;  and 
bottom  photographs  can  be  taken  on  the  same  lowering  [Spiess  and  Lonsdale,  1982], 
The  phase  measurement  technique  mentioned  above  is  feasible  but  not  implemented 
on  the  Deep  Tow  system. 

H.5.  NEW  APPROACHES 

The  systems  described  in  the  previous  section  all  operate  as  sidescan  sonars. 
As  such  they  provide  qualitative  acoustic  baekscattering  measurements  of  the  deep 
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Figure  2.2  Envelopes  of  echoes  received  by  Sea  Beam's  16  beams. 

Beams  are  equally  spaced  2  0  2/3  apart  within  ±20  ’  of  the  ship’s  vertical  axis.  In  this  display, 
returns  are  not  compensated  for  ship’s  roll.  Time  is  in  seconds  after  transmission.  Amplitudes 
have  been  corrected  for  transmission  loss  by  a  time-varied  gain  and  are  displayed  in  volts. 


seafloor.  This  section  looks  at  new  measurement  approaches  using  multibeam  echo- 
sounders  such  as  Sea  Beam  or  experimental  multifrequency  acoustic  arrays  mounted 
on  instrument  packages  such  as  Deep  Tow  |Spicss  et  al,  I98l|. 

Sidescan  sonars  typically  transmit  and  receive  with  the  same  fan-shaped 
beam  which  is  narrow  along  track  and  broad  across  track  (Figure  2.1).  The  Sea  Beam 
system  transmits  with  a  similar  geometry,  but  it  receives  with  sixteen  narrow  beams 
spaced  2  2/3  degrees  apart  athwartships  within  ±20  degrees  incidence  Each  receive 
beam  is  2  2/3  degree  wide.  With  this  geometry,  the  system  is  able  to  process  back- 
scattered  returns  (Fig  2.2)  at  discrete  angles  of  incidence  and  calculate  a  set  of  depths 
and  cross-track  distances  for  each  transmission  cycle.  This  numerical  bathymetry  is 
then  output  in  near  real-time  as  a  contoured  chart  of  the  swath  of  seafloor  surveyed. 

Sea  Beam  processes  the  acoustic  signals  to  determine  depths  but  has  no 
internal  provision  for  recording  the  actual  waveform  To  preserve  these  signals, 
requires  a  parallel  acoustic  data  acquisition  system  which  records  digitally  the  echo 
envelopes  on  magnetic  tape  for  later  processing  and  analysis  [de  Moustier,  1985a],  A 
typical  set  of  returns  corresponding  to  one  transmission  cycle  is  illustrated  in  Figure 
2.2.  Such  data  make  it  possible  to  obtain  quantitative  backscattering  measurements 
at  discrete  angles  of  incidence.  Because  variations  in  bottom  characteristics  cause 
fluctuations  in  the  acoustic  backscatter,  and  because  numerical  bathymetry  allows 
corrections  for  slope  effects  to  be  made,  the  geological  characteristics  of  the  seafloor 
(e.g.  microroughness  or  bottom  type)  can  be  inferred  in  part  from  these  measurements 
[de  Moustier,  1985b;  Patterson,  1967], 

In  a  new  application  of  the  Sea  Beam  system,  the  acoustic  returns  it  receives 
are  used  in  a  sidescan-like  mode  by  combining  echoes  on  either  side  of  vertical 
incidence.  Working  with  digitally  recorded  Sea  Beam  acoustic  data  (Fig.  2.2),  a  peak 
detection  process  is  used  to  obtain  echo  amplitudes  at  incremental  slant  ranges  to 
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Figure  2.3  Acoustic  images  and  bathymetry,  (a)  Slant-range  corrected  SeaMARC  I  image  of 
the  crest  of  the  East  Pacific  rise  at  10  °  N.  Distances  are  true  horizontal  distances,  (b) 
Uncorrected  Sea  Beam  image  of  the  same  area;  the  dashed  line  in  (a)  represents  the 
corresponding  Sea  Beam  track.  The  cross-track  dimension  is  in  seconds  from  first  arrival 
(differential  slant  range),  (c)  Corresponding  swath  of  Sea  Beam  bathymetry  at  10  m  contour 
interval.  Tick  marks  point  downhill. 
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ANGLE  OF  INCIDENCE  IN  8/3  DEGREE  STEPS 


Figure  2.4  Echo  character  of  individual  roll-compensated  beams .  Same  Sea  Beam  acoustic 
data  as  that  of  Figure  2.3.  At  each  angle  of  incidence,  vertical  grid  lines  mark  the  time  of  first 
arrival  for  each  transmission  cycle,  and  echoes  are  displayed  in  bins  .36  sec.  wide.  The  gap 
between  a  grid  line  and  the  corresponding  return  is  the  differential  slant  range.  Intensified 
returns  (darker)  observable  beyond  10°  incidence  are  indicative  of  slopping  relief  as  sketched 
in  Figure  2.1. 
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port  and  starboard  This  process  starts  at  the  first  arrival  of  each  transmission  cycle. 
An  example  of  the  resulting  acoustic  image  is  shown  in  Figure  2.3  along  with  the 
corresponding  Sea  Beam  bathymetry,  and  a  ScaMAIlC  I  image  of  the  same  region 
included  for  comparison.  In  this  example,  the  Sea  Beam  acoustic  image  is  not 
corrected  for  for  slant  range  or  for  angle  of  incidence.  However,  as  Sea  Beam  com¬ 
putes  a  cross-track  bathymetric  profile  every  transmission  cycle,  slant  range  correc¬ 
tions  can  be  performed  without  the  ambiguities  inherent  to  the  conventional  horizon¬ 
tal  plane  method  mentioned  previously.  The  correction  for  angle  of  incidence  entails 
applying  a  time-varying-gain  to  the  signals  displayed  in  order  to  compensate  for  the 
drop  in  backscatter  as  a  function  of  angle.  Such  a  correction  is  performed  in  most 
sidescan  systems  and  yields  an  acoustic  image  having  nearly  uniform  resolution  out  to 
the  edges. 

This  new  acoustic  imaging  application  of  Sea  Beam  gives  textural  informa¬ 
tion  about  the  seafloor  and  brings  out  features  not  discernible  in  the  contoured  bathy¬ 
metry.  Unlike  conventional  sidescan  images  (Fig  2.1),  acoustic  images  obtained  with 
Sea  Beam  do  not  contain  shadows  because  angles  of  incidence  are  limited  to  the  range 
±  20  degrees,  (note  that  most  sidescan  sonars  do  not  record  data  in  this  angular  sec¬ 
tor).  The  qualitative  measure  of  backscatter  given  by  these  images  is  therefore 
representative  of  seafloor  microroughness  and  bottom  slopes.  The  effect  of  bottom 
slopes  can  also  be  removed  from  the  backscatter  data  given  that  quantitative  slope 
information  is  available  in  Sea  Beam’s  bathymetry. 

The  advantage  of  Sea  Beam’s  discrete  narrow  beams  can  also  be  preserved 
by  outputting  the  echo  character  on  each  beam  in  a  gray  level  display  such  as  that  of 
Figure  2,4.  Although  this  display  is  more  difficult  to  interpret  than  conventional  sides¬ 
can  images,  it  gives  bottom  structural  details  which  are  smoothed  out  in  the  process¬ 
ing  used  to  create  the  acoustic  image  of  Figure  2.3b. 


As  mentioned  in  Section  II.3,  the  effects  of  the  seafloor  as  an  acoustic  back- 
scattering  boundary  are  difficult  to  predict,  and  one  relies  on  measured  data  to  derive 
empirical  backscattering  models  usable  by  designers  of  sonar  systems  [Weinberg,  1981 
However,  there  is  a  definite  need  for  deep  seafloor  measurements  as  the  existing  data 
base  is  very  limited.  Systems  such  as  Sea  Beam  give  access  to  quantitative  measure¬ 
ments  as  a  function  of  angle  of  incidence  within  the  limits  of  the  available  bearriwidth. 
Recent  experiments  to  obtain  quantitative  measurements  as  a  function  of  frequency 
have  also  been  carried  out  using  multifrcquency  arrays  (4.5,  9,  15,  28,  GO,  110,  1  GO 
kHz)  mounted  on  the  stern  of  the  Deep  Tow  instrument  package,  and  projecting  fan 
beams  pointing  aft  [Spiess  et  ah, 1984b 

In  general  single  frequency  measurements  give  non-unique  answers  about  the 
nature  of  the  seafloor  because  there  is  an  infinite  combination  of  roughness  and  bot¬ 
tom  type  influencing  the  backscattering  process,  and  measurements  need  to  be  vali¬ 
dated  by  direct  bottom  sampling  or  by  bottom  photographs  or  television.  On  the  other 
hand,  multiple  frequency  systems  allow  simultaneous  sampling  of  several  roughness 
scales  thereby  giving  some  indication  of  the  respective  contribution  of  roughness  and 
bottom  type  in  the  backscattering  process.  Such  measurements  may  eventually  help 
geologists  to  acoustically  differentiate  between  bottom  types,  or  determine  the  size  of 
bottom  microrelief  (e.g.  furrows,  manganese  nodules,  etc.).  Similar  multifrequency 
measurements  have  been  obtained  with  conventional  hull-mounted  echo-sounders  in  an 
effort  to  determine  the  size  of  polymetallic  nodules  [Sumitomo  Metal  Mining  Co., 
1981]. 

H.6.  CONCLUSIONS 

In  the  remote  sensing  of  the  deep  seafloor,  acoustic  backscattering  can  be 
measured  with  sidescan  sonars  or  echo-sounders.  Sidescan  sonars  usually  provide  qual¬ 
itative  measurements  of  backscattering  which  allow  geologists  to  make  regional 
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structural  analysis  and  textural  assessments  of  the  seafloor,  with  a  resolution  depend¬ 
ing  on  the  operating  characteristics  of  the  sonar.  A  new  application  of  the  Sea  Beam 
system  in  a  sidescan  mode  has  been  shown  to  yield  similar  qualitative  measurements 
whose  usefulness  is  enhanced  by  the  precision  numerical  bathymetry  normally  avail¬ 
able  with  this  system  The  Sea  Beam  also  offers  the  potential  to  obtain  quantitative 
backscattering  measurements  as  a  function  of  angle  of  incidence.  Because  of  the  com¬ 
plicated  nature  of  the  seafloor  backscattering  process,  these  measurements  do  not  give 
a  unique  answer  about  the  nature  of  the  seafloor,  but  they  constitute  a  necessary  step 
towards  the  understanding  of  the  processes  at  work,  and  towards  the  construction  of 
empirical  formulas  for  sonar  design  purposes. 

Further  insight  into  the  processes  influencing  the  seafloor  acoustic  back- 
scatter  are  obtainable  with  multifrequency  arrays  as  they  allow  simultaneous  sam¬ 
pling  of  several  bottom  roughness  scales.  However  acoustic  measurements  alone  may 
never  be  sufficient  to  give  the  exact  nature  of  the  seafloor,  and  ground  truth  from 
some  other  sensing  technique  (e.g.  bottom  sampling  or  photography)  may  still  be 
necessary. 
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Inference  of  manganese  nodule  coverage  from  Sea  Beam  acoustic 
backscattering  data 


Christian  de  Moustier* 


ABSTRACT 

Normal  incidence  reflectivity  from  a  manganese 
nodule  field  was  measured  with  a  12  kHz  Sea  Beam 
multibeam  echo-sounding  system,  aboard  the  R/V 
Thomas  Washington  and  used  to  infer  nodule  coverage. 
A  reflectivity  map  of  the  area  was  produced  using  the 
intensity  of  the  specular  return  from  each  ping.  The 
patchiness  of  the  nodule  coverage  is  evidenced  by  defi¬ 
nite  highs  and  lows  in  the  reflectivity  pattern.  Ground 
truth  was  provided  by  near-bottom  acoustic  measure¬ 
ments  and  photographs  taken  with  the  Deep  Tow  in¬ 
strument  package  of  the  Marine  Physical  Laboratory. 
Agreement  between  the  simple  nodule  coverage  predic¬ 
tions  from  Sea  Beam  acoustic  data  and  the  bottom  pho¬ 
tographs  taken  throughout  the  area  is  98  percent.  Al¬ 
though  the  Sea  Beam  system  is  limited  in  its  dynamic 
range,  this  paper  shows  that  it  can  be  used  very  ef¬ 
fectively  to  determine  both  topography  and  nodule 
coverage  in  potential  mining  areas. 


INTRODUCTION 

Economic  considerations  make  the  exploitation  of  ferro¬ 
manganese  nodules  in  the  deep  ocean  strongly  dependent  upon 
extensive  assessment  of  the  density  of  nodule  coverage  at  po¬ 
tential  mining  sites.  Depending  upon  the  method  used  to 
survey  a  site,  the  cost  in  ship  time  and  manpower  can  become 
prohibitive.  The  approach  described  here  derives  its  cost- 
effectiveness  from  the  use  of  a  Sea  Beam  echo-sounding  system 
operating  at  optimum  ship  speed  (10-12  knots).  By  taking 
advantage  of  the  narrow-beam  and  multibeam  capabilities  of 
the  Sea  Beam,  I  was  able  to  combine  detailed  topographic 
mapping  with  bottom  backscattering  characteristics  to  infer 
nodule  coverage.  The  data  presented  here  were  gathered  during 
two  cruises  in  the  northeastern  tropical  Pacific  in  May  and 
June  of  1983.  The  first  cruise  (PASCUA  leg  5)  was  a  rapid  Sea 
Beam  survey  of  an  approximately  20  x  15  mile  area  southeast 
of  Deep  Ocean  Mining  Environmental  Study  (DOMES)  site  C 
(Figure  1 ),  on  an  Ocean  Mining  Associates  (OMA)  trial  mining 
site.  This  survey  provided  the  topographic  and  acoustic  reflec¬ 


tivity  maps  which  served  as  the  base  of  a  month-long,  fine-scale 
study  of  the  same  site,  with  the  Deep  Tow  instrument  package 
of  the  Marine  Physical  Laboratory  of  the  Scripps  Institution  of 
Oceanography  (Spiess  and  Tyce,  1973;  Spiess  and  Lonsdale, 

1982) .  The  purpose  of  the  Deep  Tow  survey  (ECHO  leg  1)  was 
to  assess,  using  a  new  multifrequency  array,  the  near-bottom 
backscattering  properties  of  a  manganese  nodule  field.  In  addi¬ 
tion.  an  environmental  impact  study  of  nodule  mining  was 
carried  out  through  an  extensive  box  coring  program  (Spiess  et 
a!.,  1984).  The  choice  of  the  site  was  motivated  by  the  existence 
of  a  good  DOMES  and  OMA  data  base  in  the  general  area. 

This  paper  discusses  the  results  of  the  Sea  Beam  acoustic 
study,  and  uses  the  Deep  Tow  data  as  ground  truth. 

BACKGROUND 

Manganese  nodules  remote  sensing 

Most  techniques  used  to  prospect  for  manganese  nodules  on 
the  deep  ocean  floor  rely  on  acoustic  remote  sensing,  near¬ 
bottom  photography  and/or  television,  bottom  sampling,  or  a 
combination  of  these  methods. 

Two  general  approaches  have  been  taken  for  acoustic  remote 
sensing  of  nodules.  The  first  is  a  single-frequency  approach, 
using  high-frequency  (  >  100  kHz)  side-scan  sonars  on  a  deeply 
towed  instrument  (Spiess.  1980).  or  using  a  shipboard  subbot¬ 
tom  profiler  to  correlate  the  thickness  of  the  transparent  layer 
with  nodule  abundance  (Piper  et  al.,  1979;  Mizuno  and  Mori- 
tani.  1976).  The  second  is  a  multifrequency  approach  [Sumi¬ 
tomo  Metal  Mining  Co.  (MFES),  Magnuson  et  al.,  1981,  1982; 
Spiess  et  al.,  1984]  which  uses  the  frequency  dependence  of 
sea-floor  acoustic  reflectivity  to  infer  nodule  sizes  and  abun¬ 
dance.  Both  approaches  rely  on  the  observation  of  a  difference 
in  acoustic  backscatter  between  a  sedimentary  bottom  laden 
with  nodules  and  a  bare  one. 

Near-bottom  photography  techniques  range  from  the  con¬ 
ventional  deeply  towed  camera  systems  to  the  unmanned,  deep 
diving,  free  vehicle  Epaulard  (Duranton  et  al.,  1980;  Galerne, 

1983) .  Ground  truth  for  remote  sensing  of  manganese  nodules 
is  obtained  through  direct  sampling  of  a  nodule  field.  To  this 
end,  dredging  is  used  where  large  amounts  of  nodules  with  little 
or  no  sediments  are  sought  On  the  other  hand,  box  cores 
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where  a  is  the  exponential  attenuation  coefficient  and  R  is  the 
Rayleigh  reflection  coefficient.  Equation  (2)  assumes  that  the 
bottom  reflects  the  sound  rather  than  scatters  tt.  This  assump¬ 
tion  seems  reasonable  for  a  manganese  nodule  field  whose 
aggregate  response  was  shown  theoretically  in  Magnuson 
( 1983)  to  be  reflective  rather  than  scattered. 

For  a  plane  wave  normally  incident  on  the  boundary  be¬ 
tween  two  media  of  respective  impedance  ptct  and  p2c2,  the 
Rayleigh  reflection  coefficient  is  the  ratio  of  the  reflected  to  the 
incident  pressure  waves: 


P: ci  ~  P 

p2c2  +  p,c,' 


(3) 


For  the  purpose,  of  this  paper,  medium  l  will  be  water  and 
medium  2  will  be  sediment  or  nodules.  The  fraction  on  the  right 
side  of  equation  (2)  indicates  that  sea  water  absorbs  some  of  the 
energy,  causing  an  exponential  loss  with  distance  from  the 
source,  and  that  intensity  decreases  proportionally  to  the 
square  of  distance  due  to  geometric  spreading.  This  is  summa¬ 
rized  in  the  sonar  equation  for  an  echo  sounder  obtained  by 
taking  10  log(0  of  both  sides  of  equation  (2): 


provide  an  undisturbed  sample  of  the  bottom  whenever  an 
assessment  of  the  geologic,  biological,  or  acoustical  properties 
of  both  the  nodules  and  their  underlying  sediments  is  needed. 

Normal  incidence  acoustic  reflectivity 

There  is  extensive  literature  dealing  with  acoustic  properties 
of  the  ocean  floor  at  various  frequencies  and  angles  of  inci¬ 
dence,  and  with  their  geologic  implications.  Tyce  (1976)  and 
Parrot  et  al.  (1980)  gave  good  reviews  of  the  literature  available. 

Due  to  its  multibeam,  narrow-beam  characteristics,  the  Sea 
Beam  system  allows  the  measurement  of  the  acoustic  back- 
scatter  of  the  sea  floor  from  each  individual  beam.  An  angular 
relationship  of  acoustic  backscattering  can  then  be  derived. 
However,  in  this  paper  I  take  a  first  look  at  the  data  by 
selecting  the  beam  nearest  normal  incidence  and  leave  the 
others  for  later  analysis. 

Since  the  depth  of  this  manganese  nodule  survey  area  is  large 
(4  500  m  on  the  average),  the  radius  of  curvature  of  the  acoustic 
wavefronts  reaching  the  bottom  is  much  larger  than  Sea  Beam’s 
acoustic  wavelength  (12  cm),  and  the  wavefronts  can  be  con¬ 
sidered  nearly  planar.  Therefore,  I  use  this  approximation  in 
the  following  development. 

In  general,  the  intensity  of  a  plane  acoustic  wave  is  related  to 
the  acoustic  pressure  P  by 


where  p  is  the  density  of  the  medium  and  c  the  propagation 
velocity  in  this  medium.  The  product  pc  is  called  the  character¬ 
istic  acoustic  impedance  of  the  medium.  The  bar  indicates  that 
the  pressure  is  averaged  over  some  time  (usually  the  integration 
time  of  the  instrument)  (Urick,  1983).  If  /0  denotes  the  intensity 
of  sound  at  unit  distance  from  the  source,  the  intensity  of  the 
normal  incidence  return  /  from  the  bottom  at  depth  r  is  given 
by 


Fig.  1.  Location  of  survey  area.  Deep  Ocean  Mining  Environ¬ 
mental  Studies  (DOMES)  sites  A  B,  C  are  shown  by  dots. 
Inset:  track  chart  of  the  PASCUA  expedition  leg  5. 


EL  =  SL  -  BL  —  TL,  (4) 

where 

EL  —  10  log  i  o  /  =  echo  level, 

SL  =  10  logl0/0  =  source  level, 

BL  =  -20  log  10R  -  bottom  loss. 


TL  =  20  log, 0(2r)  +  20ar  logl0e  =  transmission  loss. 

Equation  (4)  was  used  to  determine  the  reflection  coefficients  of 
various  marine  sediment  types  as  a  method  of  classification 
(Breslau,  1967).  Likewise,  Hamilton  (1970a,  b)  did  extensive 
work  toward  predicting  in-situ  acoustic  and  clastic  properties 
of  marine  sediments.  Tyce  (1976)  focused  on  the  determination 
of  the  attenuation  coefficient  in  ocean  sediments  with  a  4  kHz 
near-bottom  profiler.  His  data  indicate  a  large  variability  in 
seafloor  reflectivity  over  short  horizontal  distances  due  to 
changes  in  local  bottom  roughness  and  composition.  Breaker 
and  Wmokur  ( 1 967 >  observed  an  increase  in  the  magnitude  of 
the  echo  level  fluctuations  with  distance  from  the  bottom. 
Considering  that  most  of  the  acoustic  reflection  comes  from  the 
first  Fresnel  zone  which  increases  in  size  with  distance  from  the 
bottom,  they  attributed  these  fluctuations  to  changes  in  the 
reflective  characteristics  of  the  bottom  over  very  short  dis¬ 
tances.  Again,  these  changes  are  most  likely  due  to  variations  in 
bottom  substrate  and  or  bottom  roughness.  The  effect  of 
roughness  is  to  change  the  phase  relationships  between  the 
different  reflectors  within  a  Fresnel  zone.  How  rough  a  surface 
appears  depends  on  the  ratio  of  the  root-mean-square  (rms) 
bottom  roughness  a  to  the  acoustic  wavelength  k  (Clay  and 
Medwin.  1977,  chap.  10).  Using  the  vertical  component  of 
wavenumber  k,  one  can  distinguish  three  cases : 

(1)  Ka  <?  L  the  surface  appears  to  be  smooth  and  the 
amplitude  of  the  backscattered  return  is  dctei  mined  by 
the  reflection  coefficient  (coherent  return); 


and 


72 


$«a  Beam  Acoustic  Bacfcacattarlng  Data 


(2)  ko  >  1,  the  surface  appears  to  be  rough  and  the 
returns  ire  mostly  incoherent , 

(3)  ko  ^  1,  intermediate  range  bridging  the  coherent 
and  incoherent  regimes. 

Clay  and  Leong  (1973)  pointed  out  that  when  ko  ^  1  and 
when  the  spatial  dimensions  of  bottom  roughness  are  much 
smaller  than  the  diameter  of  the  first  Fresnel  zone,  the  reflective 
components  add  as  the  sum  of  squares  and  the  returned  mean 
square  signal  is  proportional  to  the  number  of  reflectors  in  the 
ensonified  area.  This  means  that  bottom  loss  depends  also  on 
the  ensonified  area  and  therefore  on  the  beam  width  of  the 
measuring  system.  The  diameter  of  the  first  Fresnel  2one  is 
given  by  d  *  (2Xh)‘  2  where  \  is  the  acoustic  wavelength  and  h 
is  the  distance  from  the  sonar  to  the  reflection  plane  (Clay  and 
Medwin,  1977,  p.  50).  For  hull-mounted  sonars  such  as  Sea 
Beam,  this  distance  is  approximately  the  water  depth.  In  this 
survey  area  Sea  Beam's  12  cm  wavelength  yields  a  diameter  of 
about  33  m  for  the  first  Fresnel  zone  at  4  500  m  depth.  As 
shown  later,  the  rms  roughness  in  this  manganese  nodule  area 
is  about  2  cm,  so  that  ko  2s  1.  and  the  spatial  dimensions  of 
roughness  ( 10  cm  or  less)  are  much  smaller  than  the  dimensions 
of  the  first  Fresnel  zone. 

Other  investigators  have  tried  to  correlate  acoustic  frequency 
dependence  of  bottom  reflectivity  with  bottom  type  and  or 
bottom  roughness  (Mackenzie,  1960;  Zhitkovskii  et  al.,  1966). 
This  turns  out  to  be  a  difficult  problem  when  the  acoustic 
wavelength  and/or  the  bottom  are  such  that  subbottom  reflec¬ 
tions  can  no  longer  be  ignored.  Note  that  when  the  bottom  is 
assumed  horizontally  stratified  with  plane  interfaces  between 
layers,  subbottom  reflections  from  thin  layers  (Clay  and 
Medwin.  1977,  chap.  2)  yield  a  closed-form  solution.  However, 
in  order  to  perform  a  first-order  analysis,  simplifying  assump¬ 
tions  considering  the  ocean  bottom  as  locally  homogeneous  in 
depth  are  made.  An  ideal  ocean  bottom  can  then  be  defined 
where,  for  normal  incidence,  the  dimensions  of  the  bottom 
irregularities  compared  to  the  acoustic  wavelength  and  the 
impedance  mismatch  at  the  water-sediment  interface  are  the 
main  parameters  in  the  determination  of  bottom  reflectivity. 
Higher  reflectivity  corresponds  to  a  high  impedance  contrast 
between  sea  water  and  bottom  (generally  hard  substrate)  and  a 
locally  smooth  relief  (dimensions  smaller  than  the  acoustic 
wavelength).  As  the  ratio  of  the  size  of  the  bottom  irregularities 
to  the  acoustic  wavelength  increases,  the  bottom  irregularities 
act  as  scatterers  and  the  specular  component  of  an  echo  is 
lower  Reflectivity  is  again  lower  when  the  impedance  contrast 
is  low  (soft  bottom  with  impedance  similar  to  that  of  sea  water). 

A  manganese  nodule  field  consists  of  a  soft  bottom  (sedi¬ 
ments)  on  which  a  variable  number  of  nodules  lie.  Nodule 
acoustic  remote  sensing  is  then  possible,  provided  there  is 
enough  impedance  contrast  between  nodules  and  sediments. 

Although  no  values  are  presently  available  for  the  acoustic 
impedance  of  manganese  nodules  or  surface  sediments  in  the 
survey  area,  the  values  given  in  Hamilton  ( 1970b)  for  sediments 
of  the  Pacific  and  in  Magnuson  et  al.  (1981)  for  nodules  can  be 
used  to  compute  an  expected  acoustic  reflectivity  difference 
between  nodules  and  their  underlying  sediments.  The  sediments 
are  mostly  pelagic  silty  clay  (Spiess  et  al.,  1984)  for  which 
Hamilton  gave  an  in-situ  reflective  coefficient  of  0.1316  at 
normal  incidence  and  the  corresponding  bottom  loss  of  17.6 
dB.  For  Pacific  nodules,  Magnuson  gave  a  range  of  values  for 
both  nodule  density  and  compressionai  wave  speed.  I  use  p2  * 


1.95  g/cm3,  c2  *  2  400  m/s  for  a  nodule,  and  from  Urick  (1983) 
p,  —  1.0475  g/cm\  ct  —  1  520  m/s  fo*-  seawater  at  4  500  m 
depth.  Substituting  these  values  in  equation  (3)  yields  a  reflec¬ 
tion  coefficient  R  ■*  0.4923,  which  when  expressed  as  the 
bottom  loss  term  in  equation  (4)  becomes  BL  —  6.2  dB.  When 
the  bottom  depth  is  relatively  constant,  the  transmission  loss  of 
equation  (4)  can  be  assumed  constant.  Therefore  bottom  loss  is 
the  only  variable  in  equation  (4),  and  the  difference  in  reflec¬ 
tivity  between  a  nodule  and  the  underlying  sediment  is  simply 
the  difference  in  bottom  loss:  1 1.4  dB. 

Hamilton  (1970b)  also  indicated  there  is  little  or  no  depen¬ 
dence  of  bottom  loss  with  acoustic  frequency  for  sediments 
where  there  is  no  lower  layer  to  reflect  sound  which  interferes 
with  sea-floor  reflections.  Such  is  not  the  case  for  a  nodule  field. 
Studies  by  Magnuson  et  al.  (1981,  1982)  and  the  Sumitomo 
Metal  Mining  Co.,  where  manganese  nodules  were  modeled  as 
spheres  on  an  infinite  flat  plane,  show  that  the  bottom  reflec¬ 
tivity  of  a  nodule  field  can  be  expected  to  reach  a  maximum  for 
a  frequency  corresponding  to 

kn  -  1,  (5) 

where  k  is  the  acoustic  wavenumber  and  a  is  the  mean  nodule 
radius.  By  comparing  the  results  from  measurements  at  sea 
with  the  foregoing  models,  some  clues  arise  as  to  the  validity  of 
the  assumptions  made.  However,  before  I  can  discuss  these 
results,  I  need  to  characterize  the  data  base  and  the  measuring 
systems.  Sea  Beam  and  Deep  Tow,  used  in  this  survey 

SEA  BEAM  ACOUSTICS 

The  Sea  Beam  system  uses  a  multibeam,  narrow-beam  echo 
sounder  operating  at  a  frequency  of  12.158  kHz  with  a  pulse 
length  of  7  ms  and  an  echo  processor  to  generate,  in  near  real 
time,  contour  maps  of  the  ocean  floor  while  the  ship  is  under- 


Fig.  2.  Sea  Beam  geometry  The  transmitted  beam  pattern 
ensonifies  an  area  on  the  sea  floor  2§  degrees  in  the  fore-ait 
direction  and  54  degrees  athwartships,  in  a  vertical  plane  from 
the  ship.  The  receiving  beam  pattern  appears  on  the  sea  floor  as 
16  rectangles,  each  20  degrees  fore-aft  by  2$  degrees  athwart¬ 
ships.  The  received  acoustic  energy  therefore  comes  from  the 
intersection  of  the  two  beam  patterns  which  is  delimited  by  16 
squares  2  j  degrees  on  a  side 
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way.  Figure  2  illustrates  the  Sea  Beam  transmit/receive  geome¬ 
try.  The  transmitted  beam  pattern  spans  54  degrees  athwart- 
ships  by  2  $  degrees  in  the  fore-aft  direction.  It  is  pitch-stabilized 
to  ensure  vertical  projection.  As  a  result  of  the  Sea  Beam  beam 
forming,  the  receiving  beam  pattern  can  be  approximated  by  16 
adjacent  rectangles  2  2/3  degrees  by  20  degrees.  The  acoustic 
energy  received  at  the  ship  comes  from  the  intersection  of  the 
transmit  and  receive  beam  patterns,  that  is  16  “squares”  2  2/3 
on  a  side.1 

With  contour  mapping  as  a  goal,  the  Sea  Beam  does  not 
preserve  all  the  acoustic  information  it  receives.  As  a  result 
some  valuable  information  not  necessary  for  depth  determi¬ 
nation,  such  as  the  amplitude  and  the  shape  of  the  echo  signal, 
is  lost. 

In  order  to  preserve  the  echo  amplitude  information,  a  paral¬ 
lel  data  acquisition  system  was  built  around  an  LSI  11/23 
minicomputer  (de  Mousticr,  1985).  This  allows  digital  record¬ 
ing  on  magnetic  tape  the  16  Sea  Beam  detected,  nonroll- 
compensated  beams,  along  with  time-varied-gain  (TVG)  and 
ship  roll  information.  The  TVG  is  designed  to  compensate  for 
transmission  loss  [attenuation  and  spreading  terms  in  equation 
(2)]  and  has  already  been  applied  by  the  Sea  Beam  system  when 
the  data  are  recorded.  Figure  3  shows  a  typical  16  beam  return, 
not  roll-compensated,  where  the  envelopes  of  the  detected 
acoustic  signals  are  plotted  in  analog-to-<ligitai  (A/D)  units 
versus  time  after  transmission.  These  units  correspond  directly 
to  sound  pressure  levels  with  TVG  applied  In  this  instance,  the 
sea  floor  is  essentially  flat,  resulting  in  the  parabola  outlined  by 
the  16  returns.  The  amplitude  is  highest  in  the  near-specular 
direction  {first  arrival  in  time)  and  decreases  rapidly  from  the 

‘For  a  detailed  description  of  the  system,  sec  Rcnard  and  Allenou 
(1979). 


Center  outward  The  early,  synchronous  returns  on  the  side 
beams  represent  side  lobe  energy  from  the  specular  return  and 
are  rejected  in  both  depth  and  return  amplitude  processing. 

The  Sea  Beam  system  has  not  been  calibrated;  therefore  in 
the  following  all  the  measurements  are  presented  as  relative. 

DATA 

The  data  presented  here  were  collected  during  a  40-hour  Sea 
Beam  survey  and  a  subsequent  Deep  Tow  near  bottom  investi¬ 
gation.  Although  the  Sea  Beam  system  maps  a  swath  of  sea 
floor  three-quarters  of  the  ocean  depth  wide,  a  high  degree  of 
overlap  between  swaths  (80  percent  in  places)  was  used  to 
provtde  as  much  redundant  acoustic  data  as  possible.  Figure  4 
shows  the  topography  of  the  survey  area  contoured  from  the 
Sea  Beam  data.  The  contours  are  at  20  m  intervals,  implying  a 
relatively  flat  portion  of  sea  floor  with  less  than  200  m  of  total 
relief  in  40  km.  Aside  from  two  north-south  troughs  which 
bottom  out  at  4  600  m  on  the  edges  of  the  map.  the  mean  depth 
is  4  500  m. 

The  general  pattern  agrees  well  with  the  known  north-south 
orientation  of  sea-floor  features,  perpendicular  to  the  spreading 
direction  in  this  part  of  the  Pacific  Ocean. 

At  this  depth  the  Sea  Beam  system  transmits  every  8  s.  F'or 
each  transmission,  a  set  of  16  returns  (Figure  3)  is  digitized  and 
recorded  on  magnetic  tape.  For  a  first  look  at  these  data,  the 
postprocessing  consisted  in  selecting  the  highest  amplitude  in 
each  set.  The  highest  amplitude  generally  represents  the  near- 
specular  return,  which  in  the  case  of  a  relatively  flat  bottom,  as 
in  this  area,  comes  in  first. 

This  set  of  peak  amplitudes  ts  then  low-pass  filtered  by 
applying  a  running  mean  along  the  track,  averaging  over  the 


Fig.  3.  Acoustic  signal  envelopes  of  the  16  performed  beams  at  the  output  of  the  Sea  Beam  echo-processor  receivers.  The  jc-axis 
represents  time  after  transmission  in  seconds.  The  v-axis  represents  the  amplitude  of  the  signals  expressed  in  A  D  units  (0-4095) 
linear  scale  Such  data  are  recorded  digitally  on  magnetic  tape  every  transmission  cycle,  along  with  time-vaned-gam  and  ship’s  roll 
No  roll  correction  has  been  applied  to  the  data  at  this  stage. 
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Fig.  4  Contoured  map  of  the  survey  area  from  Sea  Beam  data.  The  contours  are  in  un corrected  meters.  The  contour  interval 

is  20  m. 
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Fig.  5.  Ship’s  track  of  the  Sea  Beam  survey  in  the  manganese  nodule  area  (PASCUA  leg  5).  The  peak  amplitude  from  each  ping 
(Figure  31  has  been  selected,  low-pass  filtered,  and  plotted  along  the  track  in  A  D  units  linear  scale  10-4095).  Five  intervals  were 
chosen  on  this  scale  and  pattern  coded  as  shown  in  the  key. 
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size  of  the  vertical  beam’s  footprint .  roughly  five  transmission 
cycles  at  1 1  knots.  The  resulting  data  are  plotted  and  pattern 
coded  along  the  ship’s  track  on  Figure  5.  The  A/D  units  linear 
scale  is  the  same  as  indicated  in  Figure  3. 

The  Deep  Tow  survey  was  concentrated  in  the  southern 
portion  of  the  area  where  the  Sea  Beam  survey  density  is 
greatest.  The  Sea  Beam  acoustic  data  served  as  a  reference  for 
the  Deep  Tow  work  during  which  camera  runs,  near-bottom 
multifrequency  (4.5,  9,  15,  30,  110,  160  kHz)  backscattering 
measurements,  and  side-looking  sonar  data  (110  kHz)  were 
collected.  As  part  of  a  mining  environmental  impact  study,  16 
box  cores  were  also  taken  throughout  the  Deep  Tow  survey 


area.  The  various  elements  of  the  Deep  Tow  work  are  summa¬ 
rized  in  Figure  6  (Spiess  et  al.,  1984). 

DISCUSSION 
Sea  Beam  reflectivity  data 

Inspection  of  the  peak-amplitude  data  reveals  a  substantial 
ping-to-ping  variability  as  well  as  definite  mean  amplitude 
highs  and  lows  along  the  ship’s  track.  This  point  is  illustrated  in 
Figure  7  where  amplitude  versus  distance  along  the  track  is 
plotted  in  conjunction  with  the  center  beam  width  profile.  The 
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Fig.  6.  Deep  Tow  track  chart  showing  camera  runs,  near-bottom  multifrequency  backscattering  measurement  runs,  box  core 
locations,  and  navigation  transponder  positions.  Coordinates  are  in  meters  east  and  north  of  14°34'N,  I25°3(yw. 
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amplitude  is  expressed  in  relative  decibels,  with  reference  to  the 
mean  amplitude  of  the  data  displayed.  The  data  shown  in  this 
figure  have  been  averaged  over  five  transmission  cycles  for 
clarity;  actual  ping-to-ping  variability  is  much  greater.  Some 
small  correlation  with  topography  can  be  seen.  However,  re¬ 
member  that  the  vertical  exaggeration  on  the  depth  profile  is 
18.  The  slopes  are  therefore  rarely  in  excess  of  2  degrees.  On  the 
left  side  of  Figure  7,  large  fluctuations  in  amplitude  are  appar¬ 
ent,  and  the  average  value  is  around  —6  dB.  Compared  with 
the  right  side  where  the  average  amplitude  is  -1-4  dB,  the 
bottom  reflectivity  has  changed  significantly  from  one  side  to 
the  other. 

Following  this  result,  I  divided  the  peak  amplitude  data  set 
into  five  intervals  on  the  A/D  linear  scale  (0—4  095)  and  pattern 
coded  each  one  as  shown  in  Figure  5.  The  data  match,  to  within 
130  A/D  units  on  the  average,  over  all  except  three  track 
intersections  where  there  are  slight  offsets.  Because  normal 
satellite  navigation  is  simply  inadequate  for  track  positioning 
at  the  level  of  Sea  Beam’s  resolution,  one  relies  upon  matching 
Sea  Beam  topography  at  crossings  and  overlaps  for  accurate 
track  positioning.  However,  accurate  adjustments  are  difficult 
to  achieve  when  the  sea  floor  surveyed  exhibits  little  relief,  as  in 
this  area. 

A  histogram  of  the  data  distribution  is  shown  in  Figure  8. 
The  histogram  suggests  that  the  data  suffer  amplitude  clipping 
at  3  000  A/D  units,  since  there  are  no  values  beyond  this  point. 
This  unfortunate  characteristic  is  an  artifact  of  the  data  acqui¬ 
sition  system  which  was  saturating  at  7.3  V.  This  problem  has 
since  been  corrected.  In  addition,  it  was  found  that  the  output 
of  the  Sea  Beam  echo  processor  receivers  is  limited  to  8.5  V 
instead  of  the  10  V  announced  by  the  manufacturer.  As  a  result, 
the  receivers’  dynamic  range  is  apparently  inadequate  to  ac¬ 
commodate  both  the  low  backscatter  levels  received  on  the  side 
beams  from  a  smooth  bottom  and  the  high  amplitude  associ¬ 


ated  with  some  strong  specular  returns.  For  this  simplified 
application,  however,  the  nodule  distribution  results  are  not 
significantly  affected,  since  I  arbitrarily  chose  to  start  the  last 
A,  D  interval  at  2  500  units.  Nevertheless,  this  amplitude  clip¬ 
ping  does  limit  resolution  for  high-amplitude  returns  and  pre¬ 
vents  discussion  of  the  abilities  of  such  a  sonar  to  characterize 
similar  nodule  fields  or  other  seafloor  types.  From  discussions 
with  the  Sea  Beam  manufacturer,  it  seems  that  a  first-order 
improvement  to  this  dynamic  range  limitation  can  be  achieved 
by  a  simple  modification  of  the  detection  amplifiers.  Obtaining 
the  full  range  needed,  however,  might  require  modification  of 
the  hydrophone  preamplifiers  which  have  been  found  to  satu¬ 
rate  occasionally  and  would  also  require  use  of  logarithmic 
detection  amplifiers. 

Drawn  from  the  data  in  Figure  5,  a  contour  map  of  reflec¬ 
tivity  in  the  southern  portion  of  the  survey  area  is  presented  in 
Figure  9.  It  uses  the  same  pattern  code  as  previously  described. 
Two  main  characteristics  stand  out  on  this  map:  a  region  of 
strong  reflectivity  in  the  center,  bordered  to  the  east  and  to  the 
west  by  strips  of  much  lower  reflectivity.  Some  of  the  low  values 
correspond  to  the  north-south  troughs  apparent  in  the  bathy¬ 
metry  (Figure  4).  However,  the  low  reflectivity  extends  3  to  4 
km  into  the  flat  center  portion.  Changes  in  relative  reflectivity 
of  12  dB  or  more  (Figure  7)  can  therefore  be  expected  as  one 
crosses  the  leveled  ground  area  in  an  east-west  direction.  Al¬ 
though  this  number  is  biased  in  trying  to  account  for  the 
clipping  of  the  high-amplitude  return  (I  added  2  dB  to  the  4  dB 
average  amplitude  shown  on  the  right  side  of  Figure  7),  it  is 
reassuring  to  note  that  it  is  close  to  the  expected  value  com¬ 
puted  in  the  “Background”  section  (11.4  dB).  Since  we  were 
surveying  a  manganese  nodule  site,  l  would  like  to  associate  the 
above  changes  with  differences  in  nodule  coverage,  with  high- 
amplitude  values  corresponding  to  a  dense  coverage  and  low 
values  indicating  absence  of  nodules.  To  do  so,  however,  re- 


Fig.  7  The  peak  amplitude  from  each  transmission  cycle  (Figure  3)  is  selected,  low-pass  JU*  T^^TomsocmSnR^enlwr'bwS^Septh 
in  decibels  normalized  by  the  mean  of  the  portion  of  data  shown.  The  mean  is  at  0  dB.  The nLth 7s  oven  in 
profile  is  plotted  with  a  vertical  exaggeration  of  18  to  give  an  indication  of  slope  along  the  sh.ps  trac  p  & 
uncorrected  meters  and  distance  in  nautical  miles. 
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quires  ground  truth  available  only  from  near-bottom  instru¬ 
ments  with  bottom  photographic  capabilities.  Therefore,  this 
Sea  Beam  reflectivity  work  was  done  in  an  area  intended  for 
Deep  Tow  coverage. 

Deep  Tow  ground  truth 

It  was  a  rewarding  experience  to  watch  Deep  Tow  camera 
runs  unveil  snapshots  of  the  sea  floor  which  matched  the  simple 
predicted  bottom  pattern.  Each  camera  run  (Figure  6)  consis¬ 
tently  output  pictures  which  agreed  with  the  simple  manganese 
nodule  coverage  inferences  made  from  the  Sea  Beam  reflectivity 
data.  Three  representative  pictures  are  shown  on  Figure  10. 
Each  picture  typifies  an  amplitude  interval.  The  dense  nodule 
coverage  of  photo  A  corresponds  to  the  dot  patterned  areas  on 
Figure  9  (2  500  A/D  units  and  greater).  The  intermediate 
nodule  coverage  seen  in  photo  B  is  found  in  areas  with  the 
brick  pattern  (2  000-2  500  A/D  units)  and  the  hatched  pattern 
(2  000-1  500  A/D  units).  The  two  remaining  patterns  (1  500 
A/D  units  or  lower)  turned  out  to  be  bare  mud  as  shown  on 
photo  C.  Inspection  of  over  3  000  Deep  Tow  bottom  photo¬ 
graphs  and  classification  according  to  the  three  types  shown 
(dense,  intermediate,  and  bare)  yielded  a  98  percent  correlation 
with  the  simple  nodule  coverage  made  from  the  Sea  Beam  data. 

During  a  camera  run,  the  towed  instrument  is  flown  approxi¬ 
mately  10  m  above  the  bottom.  This  renders  pictures  which 
cover  about  35  m2  of  sea  floor,  thus  giving  a  clue  to  nodule 
sizes.  From  the  photographs  and  from  the  I  4  m2  box  cores 
(photo  D),  the  average  nodule  was  found  to  be  4  to  8  cm  in 
diameter,  with  a  flattened  and  irregular  caulifiower-likc  shape, 
as  opposed  to  the  rounder  nodules  found  in  other  areas 
(Heezen  and  Hollister,  1971 ,  Bischoff  and  Piper,  1979). 

The  first  results  from  the  Deep  Tow  multifrequency,  near¬ 
bottom  acoustic  bickscattering  measurements  (Spiess  et  a]., 
1984)  also  reinforce  the  belief  that  Sea  Beam  is  well  suited  to 


AMPLITUDE  x  100 

Fig.  8.  Histogram  of  peak  amplitude  data  distribution  from  the 
data  shown  in  Figure  5.  20  percent  of  the  values  fall  in  the 
2  900-3  000  A/D  units  range  due  to  saturation  in  the  data 
acquisition  system. 


map  out  manganese  nodule  coverage  in  this  prospective  mining 
area.  Figure  1 1  shows  the  normal-incidence  reflectivity  differ¬ 
ence  between  a  densely  nodule-covered  patch  and  a  barren  one, 
as  measured  at  the  six  frequencies  of  the  Deep  Tow  backscatter 
system.  The  optimum  frequency  for  this  site  appears  to  be 
between  15  kHz  (16  dB)  and  9  kHz  (13  dB).  The  operating 
frequency  of  the  Sea  Beam  system  (12  kHz)  falls  exactly  halfway 
between  these  two  frequencies.  At  12  kHz  equation  (5)  yields  a 
mean  noduie  radius  of  2  cm  (1.6  cm  at  15  kHz).  Although  the 
nodules  of  this  site  are  poorly  approximated  by  spheres,  the 
foregoing  results  indicate  that  the  Sea  Beam  frequency  is  near 
optimum  for  the  area  surveyed.  However,  do  not  presume  that 
12  kHz  will  be  optimum  for  all  types  of  nodules,  or  that  any  12 
kHz  sonar  would  provide  the  same  results,  since  Sea  Beam  is  a 
roll-  and  pitch-stabilized  system  with  2$  degree  resolution.  As 
pointed  out  earlier,  this  high  resolution  has  a  significant  effect 
on  observed  backscattered  levels  and  associated  fluctuations 
due  to  the  size  of  the  observed  area.  In  comparison,  normal 
echo  sounders  typically  have  a  beam  width  between  30  and  60 
degrees  and  a  bandwidth  wider  than  Sea  Beam’s,  resulting  in  a 
substantially  lower  signal-to-noise  ratio  Both  the  Deep  Tow 
and  the  Sea  Beam  data  point  to  the  patchiness  of  the  nodule 
coverage,  some  patches  being  less  than  a  kilometer  in  extent.  It 
is  doubtful  whether  an  unstabilized  system  which  spans  an  area 
2.4  km  wide  (30-degree  beam  at  4  500  m)  can  resolve  such 
patches. 

Note  that  the  value  obtained  for  mud-to-nodule  reflectivity 
difference  with  Sea  Beam  is  lower  than  those  obtained  with 
Deep  Tow  at  9  and  15  kHz.  Even  though  the  frequency  depen¬ 
dence  of  mud-to-nodule  reflectivity  might  not  be  linear,  it  is 
reasonable  to  expect  the  value  at  12  kHz  to  be  between  those  at 
9  and  15  kHz.  To  explain  this  discrepancy  the  differences  in  the 
geometry  of  measurement  between  the  two  systems  should  be 
considered.  Indeed,  for  Sea  Beam  the  area  of  the  first  Fresnel 
zone  is  approximately  33  m  wide  at  4  500  m  depth,  which  is 
almost  one  order  of  magnitude  larger  than  the  diameter  of  the 
first  Fresnel  zone  (4  m  at  15  kHz)  observed  with  Deep  Tow 
w  hen  the  instrument  is  80  m  off  the  bottom.  However,  the  Deep 
Tow  acoustic  data  have  been  averaged  over  spatial  distances 
along  track  commensurate  with  Sea  Beam’s  normal  incidence 
footprint.  Moreover,  the  beam  width  dependence  presumably 
cancels  out  when  the  reflectivity  difference  between  nodule  and 
mud  is  measured  with  the  same  instrument.  In  the  reflective 
model  the  geometric  argument  is  then  ruled  out.  To  account  for 
the  Sea  Beam  s  mud-to-nodule  reflectivity  difference  being 
lower  than  those  of  the  two  neighboring  Deep  Tow  frequencies 
(9,  15  kHz),  I  therefore  need  to  reconsider  the  reflectivity  model. 
In  this  manganese  nodule  area,  the  local  roughness  over  a  200 
m  honzontal  extent  is  primarily  due  to  the  presence  or  absence 
of  nodules.  From  the  bottom  photographs  and  core  samples,  I 
estimate  the  rms  roughness  to  be  about  2  cm  and  spatial 
dimensions  of  roughness  to  be  10  cm  or  less.  At  9  and  15  kHz, 
80  m  off  the  bottom,  the  diameter  of  the  first  Fresnel  zone  (5.2 
and  4  m,  respectively)  is  then  still  substantially  larger  than  the 
spatial  dimensions  of  the  bottom  roughness  and  kc  *  0.75  and 
1.25,  respectively.  As  mentioned  in  the  “Background”  section, 
ko  =  1  corresponds  to  the  ill-defined  boundary  between  reflec¬ 
tion  and  scattering. 

If  I  now  take  the  difference  between  the  mud-to-nodule  re¬ 
flectivity  difference  (16  dB)  measured  at  15  kHz  and  the  com¬ 
puted  bottom  loss  for  the  mud  (17.6  dB),  I  obtain  a  bottom  lots 
for  nodules  of  1.6  dB.  In  comparison,  a  smooth  layer  of  chert  (a 
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Fig.  9.  Contour  map  of  relative  reflectivity  drawn  from  the  lower  half  of  Figure  5,  where  ship  track  density  is  greatest.  Same  pattern 
code  as  in  Figure  5.  The  dotted  line  delimits  the  Deep  Tow  survey  area  shown  in  Figure  6. 
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sedimentary  rock  having  an  acoustic  impedance  almost  4  times 
greater  than  the  estimated  impedance  of  a  nodule)  has  a  bottom 
loss  of  3  dB  (Tyce  et  al.,  1980).  The  estimate  of  the  mud  bottom 
loss  is  perhaps  questionable.  However,  it  is  reasonable  to  be¬ 
lieve  that  this  value  is  within  2  dB  of  the  real  one.  Moreover,  the 
nodules  are  porous  and  do  not  form  a  continuous  slab,  but 
rather  are  scattered  at  variable  distances  from  each  other. 
Therefore  the  bottom  loss  of  such  an  assemblage  is  expected  to 
be  greater  than  the  3  dB  given  for  a  chert  layer.  Clearly,  a 
simple  reflective  model  is  insufficient  to  explain  the  high  value 
found  at  15  kHz. 

A  more  realistic  result  is  obtained  by  considering  the  nodules 
as  a  thin  layer  separating  water  from  mud.  In  this  case,  the 
reflection  coefficient  at  normal  incidence  (Clay  and  Medwin, 
1977,  chap.  2)  is  given  by 


where  K,  is  the  water-nodule  reflection  coefficient  [as  given  by 
equation  (3)],  R2  is  the  nodule-mud  reflection  cocfficient,/is  the 
acoustic  frequency,  h  is  the  thickness  of  the  thin  layer,  Cnod  is 
the  compressional  wave  speed  in  nodules.  Using  p  =  1.37 
g/cm3  and  C  =  1  507  m/s  for  mud  (Hamilton,  1970b),  h  =  4  cm 
for  the  thickness  of  the  thin  layer,  and  the  values  given  in  the 
“Background”  section  for  the  other  parameters,  yields  a  bottom 
loss  ofJ3  dB  at  15  kHz.  Also  with  these  values,  the  reflection 
coefficient  R  of  equation  (6)  reaches  a  maximum  at  15  kHz. 
Therefore,  this  thin-layer  model  agrees  with  the  frequency  de¬ 
pendence  shown  in  the  Deep  Tow  data  (Figure  11),  but  does 
not  match  the  mud-to-nodule  reflectivity  difference  measured 
at  15  kHz,  especially  since  (as  mentioned  above)  the  nodules  do 
not  form  a  continuous  slab.  A  scattering  model  should  then  be 
considered  for  nodules.  However,  it  requires  knowledge  of  the 
ensomfied  areas  and  therefore  of  the  beam  patterns  of  the  Deep 
Tow  backscattering  assembly.  These  beam  patterns  were  not 
available  at  the  time  of  this  writing,  precluding  this  analysis. 


Fig.  10.  Photographic  ground  truth  (a)  dense  nodule  coverage  representative  of  the  dot  patterned  areas  in  Figure  9,  (b)  intermedi¬ 
ate  coverage  representative  of  the  brick  and  hatched  patterns;  (c)  bare  mud  representative  of  peak  amplitude  values  1  500  A/D 
units  and  lower;  (d)  1/4  m2  box  core  number  358  (Figure  6)  giving  clues  to  representative  nodule  sizes  for  the  area. 
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Further  analysis  of  the  data  preeented  here,  on  a  scattering 
rather  than  a  reflective  basis,  and  investigation  of  the  angular 
dependence  of  acoustic  beckscattenng  from  the  Sea  Beam  and 
the  Deep  Tow  data  should  help  explain  the  processes  at  work 
here. 

Father  ase  of  Sea  Beam's  acoustic  returns 

Using  only  the  peak  amplitude  of  the  specular  return  of  the 
16  Sea  Beam  preformed  beams  to  map  out  normal  incidence 
bottom  reflectivity  is  a  nice  first-order  application.  It  yields 
good  results  over  this  particular  manganese  nodule  field.  Clear* 
jy,  more  can  be  expected  from  the  Sea  Beam  acoustic  data. 
Work  is  currently  in  progress  to  model  the  angular  dependence 
of  Sea  Beam  acoustic  backscatter  over  various  types  of  sea  floor 
and  to  relate  the  data  to  the  geologic  characteristics  of  the 
bottom  (Patterson,  1967).  Following  work  done  by  Clay  and 
Leong  (1973)  and  recently  by  Stanton  (1984)  the  variations  in 
the  shape  of  the  echo  envelope  can  also  be  used  to  predict 
bottom  roughness.  According  to  Stanton,  his  statistical  analysis 
works  best  for  amplitudes  of  sea  floor  roughness  less  than  a 
quarter  of  the  acoustic  wavelength. 

Since  the  Sea  Beam  acoustic  wavelength  is  about  12  cm  and 
since  I  estimated  the  local  rms  bottom  roughness  to  be  about  2 
cm,  the  data  from  this  manganese  nodule  area  qualify  for  such  a 
statistical  analysis  which  could  prove  useful  for  discrimination 
between  various  bottom  substrates  with  similar  reflectivity  pat¬ 
terns. 

CONCLUSIONS 

In  a  recent  monograph  on  polymetallic  nodules  (Rapport  de 
l’Academie  des  Sciences,  Paris,  1984)  the  Sea  Beam  system  was 
deemed  an  essential  tool  for  any  large-scale  nodule  prospection 
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Fig.  1 1.  Frequency  dependence  of  relative  reflectivity,  as  mea¬ 
sured  with  the  Deep  Tow  multifrequency  acoustic  array.  Each 
point  represents  the  difference  in  reflectivity  between  a  densely 
nodule  covered  patch  and  bare  mud  at  the  acoustic  frequency 
indicated.  For  the  data  shown  in  this  figure,  the  transmit  pulse 
lengths  at  the  various  frequencies  are:  0.9  ms  (4.5  and  9  kHz); 
1.0  ms  (30  and  100  kHz),  and  1.1  ms  (15  and  60  kHz).  The  Deep 
Tow  instruction  package  is  approximately  80  m  off  the  bottom. 


since  detailed  knowledge  of  the  bathymetry  is  a  prerequisite  to 
nodule  mining. 

The  data  presented  here  bring  a  new  dimension  to  the  Sea 
Beam  multibeam  echo-sounding  system.  In  addition  to  contour 
mapping,  I  have  shown  that  it  is  possible  to  detect  changes  in 
reflectivity  along  the  track  by  using  the  amplitude  of  the  near- 
specular  acoustic  return.  This  works  particularly  well  over  a 
manganese  nodule  field  with  little  relief,  since  nodules  and  mud 
provide  a  definite  acoustic  contrast  (1 1  dB  or  more).  Although  a 
scattering  model  would  seem  more  appropriate  than  the  simple 
reflective  one  used  here,  the  normal-incidence  backscattering 
data  gathered  with  the  Deep  Tow  multifrequency  acoustic 
array  confirm  that  high  reflectivity  contrast  between  nodules 
and  mud  should  be  expected  at  Sea  Beam  operating  frequency 
(12  kHz).  This  frequency  seems  to  be  near  optimum  for  this 
manganese  nodule  site.  Bottom  photographs  taken  throughout 
the  site  agree  surprisingly  well  with  the  simple  nodule  coverage 
prediction  made  from  Sea  Beam  acoustic  data. 

In  spite  of  dynamic  range  limitations  evidenced  in  my  data,  I 
believe  that  nodule  coverage  assessment  can  be  performed  very 
effectively  with  a  Sea  Beam  system  while  contour  mapping  the 
topography  of  a  prospective  mining  site.  In  essence,  areas 
where  amplitude  saturation  occurred  correspond  to  areas  of 
denser  nodule  coverage  with  strong  reflectivity  Photographic 
ground  truth  can  be  obtained  with  an  unmanned,  deep-diving 
free  vehicle,  with  towed  vehicle  camera  and  sonar  systems,  or 
even  with  conventional  bottom  cameras. 

More  extensive  processing  of  Sea  Beam’s  detected  echoes  is 
in  progress  with  the  intent  of  assessing  whether  the  back- 
scattered  returns  separated  by  Sea  Beam’s  high  angular  resolu¬ 
tion  can  be  correlated  with  seafloor  geologic  characteristics. 
This  includes  identification  of  the  reflectivity  and  roughness 
characteristics  of  different  bottom  types  sediments,  rocks,  rock 
outcrops,  manganese  nodules,  etc.  Initial  results  have  been  very 
encouraging,  as  suggested  both  in  this  papier  and  in  ongoing 
work  on  various  typjes  of  sea  floor. 
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Chapter  IV 


BEYOND  BATHYMETRY:  MAPPING  ACOUSTIC  BACICSCATTERING 
FROM  THE  DEEP  SEAFLOOR  WITH  SEA  BEAM. 

IV.  1.  ABSTRACT 

In  its  standard  mode  of  operation,  the  multibeam  echo-sounder  Sea  Beam 
produces  high  resolution  bathymetric  contour  charts  of  the  seafloor  surveyed.  How¬ 
ever,  additional  information  about  the  nature  of  the  seafloor  can  be  extracted  from 
the  structure  of  the  echo  signals  received  by  the  system.  Such  signals  have  been 
recorded  digitally  over  a  variety  of  seafloor  environments  for  which  independent  obser¬ 
vations  from  bottom  photographs  or  sidescan  sonars  were  available. 

An  attempt  is  made  to  relate  the  statistical  properties  of  the  bottom- 
backscattered  sound  field  to  the  independently  observed  geological  characteristics  of 
the  seafloor  surveyed.  Acoustic  boundary  mapping  over  flat  areas  is  achieved  by  fol¬ 
lowing  trend  changes  in  the  acoustic  data  both  along  and  across  track.  Such  changes 
in  the  acoustics  are  found  to  correlate  with  changes  in  bottom  type  or  roughness 
structure.  The  overall  energy  level  of  a  partial  angular-dependence  function  of  back- 
scattering  appears  to  depend  strongly  on  bottom  type,  whereas  the  shape  of  the  func¬ 
tion  does  not.  Clues  to  the  roughness  structure  of  the  bottom  are  obtained  by  relating 
the  shape  of  the  probability  density  function  of  normal-incidence  echo  envelopes  to  the 
degree  of  coherence  in  the  backscattered  acoustic  field. 
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IV.2.  INTRODUCTION 

In  the  past  few  years,  multibeam  echo-sounders  have  become  available  to  the 
scientific  community,  allowing  investigators  to  map,  with  high  resolution  and  in  near- 
real  time,  a  large  swath  of  seafloor  on  each  traverse  of  the  ship.  The  bathymetric 
charts  thus  obtained  represent  a  great  improvement  over  those  drawn  from  conven¬ 
tional  single-point  depth  recording  systems.  However,  bathymetry  reveals  the  shape 
of  seafloor  features  only  to  the  resolution  of  the  sounding  system;  it  does  not  yield 
other  seafloor  characteristics  such  as  bottom  type  or  bottom  microroughness  and  their 
respective  lateral  homogeneity. 

Because  bottom  roughness  and  variations  in  bottom  substrate  cause  fluctua¬ 
tions  in  the  backscattered  acoustic  signal,  such  seafloor  characteristics  can  be  inferred 
in  part  by  analyzing  the  structure  and  the  variations  of  this  signal  over  several 
transmission  cycles  (pings).  This  analysis  is  necessarily  statistical,  and  deals  with  an 
ensemble  of  independent  samples  since,  as  the  ship  moves,  each  ping  ensonifies  a 
slightly  different  portion  of  seafloor.  In  order  to  relate  the  statistics  of  the  backscat¬ 
tered  sound  field  to  the  characteristics  of  the  seafloor  from  which  it  emanates,  two 
main  theoretical  approaches  have  been  taken.  By  analogy  with  the  Huygens-Fresnel 
principle  of  physical  optics  [Born  and  Wolf,  1970],  both  approaches  consider  that  the 
sound  field  scattered  by  the  seafloor  consists  of  elementary  waves  in  mutual  phase 
interference.  Their  difference  lies  in  the  representation  of  the  irregularities  of  the 
seafloor.  In  the  first  approach  the  seafloor  is  a  surface  whose  irregularities  are 
described  by  a  root-mean  square  (rms)  roughness  in  the  vertical  dimension,  and  a 
correlation  distance  or  area  in  the  lateral  dimension.  The  statistics  of  the  surface  are 
then  incorporated  into  the  Helmholtz-Kirchhoff  formulation  of  the  scattering  theory  so 
that  statistics  of  the  sound  field  (usually  a  Gaussian  process)  can  be  used  to  estimate 
the  characteristics  of  the  surface  [Eckart,  1953;  Beckmann  and  Spizzichino,  1963; 
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Tolstoy  and  Clay,  1966;  Clay  and  Leong,  1974;  Clay  and  Medwin,  1977;  Bass  and 
Fuks,  1979;  Brekhovskikh  and  Lysanov,  1982].  The  second  approach  describes  the 
rough  surface  as  a  random  distribution  of  point  scatterers  reradiating  independently 
in  time.  It  is  the  quasiphenomenological  model  of  Faure  1964],  O’lshevskii  1967],  and 
Middleton  [1967],  which  uses  a  Poisson  process  giving  clues  to  the  density  of  the 
scatterers.  In  this  model,  however,  no  relation  exists  between  the  statistics  of  the 
scatterers  and  those  of  the  rough  surface.  For  this  reason,  the  first  approach  men¬ 
tioned  above  is  preferred  in  this  paper. 

In  the  Helmholtz-Kirchhoff  formulation  it  is  possible  to  determine  a  measure 
of  the  degree  of  coherence  in  the  backscattered  sound-field,  and  use  it  to  relate  the 
shape  of  the  probability  density  function  (pdf)  of  echo  envelopes  to  the  characteristics 
of  the  surface  (rms  roughness,  correlation  function)  [Stanton,  1984],  Because  of  their 
inherent  narrow  beam  geometry,  multibeam  echo-sounders  are  well  suited  for  such 
analysis.  In  addition,  backscatter  measurements  with  well  defined  spatial  resolution 
are  obtainable  both  along  and  across  the  ship’s  track.  A  description  of  seafloor 
characteristics  is  then  possible  in  terms  of  acoustic  boundaries. 

This  paper  presents  the  results  of  an  experiment  to  map  seafloor  acoustic 
backscattering  boundaries  with  a  Sea  Beam  multibeam  echo-sounder  operating  at  a 
frequency  of  12  kHz  with  a  7  ms  pulse  length  and  a  2  2/3  degree  angular  resolution. 
As  this  system  does  not  retain  the  acoustic  data  it  receives,  a  parallel  data  acquisition 
system  [de  Moustier,  1985al  built  by  the  Marine  Physical  Laboratory  (MPL)  of  the 
Scripps  Institution  of  Oceanography  was  used  to  record  digitally  the  envelope  of  the 
detected  echoes  over  a  variety  of  sea  floor  environments.  These  data  are  analyzed  as 
follows: 

(1)  First-order  statistics  (mean  and  variance)  of  peak  amplitude  in  the  near- 

specular  direction,  and  of  total  energy  for  nonspecular  beams  serve  to  quan- 
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tify  spatial  homogeneity  of  the  backscattered  sound  field. 

(2)  Angular  dependence  of  total  energy  is  used  to  complement  (1)  in  the  estima¬ 
tion  of  acoustic  boundaries  over  the  swath  of  seafloor  ensonified  by  the  Sea 
Beam 

(3)  Probability  density  functions  of  peak  amplitude  in  the  near-specular  direction 
are  compared  to  a  Rician  pdf  to  estimate  the  degree  of  coherency  of  the 
backscattered  signals,  yielding  a  qualitative  estimate  of  the  roughness  struc¬ 
ture  of  the  seafloor. 

However,  the  lack  of  phase  information,  occasional  saturation  in  both  the 
data  acquisition  and  the  Sea  Beam  systems,  and  sidelobe  interference  problems  limit 
the  scope  of  this  analysis.  Therefore  this  paper  only  intends  to  show  the  potential  for 
determining  seafloor  characteristics  that  exists  in  acoustic  backscatter  measurements 
derived  from  a  multibeam  echo-sounder.  Also,  acoustic  measurements  alone  are 
insufficient  to  determine  the  exact  nature  of  the  bottom.  Consequently  ground  truth 
must  be  obtained  by  independent  remote  sensing  (e.g.  deep  sea  photography  and/or 
television,  core  or  grab  samples,  etc.).  To  this  end  most  of  the  data  presented  here  are 
supported  by  at  least  one  independent  source  of  measurements. 

IV. 3.  THEORY 

In  this  section  we  review  the  theory  of  seafloor  acoustic  backscattering  using 
the  Helmholtz-Kirchhoff  formulation,  and  consider  its  applicability  to  backscatter 
measurements  made  with  a  Sea  Beam  system.  Next  we  derive  an  expression  for  the 
degree  of  coherence  of  the  backscattered  sound  field  valid  for  all  roughness  within  the 
limits  of  the  Kirchhoff  approximation.  This  is  an  extension  of  Stanton’s  [1984]  small¬ 
roughness  derivation,  which  relates  the  degree  of  coherence  in  the  field  to  the  rough¬ 
ness  and  correlation  structure  of  the  scattering  surface.  An  estimate  of  these  parame¬ 
ters  is  then  obtainable  from  the  shape  of  the  pdf  of  normal-incidence  echo  envelopes. 
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IV.3.1  Seafloor  acoustic  backscattering 

When  measuring  acoustic  backscatter  from  the  deep  seafloor,  the  first  com¬ 
mon  observation  is  that  individual  echoes  are  not  reproducible  and  that  there  can  be 
several  dB  of  variation  in  the  amplitude  of  the  returned  signal  from  one  ping  to  the 
next.  The  backscattering  process  is  therefore  considered  stochastic,  and  its  descrip¬ 
tion  needs  to  be  statistical. 

In  its  general  form,  the  narrow-band  acoustic  backscatter  is  expressed  as  a 
random  process  U{t)  such  that: 

17(0  — -B  cos  (w0t  +V0-  (!) 

where  E  is  the  signal  envelope,  ^  its  phase  and  u0  the  carrier  frequency.  All  the  infor¬ 
mation  concerning  the  random  function  U  is  therefore  contained  in  E  and  Vs  which  are 
random  functions  as  well. 

The  function  U  represents  the  output  of  the  receiving  hydrophone  array, 
commonly  a  voltage,  which  is  proportional  to  the  sound  pressure  of  the  bottom  echo 
impinging  upon  the  array.  Expressions  for  this  sound  pressure  have  been  derived  [Tol¬ 
stoy  and  Clay,  1966;  Clay  and  Leong,  1974;  Clay  and  Medwin,  1977;  Bass  and  Fuks, 
1979;  Brekhovskikh  and  Lysanov,  1982]  using  the  Helmholtz-Kirchhoff  formulation  in 
the  bistatic  scattering  from  a  rough  surface.  In  the  backscattering  case,  omitting  the 
time  dependence  and  propagation  losses,  a  receiver  at  Q  (Fig.  4.1)  will  sense  the  pres¬ 
sure: 


-KX3 

p(Q)  =  j'B1  exp  (2 jkR0)  Jf  D0  exp  [2 j(ax  +  eg  + 

—oo 


k 

2  R0 


(x2cos  20  +  y2))}dxdy 


with 


a 

c 

Bi 


=  k  sin  9 
=  — k  cos  0 
kBR 

2 7T  R %  cos  0 


(2) 


The  notation  follows  that  of  Clay  and  Medwin  [1977,  (Appendix  10)]  where  B  is 
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proportional  to  the  source  power,  D0  is  the  source  directivity,  k  is  the  acoustic  wave 
number  and  R  is  the  reflection  coefficient.  As  shown  in  Figure  4.1,  R0  is  the  range  to 
the  center  of  the  scattering  region  and  9  is  the  angle  of  incidence.  In  their  notation 
9  =0j  with  92=  —  8l  and  03=O.  a  and  c  are  the  horizontal  and  vertical  components  of 
wave  number.  is  the  normal  deviation  of  the  scattering  surface  from  the  xy  plane 
and  depends  on  position  in  the  plane:  $“  =  £(*,2/);  it  characterizes  the  roughness  of  the 
surface.  In  the  limit  corresponding  to  a  plane  surface,  the  return  is  a  specular 

reflection  at  normal  incidence,  and  outside  of  normal  incidence  Eq.  (2)  is  a  function  of 
the  beam  pattern  of  the  acoustic  array.  At  the  other  extreme,  when  £  is  large  the 
return  comes  mostly  from  acoustic  energy  scattered  by  the  rough  surface  back 
towards  the  receiver.  So  in  the  general  case  the  return  is  a  combination  of  scattered 
(incoherent)  and  reflected  (coherent)  energy  which  varies  depending  on  the  angle  of 
incidence  and  the  beam  pattern  of  the  array. 

Eq.  (2)  is  derived  in  the  far-field  of  the  transmit/receive  system  so  that 
ranges  R  can  be  approximated  by  R0  (Fig.  4.1)  in  the  expressions  of  the  incident  and 
backscattered  pressure  fields  except  in  the  phase  terms,  where  R  is  expanded  to 
second  order  terms  to  account  for  Fresnel  zone  contributions.  This  approximation  is 
valid  for  Sea  Beam’s  12  cm  wavelength  at  ocean  depths. 

The  Kirchhoff  approximation  is  also  used.  It  assumes  that  in  the  boundary 
conditions  on  the  surface,  the  reflection  coefficient  R  can  be  used  at  every  point  on  the 
rough  surface  by  approximating  the  field  at  any  point  on  the  surface  by  that  which 
would  be  present  on  the  tangent  plane  at  that  point  (Ch.  3  of  Beckmann  and  Spizzi- 
chino  [1963]).  This  requires  that  the  radius  of  curvature  of  the  irregularities  on  the 
seafloor  be  large  compared  to  the  acoustic  wavelength  or  that  no  shadowing  effects 
occur  within  the  ensonified  area.  These  conditions  are  met  in  most  of  the  data 
presented  here. 
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Figure  4,1  Geometry  of  backseattering.  Source  and  receiver  are  at  Q(x,z)  a  distance  R  away 
from  the  element  of  surface  ds.  (x,y)  is  the  elevation  of  ds  above  the  xy  plane. 
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Eckart’s  [1953]  small-slope  approximation  was  not  used  in  this  derivation. 
Had  this  approximation  been  used,  the  angular  dependence  in  the  fraction  of  Eq. 
(2)  would  appear  as  a  cos  9  in  the  numerator  instead  of  the  denominator.  The 
discrepancy  between  these  two  results  comes  from  the  contribution  of  the  horizontal 
component  of  wavenumber  (a),  which  is  neglected  in  the  small-slope  approximation. 

Eq.  (2)  also  assumes  that  the  area  ensonified  is  small  compared  to  R0)  so  that 
the  dependence  of  R0  on  x,y  or  c  can  be  ignored.  Likewise  R,  a  and  c,  which  depend 
mostly  on  the  angle  of  incidence  0,  are  assumed  to  have  only  small  variations  within 
the  scattering  area.  These  approximations  therefore  require  that  the  acoustic  system 
have  a  small  beam  width  and  commensurate  pulse  length.  Both  assumptions  are  rea¬ 
sonable  for  the  Sea  Beam  system. 

It  follows  that  for  a  given  angle  of  incidence  0 ,  the  random  character  of  p(Q) 
in  Eq.  (2)  is  mostly  due  to  the  fluctuations  of  the  phase  term  e2;cf,  which  accounts  for 
the  irregularities  of  the  bottom.  As  a  consequence,  the  bottom  roughness  $  is  con¬ 
sidered  a  random  process,  a  statistical  description  of  the  backscattered  sound  field  can 
be  achieved  by  ensemble-averaging  returns  from  successive  pings.  Such  returns  are 
independent  from  each  other  since  the  ship  moves  between  pings  and  a  slightly 
different  portion  of  seafloor  is  sampled  every  time.  For  an  ensemble  of  such  surfaces 
over  several  pings,  the  stochastic  function  ${x,y)  is  characterized  by  a  rms  roughness 
amplitude  (assuming  $  has  a  zero  mean) 

^  (3) 

a  "correlation  function": 

c{x\yi)  =  <${x>y)  y  +  y')  >  (4) 

and  a  pdf  W(^).  Eq.  (3)  and  (4)  describe  respectively  the  vertical  and  the  lateral 
dimensions  of  the  random  rough  surface.  (Note  that  in  the  statistics  literature,  Eq.  (4) 
is  referred  to  as  a  normalized  correlation  function  or  correlation  coefficient). 
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The  (ensemble)  average  pressure  is  then: 

-too 

<P>&~  j  pW{<;)d$  =  p0  <e2jcf>  (5) 

—OQ 

where  p„  is  the  reflection  from  a  plane  interface  (?  =  0  in  Eq.  (2)),  and  <e2,e  >  is 
defined  Papoulis,  1965]  as  the  characteristic  function  of  f. 

Assuming  f  to  be  normally  distributed  with  zero  mean,  Eq.  (5)  reduces  to: 

<P>E  =  Poe_2<:2'T2>  c2  =  fc2  cos2  9  (6) 

where  the  exponential  term  is  associated  with  the  degree  of  coherence  in  the  backscat- 

tered  sound  [Eckart,  1953;  Tolstoy  and  Clay,  1966].  For  large  cV,  <p>s  tends  to 
zero  and  the  backscatter  is  incoherent.  Conversely,  as  cV2  tends  to  zero  (plane  inter¬ 
face)  <p>e  tends  to  p0  and  the  backscatter  is  coherent. 

Eq.  (6)  makes  it  possible  to  estimate  the  rms  roughness  a  directly  from  meas¬ 
urements  of  the  coherently  backscattered  signals  <p>e ■  However,  although  this  is  a 
proven  technique  [Tolstoy  and  Clay,  1966;  Clay  and  Leong,  1974;  Clay,  1966;  Proud  et 
al.,  1960  ,  it  requires  measurements  of  the  phases  of  the  echo  signal,  which  are  not 
available  in  our  data.  For  this  reason,  we  resort  to  the  mean  square  pressure: 

<P2>e  =  <PP*>e 

where  p*  is  the  complex  conjugate  of  p. 

Substitution  of  Eq.  (2)  yields: 

OO 

<P2>e  «  Btffff  Dt<ex p  [2  jc (<;-<?)}>  dxdydJ di/ 

— OO 

with  Dx  =D0  Dj  exp  2j[a(x-x')  +  — ((x‘2-^c12)  cos2  9  +  y^-y^) ]  (8) 

The  average  within  the  integral  is  the  joint  characteristic  function  of  f  and  s'. 
Under  the  assumption  of  a  normally  distributed  surface  roughness,  the  joint  charac¬ 
teristic  function  becomes  [Papoulis,  1965]: 

< exp  2jc(f-d)>=  exp(— 4c2<t2(1— C)) 
where  G  is  the  correlation  function  defined  by  Eq.  (4). 


(9) 
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The  mean  square  pressure  <p2  >£,  which  includes  the  contribution  of  both  the 
reflected  (coherent)  and  scattered  (incoherent)  components  of  the  backscatter,  is 
therefore  related  to  the  statistics  of  the  rough  surface  through  the  rms  roughness  a 
and  the  spatial  correlation  function  C.  To  estimate  these  parameters  requires  the 
determination  of  the  respective  contributions  from  the  coherent  and  incoherent  com¬ 
ponents  to  the  total  backscattered  field,  which  is  equivalent  to  obtaining  a  measure  of 
the  degree  of  coherence  in  the  field.  To  this  end,  the  total  mean  square  pressure  is 
expressed  as  the  sum  of  the  coherent  <p  and  incoherent  <S2>%  components  [Tol¬ 
stoy  and  Clay,  1966,  p.  204] 

<p2>E=<p  >|+  <S2  >E  (10) 

where  <p  >|  =  <p  >E  <p*  >E,  so  that  <S2>E  is  simply  the  mean  square  second 
moment  of  p: 


■too 

<  S2  >  E  «  ms  dx  <exp2 jc(<;  —  ^)>  —  <exp2;cf  ><  exp(— 2jc  ^S^dxdydx1  dy* .  (11) 

-*oo  L  ■* 

It  follows  that  the  ratio  of  the  coherent  to  the  incoherent  components, 


<P  >1 

<S2>e’ 


(12) 


is  a  measure  of  the  degree  of  coherence  of  the  backscattered  field.  Other  authors 
[Novarini  and  Caruthers,  1972]  have  used  the  ratio  of  the  coherent  part  <p  to  the 
total  mean  square  pressure  <p2>  as  a  measure  of  the  degree  of  coherence.  Both 
expressions  are  valid,  but  Eq.  (12)  is  the  more  useful  in  this  context  since  it  gives  a 
link  between  the  statistics  of  the  backscattered  field  and  a  Rician  pdf  [Stanton,  1984]. 


To  express  7  as  a  function  of  a  and  C,  we  first  use  Clay’s  change  of  variable 
(Appendix  of  Clay  and  Medwin  1977 j  and  Clay  and  Leong  1974]): 


X  =x"  f/2 

y  =  y"  +V/2,  y'  =  y"  -v/2  (13) 

to  integrate  Eq.  (11)  over  x11  and  y ",  using  a  Gaussian  illumination  function  for  D0\ 
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D0  =exp(— — • 


jL, 


X2  Y2 


(14) 


where  X  and  Y  are  the  semi-minor  (respectively  major)  axes  of  the  ellipse  outlined  by 
the  intersection  of  the  mainlobe  of  the  beam  pattern  and  the  scattering  surface. 


For  beamwidths  Ax  and  A0  in  the  x  and  y  directions  respectively,  we  have 


x=*sm£2c  Y=Rs[n^,  (15) 

cosv 

where  0  is  the  angle  of  incidence  of  the  beam.  We  assume  spatial  stationarity  for  the 
surface  which  means  that  the  joint  characteristic  function  of  and  <*  depends  only  on 
the  distance  between  points  on  the  surface  through  f  =  x—x*  and  r\  —  y^y1  ,  and  that 
the  correlation  function  depends  only  on  £  and  tj:  C  =  C(Z,r}).  Then,  substituting  (13) 
and  (14)  into  (11)  and  integrating  over  xu  and  y!f  yield: 

*+00 

<S2  >E  =  B2f f D%(< exp2; c (f  —  ^)>  —  < exp2/c £> < exp(— 2;c {*)>) d£dr) 


with  B2 


BfirXY 

2 


B2  R2k2XY 

87rR%cos20 


D2  =  exp 


+  a£2+/?7/2) 


1 

a  =  — 

xtpcos^e 

\ 

+4- 

and  ft  =  — 

1^4  +-L) 

2 

R? 

X2  j 

2 

R?  Y2} 

We  then  obtain  an  expression  for  =  <P>e  using  Eq.  (2),  (5),  and  (14) 

+oo 

<p>§  —  B2  //  D2<exp  2jc  f>  <  exp(— 2/c  <f)>d £dr) 

—00 


(16) 


(17) 


In  keeping  with  the  assumption  of  a- normally  distributed  surface  roughness,  we  substi¬ 
tute  the  characteristic  function  by  its  expression  given  in  Eq  (6)  and  integrate  Eq  (17) 
to  obtain: 


<?>!  =  ■ 


B2  7T 

:  V^Jexp 

Likewise,  using  (6)  and  (9),  Eq.  (16)  becomes: 


.  2  2  a 


+00 


<S2  >e  =  B2  J JD2  exp (4eVC((,?)H 


d£dri 


(18) 


(19) 
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Hence  the  result: 


7 


-i 


<S2>e 

<p>i 


/ r  t  +0° 

=  exp  a2/a\  J J  exp  —(2 aj£  +  a£2  +  flr}2)\  exp  f 4 c 2cr2 C7(^,r7 )1 —  1  d£dr}  (20) 

7T  L  J  _o0  }  \  ) 

which  relates  the  degree  of  coherence  of  the  backscattered  field  to  the  parameters  of 
the  rough  surface  <r  and  G. 


This  result  is  general  in  the  sense  that  no  limitations  have  been  imposed  on 
the  wave  length  of  the  sound  radiation  or  equivalently  on  the  length  scale  of  the  sur¬ 
face  roughness.  We  have  only  assumed  that  the  surface  satisfies  the  Kirchhoff  boun¬ 
dary  condition,  which  requires  there  be  no  sharp  edges  on  the  scattering  surface.  We 
have  also  assumed  that  the  roughness  is  normally  distributed  and  that  it  is  spatially 
stationary  by  casting  the  correlation  function  C  as  a  function  of  the  distance  between 
points  on  the  surface  ((7(£,t?)).  Although  the  applicability  of  such  statistical  properties 
to  the  ocean  floor  may  seem  questionable,  especially  with  regard  to  stationarity,  they 
are  useful  in  reducing  the  foregoing  integrals  to  more  manageable  expressions. 

Outside  of  normal  incidence,  the  coherent  component  of  backscattering  is 
likely  to  be  small  or  negligible  compared  to  the  incoherent  component,  so  that  7  will 
tend  to  zero.  Normal  incidence  is  more  interesting  in  that  both  components  are  then 
equally  likely  to  dominate,  depending  on  the  type  of  rough  surface;  and  there  is  a 
direct  relationship  between  7  and  the  shape  of  the  pdf  of  normal  incidence  echo 
envelopes. 


At  normal  incidence  (0  =0),  Eq.  (20)  becomes: 


7“*  =  //ex  p  -  (af2  +  @rj2)]  I 

—CXI  L  I 


7T 


exp 


4 C((,v) 


d(drj 


(21) 


To  integrate  Eq.  (21)  requires  the  determination  of  an  analytical  form  of  the  correla¬ 
tion  function.  Although  several  forms  of  this  function  have  been  used  in  rough-surface 
scattering  theory  [Horton,  1972],  most  of  them  are  one-dimensional  and  are  restricted 
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to  ripple-like  features.  Stanton  [1984]  introduced  a  correlation  area  by  considering  a 
two-dimensional  correlation  function  which  is  zero  on  the  average  except  for  a  small 
"elliptic  cone"  around  f  =  rf  =0  where  it  decreases  linearly  from  1  at  £  =  77  =0.  His 
application  was  a  surface  consisting  of  ripples,  but  the  approach  is  also  valid  for  other 
types  of  roughness  (nodule  fields,  bed  rocks,  etc.)  more  representative  of  our  data.  For 
ease  of  integration  we  prefer  to  use,  as  a  two-dimensional  correlation  function,  an 
elliptic  paraboloid  such  that: 


C(s,v)  = 


0 


ii 

/? 


otherwise 


(22) 


In  the  limit  (4£2cr2<Cl)  of  small  to  moderate  roughness  the  second  exponential  in  Eq. 
(21)  is  expandable  in  a  Taylor  series.  Keeping  the  first  two  terms  of  the  expansion 
and  using  Eq.  (22),  Eq.  (21)  integrates  into: 


7-1  2s  —  ^*2-  kWlxl2  (23) 

3  7 r 

where  lx  and  /2  are  the  correlation  distances  beyond  which  (7(£,??)  =0. 

It  has  been  assumed  that  for  \$<t\  and  \t]\<12  the  contribution  of  the  first 
exponential  term  in  Eq.  (21)  can  be  neglected.  This  is  equivalent  to  requiring  the 
exponential  term  to  reach  e-1  outside  of  the  range  |f| <lY  \r]\<l2l  so  that  a£2  +  ftrj2  =  1 
for  \^Xi  and  l^>/2.  Except  for  a  factor  of  2,  Eq.  (23)  is  equivalent  to  Stanton’s  Eq. 
(16)  1984].  It  shows  that  for  small  to  moderate  roughness,  the  degree  of  coherence  in 
the  specular  backscattered  sound  field  depends  on  the  relative  roughness  (£cr)  and  the 
correlation  area  (/x/2)  of  the  surface  as  well  as  on  the  beam  widths  of  the  measuring 
system  (a,/?).  Because  the  same  degree  of  coherence  is  a  parameter  which  controls  the 
shape  of  the  pdf  of  normal  incidence  echo  envelopes,  Eq.  (23)  makes  it  possible  to 
evaluate  the  parameters  of  the  rough  surface  directly  from  the  shape  of  the  pdf.  Also, 
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pending  determination  of  the  rms  roughness  from  coherent  measurements  (Eq.  (6)),  the 
correlation  area  is  readily  obtained  from  Eq.  (23). 


IV.3.2  Envelope  distributions 

As  summarized  in  Appendix  A,  the  Rician  distribution  [Rice,  1945]  is  useful  in 
describing  the  effects  of  coherent  scattering  on  the  shape  of  the  pdf  of  echo  envelopes. 
By  considering  the  backscattered  pressure  wave  as  a  sum  of  a  coherent  component  of 
rms  value  A  and  an  incoherent  component  normally  distributed  with  variance  cr^  the 
parameter  7  can  be  written  as: 


7  = 


(24) 


It  follows  that  the  Rician  distribution  for  the  corresponding  echo  envelopes  E  (Eq.  1) 
can  be  expressed  as  a  function  of  E,  its  second  moment  <E2>}  and  7: 


W(E)  =(2+7) 


E 


<E*> 


exp 


(2+7)E2  ±  7  <E2  > 


2  <E2  > 


Io 


E 


<E 2  >* 


(7(2+7))* 


(25) 


In  the  limit  where  7  tends  to  zero  (dominance  of  the  incoherent  component),  Eq.  (25) 
tends  to  the  Rayleigh  distribution  (Eq.  A. 4).  Conversely,  when  the  coherent  com¬ 
ponent  dominates  (7^1)  Eq.  (25)  tends  to  a  Gaussian-like  distribution.  A  qualitative 
estimate  of  the  degree  of  coherence  in  the  backscattered  sound  field  can  then  be 
obtained  directly  from  the  value  of  7  that  gives  the  best  fit  of  the  Rician  distribution 
to  a  histogram  of  the  normal-incidence  echo-envelope  peak  distribution. 


This  is  analogous  to  the  determination  of  the  contrast  of  a  speckle  pattern 
[Goodman,  1975]  in  laser  optics.  Speckle  patterns  are  the  result  of  constructive  and 
destructive  interference  between  light  waves  scattered  from  elementary  areas  of  a 
rough  surface.  Their  contrast  is  analogous  to  7"1  (Eq.  21)  since  an  ideal  mirror  surface 
does  not  produce  a  speckle  pattern  (contrast  =0). 

As  in  optics,  useful  quantitative  results  concerning  the  roughness  structure  of 


97 


the  surface  are  obtained  for  acoustic  backscatter  only  when  the  roughness  of  the  sur¬ 
face  is  a  fraction  of  the  acoustic  wavelength.  In  such  cases,  the  value  of  7  derived 
from  the  shape  of  the  pdf  can  be  directly  related  to  the  product  of  the  rms  roughness 
and  the  correlation  area  of  the  surface  (Eq.  23).  For  practical  purposes,  this  technique 
is  applicable  when  7>1  because  there  is  little  difference  in  the  shape  of  the  pdf’s  for 
0<7<1  (e.g.  Figure  2  of  [Stanton,  1984]). 

IV.4.  SEA  BEAM  ACOUSTIC  DATA. 

To  explain  some  of  the  peculiarities  of  the  Sea  Beam  data  presented  here,  a 
brief  description  of  the  acoustic  geometry  is  necessary. 

TV.4.1  Acoustic  geometry 

Sea  Beam  uses  two  arrays  perpendicular  to  each  other.  The  transmit  array 
consists  of  twenty  projectors  aligned  along  the  ship’s  keel.  The  outputs  of  the  projec¬ 
tors  are  amplitude-shaded  for  sidelobe  control  according  to  the  Dolph-Chebyschev 
method  [Dolph,  1946],  and  phase  corrected  to  ensure  vertical  projection  at  all  times. 
The  resulting  beam  pattern  is  2  2/3  degrees  wide  in  the  fore-aft  direction  and  54 
degrees  athwartships  (Fig.  4.2a).  These  angles  correspond  to  the  half  power  point  of 
the  beam  patterns.  This  convention  will  be  kept  in  the  following. 

The  receive  array  consists  of  forty  elements  arranged  athwartships  in  a  V- 
shaped  array  symmetrical  about  the  ship’s  keel.  The  outputs  of  these  elements  are 
amplitude-shaded  (Dolph-Chebyschev  shading)  for  sidelobe  control,  and  the  resulting 
beam  is  steered  phase-wise  to  form  sixteen  preformed  beams  spaced  2  2/3  degrees 
apart  between  +20  and  -  20  degrees  of  incidence.  Each  beam  is  2  2/3  degrees  wide  in 
the  athwartships  direction  and  20  degrees  fore-aft  (Fig.  4.2b).  The  comparatively 
wide  beam  width  in  the  fore-aft  direction  is  meant  to  accommodate  pitch  angles  of  ± 
10  degrees  for  no  pitch  correction  is  performed  on  the  receive  array. 


a)  TRANSMIT  b)  RECEIVE  c)  TRANSMIT/RECEIVE 
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Figure  4.3  Preformed  beams.  Four  out  of  the  sixteen  preformed  beams  are  displayed  to  show 
how  sidelobes  from  each  beam  "look"  into  the  direction  of  the  mainlobe  of  all  the  others.  The 
sidelobe  contributions  to  the  sixteen  bottom  echoes  received  each  ping  can  be  seen  in  Figure  4.4. 
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Note  also  that  each  of  the  preformed  beams  has  sidelobes  pointing  in  the  direction  of 
(but  not  necessarily  aligned  with)  the  mainlobe  of  all  the  other  beams  (Figure  4.3).  As 
will  be  shown  in  the  following,  this  fact  is  significant  when  dealing  with  echo  envelopes 
of  the  sixteen  preformed  beams.  For  a  broader  description  of  the  Sea  Beam  system, 
the  reader  is  referred  to  a  comprehensive  review  by  Renard  and  Allenou  [1979]. 

IV.4.2  Sea  Beam  digitized  echo  envelopes 

Figure  4.4  illustrates  a  typical  set  of  echo  envelopes  as  they  appear  at  the 
output  of  Sea  Beam’s  echo-processor  receivers.  Each  envelope  corresponds  to  the 
return  on  one  of  the  sixteen  preformed  beams  numbered  1  through  8  from  the  center 
of  the  ship  out,  on  both  port  and  starboard.  These  data  have  not  been  roll- 
compensated;  therefore  the  ship’s  center  line  does  not  necessarily  lie  on  the  true  verti¬ 
cal  The  ridge  of  synchronous  returns  corresponds  to  energy  from  a  strong  return  in 
the  near-specular  direction  leaking  into  the  sidelobes  of  all  the  other  beams  (e.g.  star¬ 
board  beam  ^  1  in  this  figure).  The  term  near-specular  is  used  here  to  indicate  that 
because  Sea  Beam  receives  with  discrete  beams,  the  return  may  be  within  2  2/3 
degrees  of  the  specular  direction.  Likewise,  in  the  following  nonspecular  beams  refer  to 
those  beams  which  are  not  within  2  2/3  degrees  of  the  specular  direction.  In  its 
current  mode  of  operation,  the  Sea  Beam  echo  processor  digitizes  these  echo  envelopes 
and  applies  ray-bending,  roll,  gain  and  sidelobe  corrections.  A  time  of  arrival  is  then 
determined  at  the  center  of  mass  of  each  of  the  corrected  echoes  and  is  converted  into 
a  depth  and  a  cross-track  distance.  These  are  in  turn  logged  on  magnetic  tape  as  well 
as  recorded  on  paper  as  a  bathymetric  contour  chart.  No  further  use  is  made  of  the 
echo  signals,  which  are  then  discarded.  In  an  effort  to  preserve  these  echo  signals  for 
analysis,  MPL  built  a  parallel  acoustic  data  acquisition  system  to  record  digitally 
data  such  as  that  of  Figure  4.4  along  with  a  time-varied  gain  and  the  ship’s  roll 
angle.  In  this  paper  some  of  the  envelope  data  collected  with  the  MPL  system  suffer 
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Figure  4.4  Acoustic  signal  envelopes  of  the  sixteen  preformed  beams  at  the  output  of  Sea 
Beam’s  Echo  Processor  Receivers .  Time  is  in  seconds  after  transmission.  Signal  amplitude  is 
in  volts  corrected  for  acoustic  transmission  loss  by  a  time-varied  gain.  In  this  display,  the 
data  have  not  been  corrected  for  ship  roll  or  receiver  gain.  The  ridge  of  synchronous  returns  is 
the  sidelobe  response  to  a  strong  return  in  the  near-specular  direction  (STBD  beam  #  1  in  this 
Figure). 
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amplitude  clipping  due  to  saturation  in  both  the  data  acquisition  system  and  the  Sea 
Beam  echo-processor  receivers.  The  saturation  in  the  receivers  occasionally  disrupts 
Sea  Beam’s  echo  processing,  resulting  in  spurious  depth  measurements.  Likewise, 
because  the  echo  processing  is  done  on  the  envelope  of  the  signals,  the  system  cannot 
differentiate  between  sidelobe  response  and  bottom  return  when  the  two  overlap  (e.g. 
beam  #  2  in  Fig.  4.4).  Improper  depth  determination  and,  bathymetric  artifacts 
result.  Such  bathymetric  artifacts  have  been  analyzed  in  detail  elsewhere,  [de  Mous- 
tier  and  Kleinrock,  in  press]  and  we  shall  not  repeat  the  discussion  here.  Nevertheless 
this  sidelobe  interference  constitutes  a  serious  hindrance  when  analyzing  signal 
envelopes  where  sidelobe  and  bottom  return  overlap.  To  tell  them  apart  requires 
phase  information  which  is  not  available  in  our  present  data  set. 

IV. 4. 3  Data  reduction 

Owing  to  the  limitations  outlined  above,  the  acoustic  data  recorded  from  Sea 
Beam  were  reduced  to  be  analyzed  in  two  ways.  The  first  approach  concentrates  on 
the  statistics  of  the  peak  amplitude  in  the  near-specular  direction;  the  second  deals 
with  both  angular  dependence  and  statistics  of  total  energy  in  nonspecular  beams.  To 
overcome  the  saturation  problem  mentioned  above,  the  mean  sidelobe  response  was 
used  to  retrieve  the  peak  amplitude  of  the  near-specular  returns  that  appeared 
clipped.  The  method  assumes  that  the  ratio  of  the  near-specular  peak  amplitude  to 
the  corresponding  mean  sidelobe  response  is  approximately  constant  from  ping  to  ping, 
and  that  most  of  the  variability  in  this  ratio  is  due  to  (I)  the  slight  misalignment  of 
sidelobes  with  the  mainlobe  of  the  near-specular  beam  (Figure  4.3)  and  (2)  the  roll  of 
the  ship.  For  simplicity,  only  data  collected  over  nearly  flat  seafloor  (slope  angles  less 
than  2  2/3  degrees)  were  used  in  this  exercise,  for  owing  to  the  acoustic  geometry,  no 
specular  backscatter  is  to  be  expected  from  bottoms  sloping  up  or  down.  To  avoid 
bias  from  bottom  returns,  the  mean  sidelobe  response  was  computed  as  the  arithmetic 
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mean  of  the  peak  values  of  sidelobe  contributions  well  separated  from  the  bottom 
returns.  As  an  example,  beams  number  port  8-3  and  starboard  4-8  would  qualify  in 
Figure  4.4.  By  working  on  non-clipped  data,  this  mean  was  then  compared  with  the 
amplitude  of  the  corresponding  near-specular  peak  by  computing  their  ratio. 

From  an  ensemble  of  over  four  hundred  such  ratios  the  mainlobe  to  sidelobe 
ratio  is  17.80  (25  dB)  with  a  standard  deviation  of  2.14.  This  value  is  in  agreement 
with  what  Dolph  [1946]  predicted  for  his  amplitude  shading  method.  It  is  however  2 
to  3  dB  lower  than  the  level  measured  by  Renard  and  Allenou  [1979]  (27-28  dB)  on  a 
different  Sea  Beam  system.  But  because  their  measurements  were  done  on  only  two 
preformed  beams,  it  is  reasonable  to  expect  the  level  obtained  by  averaging  over  ten 
beams  to  be  lower. 

The  ping-to-ping  variations  of  the  near-specular  peak  were  then  inferred  from 
the  variations  of  the  corresponding  mean  sidelobe  response  using  the  same  method. 
When  the  sidelobe  response  had  insufficient  signal-to-noise  ratio,  the  near-specular 
return  was  not  clipped  so  its  peak  amplitude  was  used  to  compute  the  mean  sidelobe 
level  by  subtracting  25  dB.  Results  of  this  method  are  illustrated  in  Figure  4.5. 

Because  of  the  time  stretching  evidenced  on  the  non  specular  returns  (Fig. 
4.4),  their  description  is  more  appropriately  based  on  total  energy  than  on  peak  ampli¬ 
tude.  To  establish  a  common  reference  between  pings,  the  sixteen  beams  are  roll- 
compensated  and  put  into  a  set  of  beams  2  2/3  degrees  apart,  one  of  which  is  centered 
on  the  true  vertical.  Each  roll-compensated  beam  is  obtained  by  linear  interpolation 
between  the  two  adjacent  unstabilized  beams.  Ideally,  a  ray-bending  correction 
should  also  be  applied.  However,  errors  resulting  from  the  omission  of  this  correction 
are  small  at  the  steep  angles  of  incidence  considered  here  and  are  inconsequential  in 
the  scope  of  this  paper.  Again  roll-compensated  beams  for  which  sidelobe  response 
and  bottom  return  overlap  are  disregarded.  The  total  energy  in  a  return  is  then 
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Figure  4.5.  Variations  in  the  amplitude  of  the  near-specular  return  which  is  clipped  in  the 
recorded  data  (a,  top)  are  recovered  from  the  corresponding  mean  sidelobe  response  (a,  bot¬ 
tom).  Each  data  point  corresponds  to  one  ping.  The  cross  covariance  between  these  two  time 
series,  excluding  the  clipped  portions,  is  shown  in  (b)  where  successive  lags  are  successive  pings. 
The  mean  sidelobe  response  was  found  to  be  25  dB  below  the  near-specular  return. 
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computed  as  the  difference  between  the  mean  square  amplitude  (signal  +  noise)  in  a 
predetermined  window  and  the  mean  square  noise  in  an  equivalent  window.  This  way 
a  fixed  window  size  can  be  used  for  all  beams. 

As  no  accurate  calibration  exists  for  the  Sea  Beam  system  from  which  this 
acoustic  data  was  recorded,  the  following  results  are  given  in  relative  units. 

IV.5.  RESULTS  AND  DISCUSSION. 

Given  that  variations  in  bottom  type  and  bottom  roughness  are  the  main 
factors  in  the  nature  of  the  fluctuations  of  the  acoustic  backscatter  from  the  seafloor, 
data  from  three  geologically  different  areas  have  been  analyzed  to  determine  what 
could  be  learned  from  the  acoustics.  Although  it  is  possible  to  tell  these  areas  apart 
from  analysis  of  the  acoustic  data  alone,  our  acoustic  data  base  is  too  small  and 
incomplete  to  identify  them,  and  we  use  ground  truth  from  independent  measurements 
to  validate  the  acoustic  measurements  made  with  Sea  Beam. 

The  three  areas  investigated  are:  (1)  a  manganese  nodule  field  in  the  North 
Eastern  Tropical  Pacific  (15 0  N,  125  0  W)  for  which  there  are  bottom  photographs 
taken  with  MPL’s  Deep  Tow  instrument  package  [Spiess  and  Lonsdale,  1982],  as  well 
as  box  core  data;  (2)  a  sedimentary  environment  in  the  North  San  Clemente  basin  (~ 
150  km  south-west  of  San  Diego,  Calif.)  with  Deep  Tow  bottom  photographs  in  the 
general  area;  and  (3)  a  lava  sheet  flow  on  the  crest  of  the  East  Pacific  Rise  around 
10  °N  with  supporting  data  from  the  Lamont-Doherty  Geological  Observatory  bea 
MARC  I  side-looking  sonar  [Kosalos  and  Chayes,  1983]  and  bottom  photographs. 

For  each  area,  a  composite  figure  (Figures  4.6,8,10)  has  been  assembled  to 
give  a  synopsis  of  the  Sea  Beam  data.  Sea  Beam  bathymetry  is  shown  in  a  three- 
dimensional  view  of  a  single  swath,  as  a  mesh  of  instantaneous  cross  track  depth 
profiles,  low-pass  filtered  along  track  and  displayed  with  a  vertical  exaggeration  of  10. 
The  acoustic  data  are  presented  as  along-track  variations  of  relative  peak  amplitude 
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in  the  specular  direction  and  relative  total  energy  for  nonspecular  beams.  As 
described  in  Section  IV. 4,  the  peak  amplitude  information  has  been  recovered  from 
sidelobe  data,  and  the  energy  has  been  calculated  only  for  nonspecular  beams  where 
bottom  return  and  sidelobe  response  were  well  separated.  In  all  displays,  the  data 
have  been  low-pass  filtered  with  a  running  mean  averaging  over  the  number  of  pings 
necessary  to  traverse  a  single  beam  footprint  (2  2/3  x  2  2/3  degree)  on  the  seafloor  at 
the  depth  considered.  A  coefficient  of  variation,  which  is  simply  the  standard  devia¬ 
tion  normalized  by  the  mean  over  the  averaging  interval,  is  also  displayed  to  give  a 
measure  of  the  variability  in  the  data.  Finally,  partial  angular  dependence  profiles  of 
total  energy  are  stacked  and  low-pass  filtered  along  track  to  show  trends  in  the  back- 
scattered  acoustic  energy  both  along  and  across  track.  The  central  portion  of  the 
angular  dependence  of  total  energy  has  been  purposely  left  out  because  of  our  inability 
to  remove  the  sidelobe  contribution  contaminating  the  returns  closest  to  the  specular 
direction  without  seriously  degrading  the  returns  themselves. 

In  the  following  discussion,  we  assess  how  much  can  be  learned  about  a  por¬ 
tion  of  seafloor  surveyed  with  Sea  Beam  by  analysis  of  the  acoustic  backscatter  it 
receives.  We  consider  successively  the  system’s  ability  to  delineate  acoustic  boun¬ 
daries,  the  use  of  an  angular  dependence  function  of  backscattering  to  differentiate 
between  various  types  of  substrate,  and  the  potential  for  estimating  the  microrough¬ 
ness  of  the  bottom. 

IV. 5.1  Mapping  acoustic  boundaries 

Because  most  statistical  analyses  of  the  backscattered  sound  field  are  based 
upon  the  assumption  of  a  homogeneous  (stationary)  scattering  surface,  it  is  important 
to  be  able  to  isolate  seafloor  areas  for  which  this  assumption  holds.  To  this  end  we 
define  an  acoustic  boundary  as  the  place  where  a  marked  change  in  trend  appears  in 
the  acoustic  data.  With  a  multibeam  system  such  trends  can  be  followed  both  along 
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and  across  tracks,  given  that  it  is  possible  to  correct  for  bottom  slope  in  both  direc¬ 
tions.  In  this  paper,  we  simplify  the  problem  by  limiting  ourselves  to  nearly  flat  por- 
tions  of  seafloor. 

IV. 5. La  Manganese  nodule  area: 

In  a  previous  paper  ide  Moustier,  1985b],  we  used  the  variations  in  amplitude 
of  the  specular  beam  alone  to  infer  manganese  nodule  coverage  over  a  well  docu¬ 
mented  nodule  mining  site.  Although  the  estimates  of  coverage  were  crude,  they  were 
in  very  good  qualitative  agreement  with  estimates  of  coverage  derived  from  Deep  Tow 
bottom  photographs  of  the  same  area.  Our  ability  to  correctly  identify  bare  mud 
patches,  areas  sparsely  covered  with  nodules,  and  areas  densely  covered  throughout 
the  mining  field  was  a  good  indication  that  acoustic  boundary  mapping  is  feasible  with 
Sea  Beam.  Here  we  extend  the  analysis  to  nonspecular  beams. 

As  shown  in  Figure  4.6,  the  portion  of  seafloor  considered  is  nearly  flat,  and 
the  average  depth  is  about  4500  m.  At  this  depth,  the  Sea  Beam  transmits  every  8 
seconds  and,  since  the  ship  was  moving  at  —  5  m/s  (10  knots),  the  spatial  sampling 
interval  is  roughly  40  m  along  track.  By  comparison,  the  diameter  of  the  vertical 
incidence  footprint  is  about  200  m,  so  that  it  takes  five  pings  to  traverse  a  footprint. 
This  number  was  used  as  the  averaging  interval  in  the  low-pass  filter  of  the  data  for 
the  area. 

With  this  averaging,  the  peak  amplitude  of  normal  incidence  (0  )  returns  is 
seen  to  fluctuate  around  a  mean  level  of  about  15  dB  between  0  and  6  km  and  then 
drop  sharply  beyond  the  6  km  mark  (Fig.  4.6).  A  similar  trend,  although  much  better 
defined,  is  seen  in  the  mean  total  energy.  For  reference,  ±  13  degrees  incidence 
corresponds  to  a  distance  across  track  of  ±  1  km  from  the  vertical  incidence  point. 
The  drop  in  level  is  clearly  seen  in  both  port  and  starboard  plots,  but  it  happens  ~  1 
km  further  along  track  (7  km  mark)  on  the  starboard  side.  It  can  be  followed  also  on 
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Figure  4.6  Manganese  nodule  area  synopsis  of  Sea  Beam  data.  Upper  left:  bathymetry  for 
one  swath  of  Sea  Beam.  The  depth  scale  in  meters  is  shown  with  a  vertical  exaggeration  of  10 
and  its  origin  is  at  4600m.  The  distance  across  and  along  track  are  in  kilometers.  Across 
track,  a  vertical  bar  indicates  the  position  of  normal  incidence.  The  terrain  has  a  very  gentle 
undulation  along  track.  Upper  right:  along-track  variations  of  peak  amplitude  in  the  near- 
specular  direction  (~  0  incidence)  and  of  total  energy  at  about  13“  incidence.  Plots  have 
been  low-pass  filtered  by  averaging  over  five  pings.  The  coefficient  of  variation  (standard 
deviation/ mean)  shows  the  variability  in  the  data  for  each  of  the  three  angles  of  incidence. 
Lower:  partial  angular  dependence  of  total  energy  and  its  variations  along-track  are  displayed 
in  both  left  and  right  view.  The  center  portion  (between  ±5°  1/3  incidence)  has  been  left 
out  because  of  sidelobe  interference  and  saturation  in  the  data.  The  notch  seen  on  the  outer 

starboard  side  is  probably  a  system-related  artifact,  for  it  is  found  to  some  extent  in  all  the 
recorded  data. 
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Figure  4.7  Representative  Deep  Tow  bottom  photographs  taken  in  the  manganese  nodule  field. 
They  show  the  contrast  between  a  relatively  dense  coverage  of  manganese  nodules  (a)  and 
bare  mud  (b). 
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the  stacked  profiles  of  total  energy,  where  it  is  seen  to  run  diagonal  to  the  ship’s 
track.  This  diagonal  therefore  defines  an  acoustic  boundary. 

Deep  Tow  bottom  photographs  taken  along  the  same  track  show  that  this 
boundary  is  associated  with  a  rapid  transition  between  dense  nodule  coverage  (Fig. 
4.7a)  and  bare  mud  (Fig.  4.7b).  In  this  case,  the  change  in  bottom  substrate  from 
nodules  to  mud,  rather  than  their  relative  roughness,  seems  to  be  the  dominant  factor 
in  the  backscattering  process.  Also,  it  is  interesting  to  note  that  in  this  area,  nonspec- 
ular  beams  energy  is  better  suited  for  acoustic  boundary  mapping  than  peak  ampli¬ 
tude  in  the  specular  return.  The  peak  amplitude  data  shown  in  Figure  4.6  seems  to  be 
more  sensitive  to  small  variations  in  bottom  slope  and  to  potential  bottom  focusing 
effects,  and  would  therefore  require  more  averaging  to  bring  out  the  underlying  trend 
readily  observable  in  the  energy  data  of  nonspecular  beams. 

IV.S.l.b  North  San  Clemente  basin 

The  portion  of  Sea  Beam  data  considered  in  Figure  4.8  corresponds  to  a  flat 
area  around  32  30'N,  118  10'W.  The  bathymetry  shows  a  small  diagonal  step  rising 
10  m  over  150  m,  and  the  average  depth  is  1900  m.  At  this  depth  Sea  Beam  transmits 
every  4  seconds.  In  this  instance,  the  ship  speed  was  about  1  m/s  (~  2  knots),  yielding 
a  sampling  interval  along  track  of  roughly  4  m,  so  that  22  pings  are  necessary  to 
transverse  the  length  of  a  vertical  incidence  footprint  (88  m).  This  number  of  pings 
was  therefore  used  as  the  averaging  interval  to  low-pass  filter  the  data. 

Because  this  averaging  interval  is  large,  the  profile  of  normal  incidence  peak 
amplitude  data  appears  relatively  smooth,  with  less  variability  than  that  of  Figure 
4.6.  No  major  trend  changes  are  observable  in  this  profile  (Fig.  4.8).  The  drop  in 
amplitude  (arrow)  around  the  1  km  mark  is  most  likely  due  to  the  small  bathymetric 
step.  The  uniformity  of  this  portion  of  seafloor  is  confirmed  by  the  mean  total  energy 
data  even  though  a  small  undulation  with  a  1.5  km  wavelength  is  noticeable  on  both 
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Figure  4.8  North  San  Clemente  Basin:  synopsis  of  Sea  Beam  data.  Layout  and  scales  are  the 
same  as  in  Figure  4.6.  The  origin  of  the  depth  scale  is  at  2000  m.  Depth  is  shown  with  a  verti¬ 
cal  exaggeration  of  10.  A  small  step  (~  10  m  high)  runs  diagonal  to  the  swath  around  the  1  km 
mark  along-track.  Plots  of  amplitude  and  total  energy  data  have  been  low-pass  filtered  by 
averaging  over  22  pings.  The  small  step  in  the  bathymetry  appears  in  the  amplitude  data  (ar¬ 
row)  as  a  marked  dip  in  amplitude.  The  gap  in  the  center  portion  of  the  angular  dependence  of 
total  energy  is  larger  in  this  figure  than  in  Figure  4.6  because  as  the  water  depth  decreases  so 
does  the  time  separation  between  arrivals  on  individual  beams.  As  a  result,  the  sidelobe  in¬ 
terference  affects  more  beams. 
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Figure  4.9  Deep  Tow  bottom  photograph  taken  in  the  North  San  Clemente  Basin.  The  sedi¬ 
ments  appear  relatively  smooth  with  occasional  animal  generated  boreholes  and  craters,  and 
numerous  worm  tracks. 
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port  and  starboard.  This  trend  is  most  marked  at  13  degrees  incidence  to  port.  Dis¬ 
tances  across  track  are  440  m  and  545  m  from  the  vertical  incidence  point  at  13  and 
16  degrees  incidence  respectively.  As  seen  in  the  stacked  profiles  of  total  energy,  the 
undulations  on  port  and  starboard  are  not  symmetric  with  respect  to  vertical 
incidence,  but  appear  to  be  offset  diagonally  from  each  other  in  roughly  the  same 
orientation  as  the  bathymetric  step.  A  tenuous  acoustic  boundary  can  then  be  defined 
on  this  Sea  Beam  swath  as  a  patch  of  higher  backscatter  extending  1.5  km  along 
track  and  trending  diagonally  across  track. 

Owing  to  the  complex  nature  of  sedimentation  patterns  in  this  area,  it  is 
difficult  to  relate  the  patchiness  observed  in  the  acoustics  to  geological  processes.  Sed¬ 
iments  of  the  San  Clemente  basin  are  of  both  turbidite  and  pelagic  origin  and  contain 
mostly  fine  grained  sand  and  muds  (biogenic  and/or  micaceous)  [Normark  and  Piper, 
1972;  Emery,  I960].  Deep  Tow  bottom  photographs  taken  in  the  vicinity  (Fig.  4.9) 
show  a  smooth  sedimentary  bottom  with  evidence  of  intense  animal  activity  (worm 
tracks  and  feces,  holes  and  craters,  etc.)  so  that  bioturbation  must  play  a  major  role 
in  the  vertical  distribution  of  sediments.  Several  speculations  can  be  made  to  account 
for  the  acoustic  boundary  observed.  It  could  be  due  to  a  patch  of  seafloor  where  sand 
has  been  bioturbated  with  the  overlying  mud,  thereby  enhancing  its  backscattering 
properties.  A  thinning  of  the  surficial  mud  layer  over  sand  would  have  a  similar  effect. 
Another  possibility  is  a  change  in  the  fine-scale  roughness  of  the  bottom  as  a  result  of 
animal  activity,  higher  backscattered  energy  corresponding  to  a  rougher  interface. 
More  data  (subbottom  profiles  and/or  cores)  are  necessary  to  determine  whether 
roughness  or  bottom  type  or  a  combination  dominates  the  backscattering  process  in 


this  area. 


Figure  4*10  Lava  sheet  flow  on  the  crest  of  the  East  Pacific  Rise:  synopsis  of  Sea  Beam  data. 
Layout  and  scales  are  identical  to  those  of  Figures  4.6  and  4.8.  The  origin  of  the  depth  scale  is 
at  2700  m.  The  acoustic  data  displayed  correspond  to  the  section  of  bathymetry  (~  2.7  km 
long)  delimited  by  the  arrows.  In  this  section,  the  bathymetry  is  nearly  flat.  For  compactness 
of  this  figure,  the  origin  (0)  of  the  vertical  axes  for  the  acoustic  data  is  not  shown;  however  the 
scale  is  the  same  as  in  Figures  4.6  and  4.8. 
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Figure  4.11  Sea  MARC  I  side  looking  sonar  image  of  the  lava  sheet  flow  (a).  The  swath 
covered  by  Sea  Beam  is  delimited  by  the  dotted  line.  On  this  image  acoustic  shadows  are  white 
and  intense  backscatter  is  black.  Fissures  and  flow  channels  observable  in  this  image  are  out¬ 
lined  in  (b)  where  distance  along-track  originates  at  the  same  location  as  the  0  km  point  of  the 
acoustic  data  in  Figure  4.10,  The  arrow  points  to  the  location  of  the  bottom  photograph  (c). 
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IV.S.l.c  Rise  crest  environment 

The  acoustic  data  presented  in  Figure  4.10  corresponds  to  the  section  of 
bathymetry  delimited  by  the  arrows.  In  those  bounds,  the  sea  floor  is  nearly  flat,  with 
an  average  depth  of  2560  m.  With  a  6  second  ping  rate  and  a  ship  speed  of  about  1 
m/s,  the  along-track  sampling  interval  is  6  m  (compared  to  a  footprint  120  m  in 
diameter),  and  the  averaging  interval  used  for  low-pass  filtering  of  the  data  is  20 
pings. 

The  same  smoothing  effect  of  large  averaging  intervals  noted  in  the  San 
Clemente  basin  data  is  seen  here.  However,  peak  amplitude  variations  are  relatively 
higher  indicating  a  greater  variability  in  the  raw  data.  Small  trend  changes  over 
along-track  distances  1  km  or  less  are  noticeable  in  the  peak  amplitude  as  well  as  in 
the  total  energy  data.  They  can  also  be  followed  in  the  stacked  profiles  of  total 

i 

energy,  which  show  no  along-track  symmetry  about  vertical  incidence.  In  this  case, 
relation  of  the  acoustic  data  to  the  geological  processes  is  facilitated  by  the  unique¬ 
ness  of  the  area  and  the  availability  of  independent  measurements  made  simultane¬ 
ously  with  the  Sea  MARC  I  system. 

This  region  of  the  East  Pacific  Rise  around  lCTOS'N  is  characterized  by 
nearly  flat-lying  basalt  sheet  flows  covering  an  area  over  10  km  long  and  3  km  wide  on 
the  rise  axis  [Rise  Axis  Tectonic  Team,  19831  A  5  km  segment  from  this  area  imaged 
with  the  Sea  MARC  I  side-looking  sonar  system  is  shown  in  Figure  4.11a.  On  this 
image,  acoustic  shadows  are  white  and  intense  backscatter  is  black.  A  line  drawing 
(Fig  4.11b)  of  this  image  indicates  fissures  and  flow  channels  observable  in  the  sidescan 
data  as  well  as  in  the  area  covered  by  the  Sea  Beam  acoustic  data  considered  here 
(dotted  line).  The  basalts  are  fresh  as  evidenced  by  their  glass  coatings  and  the  pau¬ 
city  of  sediment  seen  in  bottom  photographs  in  the  area  (Fig.  4.11c).  Patches  of 
benthic  organisms  seen  in  the  bottom  photographs  and  temperature  measurements 
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also  indicate  that  this  region  of  the  rise  axis  is  hydrothermally  active  [Rise  Axis  Tec¬ 
tonics  Team,  1983]. 

The  distribution  of  flow  channels  outlined  inside  the  dotted  line  (Fig.  4.11b) 
coincides  remarkably  well  with  the  pattern  seen  in  the  total-energy  profiles  (Fig.  4.10) 
at  ±  13  degrees  incidence  (~  590  m  across  track  from  the  point  of  vertical  incidence). 
Total  energy  is  high  between  flow  channel  zones  and  drops  upon  crossing  the  zones. 
These  flow  channels  appear  in  the  side-scan  image  (Fig.  4.11a)  as  a  lighter  shade  of 
grey  than  their  surroundings,  indicating  reduced  backscattering  properties.  Because 
this  portion  of  seafloor  is  uniform  in  type  (basalts)  and  nearly  flat,  the  change  in  back- 
scattering  properties  is  most  likely  due  to  a  change  in  the  microroughness.  We  there¬ 
fore  conclude  that  in  this  area  the  acoustic  backscatter  is  dominated  by  bottom 
microroughness. 

IV*5.2  Angular  dependence 

Having  isolated  acoustically  homogeneous  areas  of  the  seafloor  by  defining 
acoustic  boundaries,  one  would  like  to  use  the  acoustic  data  to  identify  the  type  of 
bottom  within  each  area.  Because  Sea  Beam  measures  acoustic  backscatter  simul¬ 
taneously  at  sixteen  angles  of  incidence  spaced  2  2/3  degrees  apart  (usually  over  ±  20 
degrees  from  vertical  when  the  ship  is  not  rolling),  a  discrete  angular-dependence  func¬ 
tion  is  readily  obtainable.  This  function  is  a  potential  criterion  for  differentiating 
between  bottom  types. 

As  mentioned  previously,  it  has  not  been  possible  to  obtain  the  complete 
angular  dependence  from  the  data  presented  here  because  of  saturation  and  sidelobe 
interference.  As  a  result,  only  the  tails  of  the  function  are  shown  in  Figure  4.12  for 
the  three  types  of  sea  floor  considered.  In  this  figure,  the  levels  indicated  correspond 
to  the  relative  total  energy  measured  at  each  angle  over  flat,  acoustically  homogene¬ 
ous  regions  and  corrected  for  transmission  loss. 
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Figure  4.12  Angular  dependence  of  total  backscattered  energy.  The  tails  of  the  angular  depen¬ 
dence  function  shown  here  for  the  three  types  of  seafloor  considered  are  remarkable  for  the 
differences  in  their  energy  level.  However,  their  shapes  are  too  similar  to  tell  them  apart  on 
that  criterion  alone.  Vertical  bars  represent  1  standard  deviation  about  the  mean  value. 
The  label  " nodule/ mud”  corresponds  to  dense  nodule  coverage  over  a  mud  substrate,  "sand 
and  mud"  indicates  hemipelagic  sediments  found  in  the  North  San  Clemente  basin,  and 
"  rock"  refers  to  basalts  from  the  lava  sheet  flow  on  the  East  Pacific  Rise  at  10  *  N. 
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The  tails  of  these  three  angular-dependence  functions  are  mostly  remarkable 
for  the  differences  in  their  relative  energy  levels.  Hemipelagic  sediments  (sand  and 
mud)  appear  to  be  about  10  dB  above  nodules  and  10  dB  below  basalts.  Although 
very  few  deep  sea  data  in  the  kilohertz  range  exist  in  the  literature,  these  relative  lev¬ 
els  are  in  general  agreement  with  comparable  measurements  in  coastal  locations  sum¬ 
marized  by  Urick  [1983].  In  spite  of  the  very  different  nature  of  the  three  areas  con¬ 
sidered,  the  shape  of  their  partial  angular-dependence  function  is  similar  (Fig.  4.12), 
and  therefore  cannot  be  used  to  identify  them.  Consequently,  in  our  data  the  overall 
difference  in  levels  between  the  bottom  types  is  the  only  tangible  criterion  available  on 
which  to  separate  them.  However,  we  are  missing  the  information  contained  in  the 
beams  near  normal  incidence,  which  when  expressed  as  the  ratio  of  the  specular  back- 
scattered  energy  to  that  of  the  adjacent  nonspecular  beams  could  prove  a  good  indica¬ 
tor  of  the  nature  of  the  bottom. 

To  confirm  this,  a  modification  of  the  MPL  Sea  Beam  acoustic  data  acquisi¬ 
tion  system  is  in  progress  to  record  both  amplitude  and  phase  of  the  backscattered 
signals.  Sidelobe  interference  can  then  be  removed  without  degrading  the  bottom 
return  to  produce  a  complete  angular  dependence  function.  Also,  the  availability  of 
phase  information  gives  access  to  coherent  scattering  measurements  as  described  in 
Section  IV. 3,  and  yields  a  direct  estimate  of  the  bottom  microroughness  (Eq.  6). 

IV. 5. 3  Estimates  of  Surface  Statistics 

In  the  present  data  set,  we  are  restricted  to  echo  envelope  statistics  for 
estimating  parameters  of  the  rough  surface  such  as  rms  roughness  and  correlation 
area.  For  small  to  moderate  relative  roughness  (4£2a2<d),  these  parameters  are 
directly  related  to  the  shape  of  the  pdf  of  normal-incidence  echo  envelopes  through  the 
ratio  q  of  coherent  to  incoherent  energy  in  the  echo  (Eq.  23  and  24).  Estimates  of  q 
are  obtained  by  fitting  a  Rician  distribution  (Eq.  25)  to  the  histogram  of  normal 
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NORMALIZED  AMPLITUDE 


Figure  4.13  Probability  density  functions  of  peak  amplitude  in  the  near-speeular  returns  for  the 
three  types  of  seafloor,  (a)  manganese  nodule  area,  (b)  sediments  in  the  North  San  Clemente 
basin  and  (c)  lava  sheet  flow  on  the  crest  of  the  East  Pacific  Rise.  Histograms  are  from  Sea 
Beam  data.  Superimposed  curves  represent  the  corresponding  **  best"  fit  of  a  Rician  pdf  with 
the  parameter  7  holding  the  values  indicated  in  each  case. 
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incidence  echo  peak  amplitude.  This  method  has  the  advantage  of  being  independent 
of  system  calibration.  Therefore  peak-amplitude  variations  recovered  from  sidelobe 
levels  (e.g.  Fig.  4.5)  can  be  used  to  produce  a  histogram.  Three  such  histograms  for 
the  three  areas  under  consideration  are  shown  in  Figure  4.13.  The  overlying  curves 
represent  the  "  best"  fit,  in  a  chi-square  goodness  of  fit  [Bendat  and  Piersol,  1971] 
sense,  of  a  Rician  distribution  to  the  corresponding  values  of  7.  The  fit  is  not  good  in 
all  three  cases  and  a  different  distribution  may  yield  a  better  fit,  however  the  Rician 
distribution  is  the  only  one  which  relates  to  physical  parameters.  It  is  used  here  for 
this  reason.  The  three  areas  yield  markedly  different  values  of  7,  indicating  differences 
in  their  roughness  structure.  We  can  estimate  the  rms  roughness  of  each  area  from 
the  bottom  photographs  and  calculate  an  approximate  correlation  area  by  using  Eq. 
23  and  the  values  of  7  given  in  Figure  4.13.  With  Sea  Beam  k  =2tt/\c^  51  and 
A <f>  =  Ax  =2  2/3  degrees  so  that  at  normal  incidence  a  =  fl~2.8  in  Eq.  23. 

In  the  manganese  nodule  area,  7=0,  indicating  an  infinitely  rough  surface 
according  to  Eq.  23.  In  fact,  the  rms  roughness  estimated  from  box  cores  and  bottom 
photographs  is  about  2  cm  [de  Moustier,  1985b],  and  4£2<t2  —  4.15  >  1  so  that  the 
assumptions  leading  to  Eq.  23  are  violated.  The  surface  can  therefore  be  considered 
rough  for  our  purposes.  By  comparison,  7  =  4  for  the  North  San  Clemente  basin  data 
and  7  =  18  for  the  basalt  sheet  flow,  indicating  relatively  smoother  surfaces.  With 
these  two  values  of  7,  Eq.  23  yields  <t2/1/2  =2.03  lO^m4  and  4.5  10^m4  respectively.  If  <r 
were  the  same  for  both  regions,  the  correlation  areas  (/1/2)  would  differ  by  a  factor  of 
5,  making  it  possible  to  tell  the  two  types  of  bottom  apart.  Using  the  roughness  (2 
cm)  in  the  nodule  area  as  a  photographic  reference  scale,  it  is  reasonable  to  assume  a 
<1  cm  for  the  two  other  regions.  This  yields  a  lower  bound  on  the  correlation  areas 

2  2 

of  .203m  in  the  San  Clemente  basin  and  .045  m  on  the  basalt  sheet  flow. 


Both  correlation  areas  support  the  assumptions  which  lead  to  Eq.  23.  In  the 


122 


San  Clemente  basin,  the  Gaussian  assumption  is  justified  because  the  surface 
microroughness  is  due  mostly  to  animal  activity  and  is  therefore  randomly  distributed 
with  many  irregularities  within  the  ensonified  area.  These  irregularities  are  most 
likely  isotropic  so  that  /1^/2^45cm.  In  the  rise-crest  data,  the  irregularities  are 
presumably  anisotropic.  By  analogy  with  current  generated  ripples,  roughness  in  the 
direction  perpendicular  to  the  flow  of  lava  probably  has  a  longer  correlation  length 
than  roughness  parallel  to  the  flow  which  is  typically  characterized  by  linear  wrinkles 
(Fig.  4.11c).  From  bottom  photographs  taken  in  the  general  area,  the  spacing  between 
wrinkles  appears  random,  giving  some  justification  for  the  Gaussian  assumption.  The 
limiting  values  required  to  satisfy  alf  +  01$  =1  in  Eq.  21  given  that  lxl2  =  .045m2 
are  l^&Ocm  and  /2~7.5cm.  Although  the  value  of  the  correlation  area  (/1/2)  is  based 
on  an  empirical  estimate  of  the  rms  roughness  <j)  the  values  obtained  for  / x  and  /2  are 
consistent  with  the  general  roughness  character  observed  in  the  area.  It  follows  that 
in  the  limit  4 &2<t2<CL  (rms  roughness  of  the  order  of  1  cm  or  less  at  12  kHz),  two 
seafloor  environments  with  the  same  rms  roughness  can  be  differentiated  by  their 
correlation  area. 

IV.6.  CONCLUSIONS 

Deep  seafloor  acoustic  backscatter  data  measured  with  a  Sea  Beam  system 
have  been  shown  to  hold  a  wealth  of  information  on  the  nature  of  the  seafloor  sur¬ 
veyed.  Acoustic  boundaries  are  mapped  over  flat  areas  by  following  trend  changes  in 
acoustic  data  (peak  amplitude  at  normal  incidence  and  total  energy  otherwise)  both 
along  and  across  track.  Clues  to  the  nature  of  the  bottom  are  found  in  the  overall 
energy  level  of  a  partial  angular-dependence  function  of  backscattering  as  well  as  in 
the  shape  of  the  pdf  of  normal-incidence  echo  envelopes,  which  is  related  to  the  degree 
of  coherence  in  the  backscattered  acoustic  field. 
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Data  from  three  geologically  different  environments  (a  manganese  nodule 
field,  a  hemipelagic  sedimentary  basin,  and  a  rise-crest  basalt  sheet  flow)  have  been 
analyzed  and  validated  with  independent  measurements  (bottom  photographs  and 
side-looking  sonar  data).  The  backscattering  process  seemed  to  be  more  sensitive  to 
bottom  type  in  the  manganese  nodule  area,  to  bottom  roughness  in  the  rise-crest  data, 
and  to  a  combination  of  both  in  the  sedimentary  basin.  Total  energy  in  the  partial 
angular-dependence  function  was  highest  for  basalts  on  the  rise  crest;  it  was  roughly 
10  dB  lower  for  hemipelagic  sediment  in  San  Clemente  basin  and  another  10  dB  lower 
for  manganese  nodules. 

The  shape  of  the  pdf  of  echo  envelopes  was  indicative  of  a  rough  surface  in 
the  manganese  nodule  area,  and  of  smoother  surfaces  in  the  two  other  environments. 
This  was  confirmed  by  bottom  photographs.  When  given  the  same  roughness,  the  sedi¬ 
mentary  basin  and  the  rise-crest  sheet  flow  were  found  to  have  quite  different  correla¬ 
tion  areas.  For  small  relative  roughness  (4fc2er2<CL)  the  correlation  area  is  therefore  a 
useful  parameter  to  differentiate  between  seafloor  types  having  the  same  roughness. 
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IV.8*  APPENDIX  :  Envelope  distributions 

For  surfaces  rough  compared  to  an  acoustic  wavelength,  the  phase  of  the 
backscattered  echo  is  to  a  good  approximation  distributed  uniformly  over  the  interval 
(0,2tt),  and  the  set  of  pressures  received  can  be  assumed  to  be  normally  distributed 
with  variance  cr2.  From  the  central  limit  theorem,  this  approximation  is  better  as  the 
set  gets  larger  (^30).  Under  such  conditions,  the  amplitude  E  of  the  echo  as  defined  in 
Eq.  (1)  in  the  text,  has  been  shown  [Ol’shevskii,  1967;  Proud  et  ah,  1960]  to  be  Ray¬ 
leigh  distributed  with  a  pdf: 


w{E)  =  4  exp(-fr)  (A1) 

CP  2°p 

The  various  moments  of  E  are  by  definition: 

+oo 

<En  >=  J  En  W{E)dE  (A2) 

-00 

from  which,  after  substitution  of  Eq.  (Al),  we  obtain  the  second  moment: 

<E2>  =  2al  (A3) 

The  Rayleigh  pdf  can  then  be  expressed  in  terms  of  the  amplitude  of  the  echo  envelope 

E  and  its  mean  square  value  <E 2  >two  quantities  readily  measurable  in  our  data. 


W{E)  =  exp(  - 

<E2>  V  <E2> 


(A4) 


The  generalized  Rayleigh  distribution  or  Rice  distribution  [Rice,  1945]  has  been  used 
also  to  describe  the  effects  of  coherent  scattering  on  the  pdf  of  the  echo  envelope 
[Obshevskii,  1967;  Stanton,  1984].  In  this  case,  the  backscattered  pressure  wave  is 
considered  as  the  sum  of  a  coherent  component  of  rms  value  A  and  an  incoherent  com¬ 
ponent  normally  distributed  with  variance  c2,  respectively  analogous  to  the  sine  wave 
and  the  narrow-band  Gaussian  noise  of  the  original  Rician  distribution. 
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The  corresponding  echo  envelope  E  is  distributed  with  a  pdf: 


W(E)  = 


E  ,  E+A 2 .  . 


0  otherwise 

where  /„  is  the  zeroth  order  modified  Bessel  function 
moments  of  E  are  given  by 


E  >  0 


(A5) 


Remembering  Eq.  (A2),  the 


<£»>=( 2,2)2  r(-2-+l) 


(A6) 


For  even  order  moments,  the  confluent  hypergeometric  function  1Fl  (-n;  l,nr)  is  equal 
to  the  ntk  Laguerre  polynomial  [Helstrom,  1968’ 


Ln(r)  =  E 


(“O' 


(«!)2(n— t)! 


Therefore  the  second  moment  of  E  is 


(A7) 


<E2>  =  2<r  2  +  A2  (A8) 

To  express  the  Rician  distribution  in  terms  of  readily  measurable  quantities,  it  is  con¬ 
venient  to  define  a  parameter 


7  = 


which  is  the  coherent  to  incoherent  power  ratio  in  the  echo, 
and  (A9)  into  (A5)  yields: 


(A9) 

Substitution  of  Eq.  (A8) 


VF(£)=(24q)-f-  exp 

<£r> 


{2±2  )£2  +7<-E2> 

2  <E2> 


h 


E 

<E2>* 


(7(2-py))^ 


(A10) 


when  7<C1,  Eq.  (A10)  tends  to  the  Rayleigh  distribution 
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Eq.  (A4)  conversely,  when  q^>l  Eq.  (A10)  tends  to  a  Gaussian-like  distribution: 


W{E)  ~ 
with  X  = 


IX 

Vi 

7(*-l)2' 

2  x<E2> 

exp 

2 

E 


<E2>* 

since  /„(Z)  tends  to  (2irZ)~li  ez  when  Z  tends  to  infinity. 


(All) 
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Eq.  (A4)  conversely,  when  q»l  Eq.  (A10)  tends  to  a  Gaussian-like  distribution: 


W(E)  ~ 
with  X  = 


’  XI X 

Vi 

7(X-1)2] 

2  n<E2> 

exp 

2 

E 


<E2>* 
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A  SEA  BEAM  ACOUSTIC  DATA  ACQUISITION  SYSTEM 

A.1  INTRODUCTION 

In  the  last  quarter  of  1981,  a  Sea  Beam  bathymetric  survey  system  was 
installed  aboard  the  R/V  Thomas  Washington  of  the  Scripps  Institution  of  Oceanogra¬ 
phy  (SIO).  The  system  is  a  multibeam  echo-sounder  manufactured  by  General  Instru¬ 
ment  Corporation  and  its  main  function  is  to  generate  contour  maps  of  the  ocean  floor 
in  near-real  time. 

It  is  believed  that  more  information  is  contained  in  the  acoustic  signals 
returned  from  the  ocean  floor  than  is  necessary  for  contour  mapping.  However,  these 
signals  are  discarded  by  the  system  once  depth  determination  has  been  achieved.  In 
order  to  take  greater  advantage  of  Sea  Beam’s  capabilities  SIO’s  Marine  Physical 
Laboratory  devised  a  parallel  acoustic  data  acquisition  system  which  saves  the  Sea 
Beam  detected  echoes  before  they  reach  the  processing  stage.  This  paper  describes 
briefly  the  Sea  Beam  system  and  presents  the  method  used  to  record  the  echo 
envelopes  it  detects. 

A. 2  THE  SEA  BEAM  SYSTEM 

The  Sea  Beam  bathymetric  survey  system  comprises  two  major  subsystems:  a 
narrow  beam  echo-sounder  (NBES),  and  an  echo-processor  (EP).  The  narrow  beam 
echo-sounder  uses  sixteen  2  2/3  degree  beams  to  measure  water  depth  across  the 
ship’s  track.  The  swath  thus  obtained  on  the  seafloor  has  a  width  roughly  equal  to 
three  quarters  of  the  water  depth.  The  echo-processor  uses  a  mini-computer  as  a  cen¬ 
tral  control  and  processing  element,  providing  for  automatic  control  of  the  system. 
Figure  A  1  is  a  block  diagram  describing  the  functional  interfaces  of  the  Sea  Beam. 
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of  the  MPL  data  acquisition  system  is  shown  for  reference. 
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The  sonar  signal  originates  in  the  transmitter  signal  generator  which  produces  two 
low  impedance  12.158  kHz  sinewaves  in  quadrature.  A  gated  key  pulse  7ms  wide  is 
produced  by  the  timing  unit  in  the  echo-processor  module.  A  Universal  Graphic 
Recorder  (UGR)  usually  determines  the  ping  period.  The  signal  is  fed  into  a  digital 
pitch  compensator  which  steers  the  transmitted  beam  pattern  by  varying  the  phase  of 
the  signal  fed  into  each  of  the  twenty  power  amplifiers  relative  to  a  pitch  angle  from 
the  vertical  reference  gyroscope.  This  insures  vertical  projection  of  the  sonar 
transmission  within  the  limits  of  ±10  degrees  of  pitch  with  a  quarter  of  a  degree  of 
accuracy. 

Each  power  amplifier  drives  one  of  the  twenty  transducers  in  the  projector 
array.  This  array  is  mounted  along  the  ship’s  keel  and  is  roughly  2.8  m  long  by  0  16 
m  wide.  Figure  A. 2a  shows  that  the  corresponding  theoretical  transmit  beam  pattern 
is  54  degrees  wide  athwartship  and  2  2/3  degrees  wide  in  the  fore-aft  direction  The 
sidelobes  are  attenuated  30  dB  below  the  main  lobe  by  adjusting  the  power  output  of 
the  amplifiers.  The  receiving  unit  is  a  V-shaped  40  element  line  hydrophone  array 
centered  on  the  ship’s  keel.  It  measures  about  2.8  m  athwartships  and  0.4  m  fore-aft. 
The  corresponding  theoretical  beam  pattern  is  2  2/3  degrees  wide  athwartships  by 
20  degrees  in  the  fore-aft  direction  (Fig.  A.2b).  The  preamplifiers  associated  with 
each  hydrophone  create  four  quadrature  outputs  (45,  135,  225  and  315  degrees).  The 
resulting  160  signals  are  fed  into  a  resistors  matrix  which  steers  the  2  2/3  by 
20  degree  beam  in  the  athwartships  plane  to  form  sixteen  beams  within  +  20  and 
-  20  degrees  of  the  ship’s  vertical  axis.  The  sidelobes  are  attenuated  30  dB  below  the 
mainlobe  by  adjusting  the  signal  amplitudes  in  the  resistors  matrix. 

The  combination  of  the  projector  and  receiver  array  beam  patterns  outlines 
on  the  seafloor  sixteen  squares  2  2/3  degrees  on  a  side,  in  a  vertical  plane  through  the 
ship’s  fore-aft  axis,  but  roll  dependent  athwartships  (Fig.  A. 2c). 
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Figure  A.2  Sea  beam  transmit/ receive  geometry .  Computed  beam  pattern  cross-sections  in  the  athwartships  vertical  plane  cen¬ 
tered  on  the  array  and  in  the  vertical  plane  passing  through  the  ship’s  fore-aft  axis  are  shown  for  (a)  the  projector  array,  (b) 
the  receiver  array.  The  effect  of  Dolph-Chebyshev  amplitude  shading  is  also  illustrated,  (c)  is  a  summary  cartoon  showing  the 
angular  relationship  between  the  mainlobe  of  the  transmitted  beam  pattern  and  those  of  the  sixteen  preformed  beams. 


133 


Beams  are  numbered  from  the  keel  out:  1  to  8  for  both  port  and  starboard.  The  six¬ 
teen  beams  are  sent  on  one  hand  to  the  echo-processor  receivers,  via  the  beam  line 
drivers  which  act  as  buffers  between  the  sonar  and  the  receivers;  and  on  the  other 
hand  to  the  NBES  roll-compensator. 

In  response  to  a  roll  signal  from  the  vertical  reference  gyroscope,  within  the 
limits  of  ±  20  degrees  and  with  an  accuracy  of  a  quarter  of  a  degree,  sixteen  beams 
are  compensated  for  roll.  As  an  example  of  the  roll  effect,  a  20  degree  roll  to  port 
results  in  port  beam  number  8  looking  directly  below  the  ship.  This  beam  will  be  the 
vertical  beam,  and  no  echoes  will  show  on  port  side. 

The  roll-stabilized  echo  signals  go  to  the  NBES  receivers.  These  are  three 
identical  receivers  (port,  center,  starboard)  whose  function  is  to  detect,  filter,  amplify 
and  compress  the  dynamic  range  of  the  echo,  while  maximizing  the  signal  to  noise 
ratio  with  band-pass  filters. 

The  control  unit  of  the  echo-processor  sets  the  gain  of  these  receivers,  and 
gates  them  off  (rec.  blank,  in  Figure  A.l)  during  a  sonar  ping  and  its  reverberation.  It 
is  also  possible  to  output  a  depth  profile  on  the  UGR  by  selecting  one  or  a  combina¬ 
tion  of  these  receivers. 

In  the  echo-processor,  the  signals  (not  roll-compensated)  coming  from  the 
sonar  line  drivers  are  fed  into  sixteen  identical  receivers.  Each  receiver  has  a  narrow 
band  amplifier  centered  at  12.158  kHz,  a  detector  and  a  DC  amplifier.  The  echo- 
processor  receivers  output  linear  detected  signals  which  are  sent  to  the  analog-to- 
digital  (A/D)  converter  in  the  computer. 

The  computer  also  receives  roll  and  heading  inputs  through  a  synchro-to- 
digital  (S/D)  converter  and  performs  the  various  corrections  such  as:  ray  bending, 
receiver  gain  correction,  roll-compensation,  slant-range  calculation,  etc.  The  com¬ 
puter  does  the  depth  determination  and  keeps  track  of  the  bottom  by  setting  up 
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self-adjusting  gates  above  and  below  the  instantaneous  cross-track  bottom  profile. 
The  final  product  is  a  set  of  depths  and  horizontal  distances  across  the  ship’s  track 
which  are  subsequently  used  to  generate  contour  maps.  More  extensive  descriptions  of 
the  Sea  Beam  system  and  its  echo-processing  methods  may  be  found  in  Renard  and 
Allenou  (1979),  in  Farr  (1980)  and  in  de  Moustier  and  Kleinrock  (1985). 

A.3  A  Parallel  Acoustic  Data  Acquisition  System 

Because  Sea  Beam  does  not  preserve  the  bottom  backscattered  acoustic  sig¬ 
nals  it  receives,  and  because  these  signals  contain  information  about  the  nature  of  the 
seafloor,  it  was  necessary  to  build  a  parallel  acoustic  data  acquisition  system  to  record 
them  for  later  analysis. 

The  approach  taken  consists  in  recording  the  detected  echo  envelopes  of  Sea 
Beam’s  sixteen  beams  at  the  output  of  the  EP  receivers  (Fig.  A.l).  The  data  acquisi¬ 
tion  system  is  build  around  a  DEC  LSI  11/23  minicomputer  with  at  least  64  kw  of 
memory.  The  various  components  of  the  system  and  their  interface  with  Sea  Beam 
are  illustrated  in  the  block  diagram  of  Figure  A.3. 

The  analog  detected  signal  envelopes  are  tapped  at  the  output  of  the  EP 
receivers  along  with  a  time-varying  gain  (TVG)  which  is  applied  to  the  receivers  to 
compensate  for  acoustic  transmission  loss  through  the  water  column.  The  seventeen 
channels  are  sent  via  buffer  amplifiers  (Fig.  A. 4)  to  a  multichannel  12  bit  A/D  con¬ 
verter  with  direct  memory  access  (DMA).  The  buffer  amplifiers’  function  is  to  minim¬ 
ize  interferences  with  the  Sea  Beam  system.  In  addition,  the  ship’s  instantaneous  roll 
angle  is  sampled  from  the  vertical  reference  gyroscope  by  a  S/D  converter.  The 
resulting  eighteen  digitized  channels  (16  beams,  TVG  and  roll)  are  accumulated  in  a 
memory  buffer  until  the  bottom  echo  reception  cycle  is  completed.  When  the  recep¬ 
tion  cycle  is  completed,  the  buffer  is  written  on  a  standard  9-track  magnetic  tape  in  a 
DMA  operation  Meanwhile,  the  digitized  data  from  one  of  the  eighteen  channels  is 
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ECHO  PROCESSOR 


Figure  A.3  Block  diagram  of  the  MPL  Sea  Beam  acoustic  data  acquisition  system. 
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transferred  by  DMA  into  a  digital-to-analog  converter  which  drives  an  xy  display 
oscilloscope.  Outputting  a  selected  channel  on  the  oscilloscope  allows  the  operator  to 
monitor  the  proper  functioning  of  the  data  acquisition  system  and  the  performance  of 
its  bottom  tracking  algorithm.  Should  an  external  sound  source  interfere  with  Sea 
Beam's  reception  cycle,  the  interference  would  be  observable  in  the  display  so  that  the 
operator  would  know  when  corrective  action  (e.g.  phase  the  interfering  sound  source 
out  of  Sea  Beam’s  reception  cycle)  is  required. 

In  Figure  A. 5  a  simplified  timing  diagram  indicates  the  sequence  of  events 
taking  place  during  data  acquisition  The  timing  reference  for  this  sequence  is  the 
sonar  key  pulse  which  is  taken  from  the  echo-processor  in  Sea  Beam.  Upon  a  sonar 
key  pulse,  the  timer  of  an  accurate  (0.01%)  real-time  clock  is  loaded  with  a  delay 
corresponding  to  the  time  required  for  sound  to  travel  from  the  ship  to  the  bottom 
and  back.  When  the  timer  has  counted  down  to  zero,  it  sends  an  interrupt  to  the  LSI 
11/23  central  processing  unit  (CPU)  which  triggers  the  digitizing  sequence.  The  digit¬ 
ization  rate  is  about  10  kHz  for  eighteen  channels  so  that  for  any  one  channel  samples 
are  approximately  1.8  ms  apart.  The  CPLT  is  interrupted  each  A/D  conversion  scan 
cycle  to  take  in  the  ship’s  roll  data  from  the  S/D  converter. 

At  the  end  of  the  digitizing  sequence,  bottom  tracking  is  performed  by 
detecting  the  first  digitized  sample  which  exceeds  a  preset  threshold  in  the  set  of  six¬ 
teen  beams.  The  position  of  this  sample  inside  the  digitization  window  is  compared 
with  the  ideal  position  for  which  all  of  the  returns  fall  within  the  window.  Discrepan¬ 
cies  between  these  two  positions  are  resolved  by  adjusting  the  delay  time  between  the 
sonar  key  pulse  and  the  onset  of  digitization. 

The  digitized  data  is  then  written  to  tape  while  one  of  the  channels  is  also 
loaded  into  the  xy  display.  Both  operations  are  done  DMA.  The  CPU  spends  the 
remaining  time  waiting  for  the  next  sonar  key  pulse  and  the  above  sequence  repeats. 
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Figure  A.4  Buffer  amplifier  schematic. 
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Figure  A.5  Data  acquisition  timing  diagram. 
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Double  buffering  in  memory  is  used  to  avoid  overwriting  buffers  during  a  transfer  to 
tape.  By  this  method  one  buffer  fills  up  with  data  during  digitization  while  the  other 
is  being  written  on  tape. 

The  data  are  written  on  tape  sequentially.  Each  tape  record  corresponds  to 
one  Sea  Beam  reception  cycle  and  contains  14408  words  (8  word  header  +800  words 
x  18  channels).  The  record  header  is  formated  as  follows: 


word  0: 

number  of  words  in  buffer 

word  1: 

record  number  (sequential  from  start  of  tape) 

word  2: 

tape  error  status  for  previous  record 

word  3: 

digitization  rate  (bit  15-12),  conversion  gain  (bit  11-9), 

first  A/D  multiplexer  channel  to  be  digitized. 

word  4: 

delay  from  key  pulse  to  start  of  digitizing  in  increments  of  10  ms 

word  5: 

seconds  of  start  of  digitization 

word  6: 

hour-minute  of  start  of  digitization 

word  7: 

date  (day-month-year). 

Word  5  through  7  are  read  into  memory  from  a  time-of-day 

clock  once  per  Sea  Beam  reception  cycle. 

word  8-23 

— 16  beams 

word  24 

=  TVG 

word  25 

=  roll 

word  26-42 

=  16  beams 

word  43 

=  TVG 

word  44 

=  roll 

_ 

for  800  iterations. 

With  this  format  approximately  2950  records  fit  on  a  2400  ft.  reel  of  9-track  standard 
magnetic  tape  recorded  at  3200  bits  per  inch  (BPI).  Recording  time  per  tape  varies 
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with  the  key  pulse  repetition  rate  which  depends  on  the  water  depth.  Provision  has 
been  made  in  the  data  acquisition  software  to  allow  a  certain  number  of  reception 
cycles  not  to  be  recorded  on  tape  if  desired.  The  tape  consumption  may  therefore 
vary  depending  on  the  number  of  reception  cycles  skipped. 

Initial  post-processing  operations 

The  Sea  Beam  acoustic  data  recorded  on  tape  are  referenced  to  the  ship’s 
vertical  whose  orientation  varies  with  the  roll  of  the  ship.  To  compare  data  from  suc¬ 
cessive  reception  cycles,  it  is  therefore  convenient  to  roll-compensate  the  recorded 
data  so  as  to  reference  them  to  the  true  vertical.  In  addition,  because  the  EP 
receivers  have  inherent  gain  differences  between  them,  beam  to  beam  comparisons 
require  correcting  the  amplitudes  of  the  recorded  data  by  establishing  a  gain  reference 
common  to  the  sixteen  beam  channels.  This  gain  correction  and  the  roll-compensation 
are  common  to  most  post-processing  applications  of  the  recorded  data.  For  this  rea¬ 
son  they  are  outlined  in  the  following. 

Every  half  hour  the  Sea  Beam  system  calibrates  the  relative  gains  of  the  EP 
receivers  by  injecting  a  signal  from  the  test  signal  generator  (Fig.  A.l).  This  is  not  an 
absolute  calibration  for  the  echo-processor  only  seeks  to  compensate  for  gain 
differences  among  the  individual  receivers.  The  test  signal  is  not  constant  through 
time  at  the  input  to  the  receivers  because  it  is  multiplied  by  the  inverse  of  the  TVG 
which  compensates  the  detected  echoes  for  transmission  loss  during  normal  echo¬ 
sounding  operation.  This  results  in  a  nearly  constant  signal  level  at  the  output  of  each 
individual  receiver.  For  this  reason  it  is  possible  to  approximate  the  relative  gain  on 
any  one  digitized  receiver  output  by  the  mean  amplitude  observed  over  the  digitizing 
time  window.  The  data  recorded  on  tape  has  been  digitized  by  a  12  bit  A/D  con¬ 
verter  scaled  in  the  range  0  to  10  volts.  An  arbitrary  level  of  2.5  volts  (1024  A/D 
units)  was  selected  as  the  reference  level  to  which  individual  receiver  gains  would  be 
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adjusted.  2.5  volts  corresponds  roughly  to  the  mean  amplitude  observed  across  all 
sixteen  beams  during  a  calibration  cycle.  For  each  digitized  beam  channel,  a  receiver 
gain  correction  coefficient  is  calculated  as  the  ratio  of  the  reference  level  (2.5  volts) 
to  the  mean  amplitude  observed  on  that  channel  during  a  calibration  cycle.  In  subse¬ 
quent  processing  of  the  recorded  acoustic  data,  the  digitized  samples  for  each  of  the 
sixteen  beams  are  multiplied  by  the  corresponding  gain  correction  coefficients  thereby 
giving  them  a  common  reference. 

Gain  homogeneity  across  all  sixteen  beam  channels  is  important  during  the 
roll-compensation  process  because  this  process  involves  interpolation  between  beams. 
Ideally,  before  roll-compensation  a  refraction  correction  should  be  applied  to  deter¬ 
mine  the  true  angles  of  arrival  of  the  received  beams  with  respect  to  the  ship’s  verti¬ 
cal.  Because  of  the  computational  advantage  of  constant  received  beam  spacing  and 
because  angular  variations  due  to  refraction  are  small  for  small  angles  of  incidence, 
the  refraction  correction  is  omitted  in  the  following. 

Figure  A. 6  illustrates  the  roll-compensation  process.  The  sixteen  received 
beams  are  fixed  with  reference  to  the  ship’s  vertical  axis,  labeled  "keel”  in  Figures 
A. 6,  and  the  angle  between  this  axis  and  the  true  vertical  is  the  ship’s  roll  angle. 
These  beams  are  placed  in  an  array  of  fixed  stabilized  beams  spaced  2  2/3  degrees 
apart,  one  of  which  points  vertically  down.  Each  stabilized  beam  is  created  by  linear 
interpolation  between  the  two  received  beams  closest  to  its  direction.  With  a  zero  roll 
angle,  such  a  scheme  creates  fifteen  stabilized  beams:  one  vertical  beam  and  seven  on 
either  side. 

For  computing  purposes,  received  beams  are  numbered  1  to  16  from  port  to 
starboard.  In  the  example  of  Figure  A.6,  the  roll-stabilized  array  contains  twenty 
three  beams  numbered  1  to  23  from  port  to  starboard,  and  beam  12  lies  on  the  verti¬ 
cal.  By  comparison  with  the  fifteen  stabilized  beams  mentioned  above, 
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Figure  A.B  Roll-compensation  diagram. 
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this  corresponds  to  four  additional  beams  on  either  side  of  vertical,  allowing  for  roll 
angles  of  4  x  2-2/3  =  10-2/3  degrees  to  port  or  starboard.  Greater  roll  angles  can 
be  accommodated  by  adding  more  beams  to  the  stabilized  array.  The  terms  entering 
in  the  roll-compensation  algorithm  are  defined  as  follows: 

0  —  nominal  beam  spacing  (2  2/3  degrees) 

roll  =  ship’s  roll  angle  measured  by  the  vertical  reference  gyroscope 

B(I),  I  =1,  16  =  amplitudes  of  received  beams 

SB(J),  J  =1,  2M+15  =  interpolated  amplitudes  of  stabilized  beams 

M  =  number  of  additional  beams  on  either  of  vertical  in  the 

stabilized  array  (M  =4  in  Fig.  A. 6) 

Referring  to  Figure  A. 6,  the  offset  between  a  given  received  beam  and  its  stabilized 
counterpart  is  expressed  as  an  integer  number  of  beam  spacings: 


N  =  integer 


roll  + 0/2 


/» 


plus  a  fractional  remainder: 


AO  = 


NO 


roll  +  0/2  j  -  A 

Depending  on  the  direction  of  roll,  received  beams  are  shifted  right  or  left  relative  to 
the  stabilized  array.  The  offset  in  the  stabilized  beams  indexing  is  determined  accord¬ 
ingly. 

A  starboard  roll  corresponds  to  an  index  offset  of  K  =  M  -  N  in  the  stablized 
array.  The  stabilized  beams  are  then  calculated  by  linear  interpolation: 


SB{K+I)  =  B{I) 


e  -  A0 


6 


+  B{I+ 1) 


A0 


,  1  =  1,  15 


SB(J)  =0,  7=1,  K)  K*  0 


SB(J)  =0,  7  =  A +16,  2M+15;  K  *  2A/-1 
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Likewise  for  a  port  roll,  the  index  offset  is  K  =  M  +  N  and  the  linear  interpolation 
yields: 


SB(K+I)  =  B(I) 


Ae 


+  B(I+1) 


0  -  A  9 


0 


I  =  1,  15 


SB{J)  =0,  J  =  1,  K;  K  t*  0 


SB{J)  =0,  J  =  K+IQ,  2M+15;  K  *  2M-1 
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